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Abstract 

We study theoretical and experimental facets of mass determination of new particle 
states. Assuming supersymmetry, we update the quark and lepton mass matri- 
ces at the grand unification scale accounting for threshold corrections enhanced by 
large ratios of the vacuum expectation value of the two supersymmetric Higgs fields 
v u /vd = tan/3. From the hypothesis that quark and lepton masses satisfy a classic 
set of relationships suggested in some Grand Unified Theories (GUTs), we predict 
tan (3 needs to be large, and the gluino's soft mass needs to have the opposite sign 
to the wino's soft mass. Existing tools to measure the phase of the gluino's mass 
at upcoming hadron colliders require model-independent, kinematic techniques to 
determine the masses of the new supersymmetric particle states. The mass deter- 
mination is made difficult because supersymmetry is likely to have a dark-matter 
particle which will be invisible to the detector, and because the reference frame 
and energy of the parton collisions are unknown at hadron colliders. We discuss 
the current techniques to determine the mass of invisible particles. We review the 
transverse mass kinematic variable Mt2 and the use of invariant-mass edges to find 
relationships between masses. Next, we introduce a new technique to add additional 
constraints between the masses of new particle states using Mt2 at different stages 
in a symmetric decay chain. These new relationships further constrain the mass 
differences between new particle states, but still leave the absolute mass weakly 
determined. Next, we introduce the constrained mass variables M2c,lb, M2c,ub, 
M^clBi Mzcub to provide event-by-event lower-bounds and upper-bounds to the 
mass scale given mass differences. We demonstrate mass scale determination in 
realistic case studies of supersymmetry models by fitting ideal distributions to sim- 
ulated data. We conclude that the techniques introduced in this thesis have precision 
and accuracy that rival or exceed the best known techniques for invisible-particle 
mass-determination at hadron colliders. 
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Chapter 1 

Introduction 



In the mid-seventeenth century, a group of Oxford natural philosophers, including Robert Boyle 
and Robert Hooke, argued for the inclusion of experiments in the natural philosopher's toolkit 
as a means of falsifying theories [3J. This essentially marks the beginning of modern physics, 
which is rooted in an interplay between creating theoretical models, developing experimental 
techniques, making observations, and falsifying theories^ This thesis is concerned with mass 
determination of new particle states in these first two stages: we present a new theoretical 
observation leading to predictions about particle masses and their relationships, and we develop 
new experimental analysis techniques to extract the masses of new states produced in hadron 
collider experiments. The remaining steps of the cycle will follow in the next several years: New 
high-energy observations will begin soon at the Large Hadron Collider (LHC)[^] and history will 
record which theories were falsified and which (if any) were chosen by nature. 

No accepted theory of fundamental particle physics claims to be 'true', rather theories claim 
to make successful predictions within a given domain of validity for experiments performed 
within a range of accuracy and precision. The Standard Model with an extension to include 
neutrino masses is the best current 'minimal' model. Its predictions agree with all the highest- 
energy laboratory-based observations (1.96 TeV at the Tevatron) to within the best attainable 
accuracy and precision (an integrated luminosity C of about J dtC ~ 4 fb _1 ) [7][^J The Standard 
Model's agreement with collider data requires agreement with quantum corrections [HI Ch 10]. 
The success of the Standard Model will soon be challenged. 

1 There are many philosophy of physics subtleties about how physics progress is made 5 or how theories 
are falsified [6] which we leave for future philosophers. Nevertheless, acknowledging the philosophical interplay 
between experiment and scientific theory is an important foundational concept for a physics thesis. 

2 This does not mean that the LHC will observe new particles, but that the LHC will perform experiments 
with unexplored energies and luminocities that will constrain theories. 

3 Grossly speaking in particle physics, the domain of validity is given in terms of the energy of the collision, 
and the precision of the experiment is dictated by the integrated luminosity J dtC. Multiplying the J dtC with 
the cross-section for a process gives the number of events of that process one expects to occur. The larger the 
integrated luminosity, the more sensitive one is potentially to processes with small cross sections. 
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At about the time this thesis is submitted, the Large Hadron Collider (LHC) will (hopefully) 
begin taking data at a substantially larger collision energy of 10 TeV and soon after at 14 TeV 
[9]. The high luminosity of this collider (10 34 cm -2 sec -1 ) [9j enables us to potentially measure 
processes with much greater precision (hopefully about J dtC = 300 fb" 1 after three years). The 
LHC holds promise to test many candidate theories of fundamental particle physics that make 
testable claims at these new energy and precision frontiers. 

The thesis is arranged in two parts. The first regards theoretical determination of masses 
of new particle states, and the second regards experimental determination of masses at hadron 
colliders. Each part begins with a review of important developments relevant to placing the 
new results of this thesis in context. The content is aimed at recording enough details to enable 
future beginning graduate students to follow the motivations and, with the help of the cited 
references, reproduce the work. A full pedagogical explanation of quantum field theory, gauge 
theories, supersymmetry, grand unified theories (GUTs) and renormalization would require sev- 
eral thousand pages reproducing easily accessed lessons found in current textbooks. Instead 
we present summaries of key features and refer the reader to the author's favored textbooks or 
review articles on the subjects for refining details. 

Part I of this thesis contains Chapters [2](3} Part II of this thesis contains Chapters |4](8} 

Chapter [2] outlines theoretical approaches to mass determination that motivate the prin- 
ciples behind our own work. Exact symmetry, broken symmetry, and fine-tuning of radiative 
corrections form the pillars of past successes in predicting mass properties of new particle states. 
By mass properties, we mean both the magnitude and any CP violating phase that cannot be 
absorbed into redefinition of the respective field operators. We highlight a few examples of the 
powers of each pillar through the discoveries of the past century. We continue with a discussion 
of what the future likely holds for predicting masses of new particle states. The observation of 
dark-matter suggests that nature has a stable, massive, neutral particle that has not yet been 
discovered. What is the underlying theoretical origin of this dark matter? The answer on which 
we focus is supersymmetry. SUSY, short for supersymmetry, relates the couplings of current 
fermions and bosons to the couplings of new bosons and fermions known as superpartners. The 
masses of the new SUSY particles reflect the origin of supersymmetry breaking. Many theories 
provide plausible origins of supersymmetry breaking but are beyond the scope of this thesis. We 
review selected elements of SUSY related to why one might believe it has something to do with 
nature and to key elements needed for the predictions in the following chapter. 
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Chapter [3] which is drawn from work published by the author and his supervisor in |10] , 
presents a series of arguments using supersymmetry, unification, and precision measurements of 
low-energy observables to suggest some mass relationships and phases of the yet undiscovered 
superpartner masses^} In the early 1980s, Georgi and Jarlskog [11] observed that the masses 
of the quarks and leptons satisfied several properties at the grand unification scale. Because 
these masses cover six orders of magnitude, these relationships are very surprising. We discover 
that with updated experimental observations, the Georgi Jarlskog mass relationships can now 
only be satisfied for a specific class of quantum corrections enhanced by the ratio of the vacuum 
expectation value of the two supersymmetric Higgs v u /vd = tan (3. We predict that tan (5 be 
large (> 20) and the gluino mass has the sign opposite to the sign of the wino mass. 

Chapter [4] reviews the existing toolbox of experimental techniques to determine masses and 
phases of masses. If the model is known, then determining the masses and phases can be done 
by fitting a complete simulation of the model to the LHC observations. However, to discover 
properties of an unknown model, one would like to study mass determination tools that are 
mostly model-independentj^} Astrophysical dark-matter observations suggests the lightest new 
particle states be stable and leave the detector unnoticed. The dark-matter's stability suggests 
the new class of particles be produced in pairs. The pair-produced final-state invisible particles, 
if produced, would lead to large missing transverse momentum. Properties of hadron colliders 
make the task of mass determination of the dark-matter particles more difficult because we do 
not know the rest frame or energy of the initial parton collision. We review techniques based on 
combining kinematic edges, and others based on assuming enough on-shell states in one or two 
events such that the masses can be numerically solved. We continue with the transverse mass 
Mt which forms a lower-bound on the mass of the new particle which decays into a particle with 
missing transverse momentum, and was used to measure the mass of the W. The transverse 
mass Mt has a symmetry under boosts along the beam line. Next we review Mt2 [12] [13] which 
generalizes this to the case where there are two particles produced where each of them decays 
to invisible particle states. 

Chapter [5] begins the original contributions of this thesis towards more accurate and robust 
mass determination techniques in the presence of invisible new particle states. We begin by 

4 The phase of a Majorana fermion mass contributes to the degree to which a particle and its antiparticle 
have distinguishable interactions. The conventions is to remove the phase from the mass by redefining the fields 
such as to make the mass positive. This redefinition sometimes transfers the CP violating phase to the particle's 
interactions. 

5 By model-independent, we mean techniques that do not rely on the cross section's magnitude and apply 
generally to broad classes of models. 
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introducing novel ways to use the kinematic variable Mt2 to discover new constraints among 
the masses of new states. We also discuss the relationship of Mt2 to a recent new kinematic 
variable called Mqt- This work was published by the author in Ref. [14] . 

Chapter [6] introduces a new kinematic variable M2c(M_). Most model-independent mass 
determination approaches succeed in constraining the mass difference M_ = My — M/v between 
a new states Y and the dark-matter particle N, but leave the mass scale M+ = My + 
more poorly determined. M2c(M_) assume the mass difference, and then provides an event- 
by-event lower bound on the mass scale. The end-point of the M2C distribution gives the mass 
My which is equivalent to the mass scale if M_ is known. In this chapter we also discover a 
symmetry of the M<ic distribution which for direct Y pair production, makes the shape of the 
distribution entirely independent of the unknown collision energy or rest frame. Fitting the 
shape of the distribution improves our accuracy and precision considerably. We perform some 
initial estimates of the performance with several simplifying assumptions and find that with 250 
signal events we are able to determine M/v with a precision and accuracy of ±6 GeV for models 
with My + Mj\r/(My — Mtv) ~ 3. This chapters results were published in Ref |15j . 

Chapter [7j which is based on work published by the author with Alan Barr and Graham 
Ross in Ref. [16J, extends the M2C kinematic variable in two ways. First we discover that 
in the presence of large upstream transverse momentum (UTM), that we are able to bound 
the mass scale from above. This upper bound is referred to as M2c,ub- Here we perform a 
more realistic case study of the performance including backgrounds, combinatorics, detector 
acceptance, missing transverse momentum cuts and energy resolution. The case study uses 
data provided by Dr Alan Barr created with HERWIG Monte Carlo generator [T71 HH1 QI5] to 
simulate the experimental data. The author wrote Mathematica codes that predict the shape 
of the distributions built from properties one can measure with the detector. We find the mass 
Mtv by fitting to the lower-bound distribution M^c and the upper-bound distribution M2CUB 
distribution shapes. Our simulation indicates that with 700 events and all anticipated effects 
taken into consideration that we are able to measure the Mpf to ±4 GeV for models with 
My + M]\r/(My — M^v) ~ 3. This indicates that the method described is as good as, if not 
better than, the other known kinematic mass determination techniques. 

Chapter [8] which is based on work published by the author with Alan Barr and Alex Pinder 
in Ref. [20] , extends the constrained mass variables to include two mass differences in the variable 
M3C. We discuss the properties of the M%c distribution. We observe that although the technique 
is more sensitive to energy resolution errors, we are still able to determine both the mass scale and 
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the two mass differences as good if not better than other known kinematic mass determination 
techniques. For SPS la we forecast determining the LSP mass with ±2.4 GeV with about 3600 
events. 

Chapter [9] concludes the thesis. We predict the wino's mass will have the opposite sign of 
the gluino's mass. We develop new techniques to measure the mass of dark-matter particles 
produced at the LHC. Our techniques work with only two or three new particle states, and 
have precision and accuracy as good or better than other known kinematic mass determination 
techniques. 

To facilitate identifying the original ideas contained in this thesis for which the author is 
responsible, we list them here explicitly along with the location in the thesis which elaborates 
on them: 

• Chapter [3] The updated values of the strong coupling constant, top quark mass, and 
strange quark mass lead to quantitative disagreement with the Georgi-Jarslkog mass re- 
lationships unless one uses tan/3 enhanced threshold corrections and fix the gluino's mass 
to be opposite sign of the wino's mass published in Ref [TO] . 

• Chapter [5] The relationship between Tovey's Mqt and Lester, Barr, and Summer's Mt2 
variables published in Ref [14] . 

• Chapter [5] Using Mt2 along intermediate stages of a symmetric cascade decay to discover 
new constraints among the masses published in Ref [14 . 

• Chapter [6] and [7] The definit ion and use of the variable M^c to determine mass of 
dark-matter-like new particle states published in Ref [151 and [16] . 

• Chapter [7] The ability to obtain an event-by-event upper bound on the mass scale given 
the mass difference when the mass difference is known and the event has large upstream 
transverse momentum published in Ref [16]. 

• Chapter [8] Complications and advantages to the logical extension variable M^c and their 
use as a distribution which is published in Ref [20j . 

• A set of Mathematica Monte Carlo simulations, and Mt2 calculators used to test the above 
contributions for M2C1 M%c and some C++ codes for M2C and M2CUB which has not yet 
been published. The author will be happy to share any codes related to this thesis upon 
request. 
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Chapter 2 



Mass Determination in the Standard 
Model and Supersymmetry 



Chapter Overview 

This chapter highlights three pillars of past successful new-particle state prediction and mass 
determination, and it shows how these pillars are precedents for the toolbox and concepts em- 
ployed in Chapter [3fs contributions. The three pillars on which rest most demonstrated success- 
ful theoretical new-particle-state predictions and mass determinations in particle physics are: 
(i) Symmetries, (ii) Broken symmetries, and (hi) Fine Tuning of Radiative Corrections. We use 
the narrative of these pillars to review and introduce the Standard Model and supersymmetry. 



The chapter is organized as follows. Section 2.1 gives a few historical examples of these 
pillars: the positron's mass from Lorentz symmetry and from fine tuning, the f2~'s mass from 
an explicitly broken flavor symmetry, the charm quark's mass prediction from fine tuning, and 
the W ± and Z° masses from broken SU{2)l x U(1)y gauge symmetries of the Standard Model. 
These historical examples give us confidence that these three pillars hold predictive power and 
that future employments of these techniques may predict new particle states and their masses. 



Section 2.2 introduces astrophysical observation of dark-matter which suggests that nature has 



a stable, massive, neutral particle that has not yet been discovered. Section 2.3 introduces 
key features of Supersymmetry, a model with promising potential to provide a dark-matter 
particle and simultaneously address many other issues in particle physics. We discuss reasons 
why Supersymmetry may describe nature, observational hints like the top-quark mass, anomaly 
in the muon's magnetic moment (g — 2)„ ; and gauge coupling unification. We also review the 
classic SU(5) Georgi-Jarlskog mass relations and tan /3-enhanced SUSY threshold corrections. 
Chapter [3] will discuss how updated low-energy data and the Georgi-Jarlskog mass relationships 
may provide a window on predicting mass relationships of the low-energy supersymmetry. 
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2.1 Past Mass Determination and Discovery of New Particle 
States 

2.1.1 Unbroken Lorentz Symmetry: Positron Mass 

Our first example uses Lorentz symmetry to predict the existence of and mass of the positron. 
Lorentz symmetry refers to the invariance of physics to the choice of inertial frame of reference. 
The Lorentz group, which transforms vectors from their expression in one inertial frame to 
another, is generated by the matrix M al3 . For a spin-1 object (single index 4- vector), the group 
elements, parameterized by a antisymmetric 4x4 matrix 6^ u , are given by 

A(0) = exp (iM^e^) (2.1) 

where the antisymmetric matrix of generators M^ u satisfies 

[M piy , M pa ] = i{gV p M ua - g vp M^ - g pu M vp + g ua ' M pp ) (2.2) 

and g pv = diag(l, —1,-1,-1) is the Lorentz metric. We emphasize that each entry of M^ v is 
an operator (generator), where each entry of 6^ u is a number. The generators of the Poincare 
symmetry are the generators of the Lorentz symmetry supplemented with generators for space- 
time translations P^ satisfying [P^,P U ] = and 

[M^, P p ] = -i(g w P v - gu P Pn). (2.3) 

Supersymmetry, introduced later in this thesis, is a generalization of the Poincare symmetry. 

In 1928 designing a theory with Lorentz symmetry was a major goal in early quantum 
mechanics. As Dirac described in his original paper, two problems persisted, (1) preservation of 
probability in quantum mechanics requires a first order linear equation for time evolution and 
(2) the presence of positive and negative energy solutions. The Klein-Gordon equation is 

(D^D P - m 2 )(p = (2.4) 

where = — ieA^. It is invariant under Lorentz transformations but suffers from both 
problems: the equation is second order in dt, and it has solutions proportional to exp(iu;t) and 
exp(— iuit). 

Dirac's 1928 paper on the Dirac equation [21] claims only to solve problem (1). Because 
problem (2) was not solved, Dirac claimed "The resulting theory is therefore still only an ap- 
proximation". However, the paper shows how to do Lorentz invariant quantum mechanics of 
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spin 1/2 fields. Although in different notation, Dirac's paper showed that if one has a set of four 
matrices 7 M that satisfy 

-f-f + 7 V = 2g^ (2.5) 

then the system of equations 

{i^Df, - m)ip = (2.6) 

transforms covariantly under Lorentz transformations if we define a new class of transforms 
-A-i/2 = exp(i0 /JI/ (M 1/ / 2 ) Ml/ ) where the group is generated by 

(M 1/2 r=J( 7 Y-7Y). (2-7) 



The new generators M 1 / 2 satisfy Eq(2.2) so they form a representation of the Lorentz group 
specific to spin 1/2. The Dirac field transforms as ip 1 = h\j2'4> and the 7^ matrices transform as 
A^ = A-/ 27 ^A 1/2 . 

Dirac interprets that the negative-energy solutions to the equations will behave as if they 
have the opposite charge but same mass in the presence of a background electromagnetic field. 
The formation of a relativistic quantum mechanics which still possesses negative energy solutions 
suggests that this alternative solution with negative energy and the opposite charge and same 
mass may be a prediction of relativistic quantum mechanics. Indeed, Anderson observed the 
negative energy version of the electron in 1933 (22 J Q The positron has the opposite charge of 
the electron but the same mass. 

Today the Dirac equation is interpreted in terms of quantum field theory where all energies 
are always considered positive. With hindsight we see that Dirac's motivation was partially 
wrong, and the Klein-Gordon equation provides just as good an argument for antiparticles. 

2.1.2 Renormalization, Fine Tuning and the Positron 

The next example is not the historical origin of the positron's prediction, but could have been 
if physicists in the early 1900s understood the world with today's effective field theory tools. 
The electron's self-energy poses another problem which effective quantum field theory and the 
existence of the positron solve [23] . 

If we divide up the electron's charge into two piecea^J we ask how much energy is stored in 
the system in the process of bringing the two halves together from infinity? At what distance 



1 Anderson's first paper makes no reference to Dirac. However, he does introduce the term positron and suggests 
renaming the electron the negatron. 

2 The choice of two pieces is arbitrary and simplifies the calculations to illuminate the fine-tuning concept 
trying to be communicated. 
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apart will the 'self energy' of the electron be equal to the mass energy of the electron? This 
is approximately the classical electron radius^J and the answer is around ten times bigger than 
atomic nuclei at around r e = 7 x 10 -14 m. Electron-electron scattering can probe closer than 
this for collisions with y^s ~ 2 MeV. Also electrons are emitted from nuclei during (3 decay 
suggesting the electron must be smaller than the nucleus. We now break the mass up into two 
quantities: the bare mass m eo and the 'self-energy' mass 5m with the observed mass equalling 
m e = m eo + 5m e . Phrasing self energy in terms of a cut-off one finds 

5m e « j (2.8) 

and the cut off A indicates the energy needed to probe the minimum distance within which two 
pieces of the electrons mass and charge must be contained. At the Plank scale, this requires a 
cancelation between the bare mass m e>0 and the self energy bm e to more than 22 decimal places 
to end up with the observed mass m e = 0.511 MeV! Fine tuning is where large cancelations 
between terms are needed to give an observable quantity. This large cancelation could be 
perceived as an indication of missing physics below the scale of about 2 MeV where we would 
have had a cancelation of the same size as the observable quantity m e . 

The effective quantum field theory used to describe electromagnetism introduces the positron 
with the same mass as the electron. The positron acts to moderate this "self-energy". In QFT, 
the self energy of the electron is partly due to electron interaction with electron, but also partly 
due to electron interaction with a virtual positron. This is because the two interaction vertices 
x and y are integrated over j d xd 4 y so the interaction The resulting self energy, in terms of a 
cut off A, is 

^--ln^ (2.9) 
m e 4vr mj, 

where a = e 2 /4-7r. Weisskopf and Furry were the first to discover that quantum field theory 
with both positrons and electrons leads to only a logarithmic divergence [24, 25J. Now we see 
that taking the cut off A to the Plank scale only gives a 6% correction. There is no longer a 
cancelation of two large terms, an issue solved by introducing new physics below the scale at 
which the low-energy effective theory became fine tuned. 

2.1.3 Broken Flavor Symmetry: VL~ Mass 

The next example brings us to the early 1960s when Murray Gell-mann [26] and Yuval Ne'eman 
|27j |^]were both studying broken SU(3) flavor symmetry as a way to understand the zoo of par- 
3 The classical electron radius is 4x this value. 

4 Yuval Ne'eman, like the author of this thesis, was a member of the military while studying for his PhD [28] . 



tides being discovered in the 1950s and 1960s. Their study led to the theoretical determination 
of mass of the £l~ baryon before it was observed. 

To place the SU(3) flavor symmetry in context, we begin first with the SU(2) isospin sym- 
metry. Isospin symmetry postulates that protons p and neutrons n are indistinguishable if one 
ignores electromagnetism. Therefore, the equations and forces governing p and n interactions 
should have a symmetry on rotating the proton's complex field p into the neutron's complex 
field n with an SU{2) rotation. The symmetry is explicitly broken by the electromagnetic force 
and to a lesser extent the quark masses, and a special direction is singled out. We can label the 
eigenvalues along this broken direction. The eigenvalues of ^3 isospin generator label the states 
along the broken isospin axis: I3 = 1/2 denotes a proton, and I3 = —1/2 denotes a neutron. 
Today we can trace the isospin symmetry to the presence of up and down quarks with nearly 
the same mass. Isospin symmetry is broken both by electromagnetism and by the up and down 
quark mass difference. 

Next we wish to understand how SU(3)f symmetry predicted the mass of the . The 
SU(3) flavor symmetry is an extension of isospin symmetry. It can be understood in modern 
language as the existence of three light quarks. The symmetry is broken because the strange 
quark is considerably more massive ( ~ 103 ± 12 GeV ) than the up and down quark ( m u ~ 
2.7 db 0.5 GeV, and = 5.7 ± 0.5 GeV) [^J The group SU(3) has isospin SU{2) as a supgroup 
so I3 remains a good quantum number; in addition SU(3) states are also specified by the 
hypercharge Y. Quarks and anti-quarks are given the following charge assignments (i3,Y): 
u = (1/2,1/3), d = (-1/2,1/3), s = (0,-2/3), u = (-1/2,-1/3), d = (1/2,-1/3), s = 
(0,2/3). Representations of SU(3) are formed by tensors of the fundamental its conjugate with 
either symmetrized or antisymmetrized indices and with traces between upper and lower indices 
removed. Representations are named by the number of independent components that tensor 
possesses and shown by bold numbers like 3, 8, 10, 27 etc. 

Gell-mann and Ne'eman did not know what representation of SU(3) was the correct one to 
describe the many baryons and mesons being discovered; to describe the spin 3/2 light baryons, 
they were each considering the 10 and the 27. The representation 10 is formed by B abc where 
a, b, and c are indexes that run over u, d, and s, and where the tensor B abc is symmetric on 



interchanges of the three indices. The states are displayed as the red dots in Fig. 2.1 where the 
conserved quantum numbers of I3 is plotted against hypercharge Y. The 27 is given by a tensor 
B ab cd where (ab) is symmetrized, (cd) is symmetrized, and the 9 traces B ac cd = are removed. 



5 Quark masses are very difficult to define because they are not free particles. Here we quote the current quark 
masses at an M S renormalization scheme scale of fi — 2 GeV as fit to observations in Chapter |3] 
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J p =(3/2) + Baryons 




Figure 2.1: The states of the SU(3) 10 of J p = 3/2+ baryons in 1962. Also shown is the SU(3) 
27. 



The 27 is shown by the smaller blue dots. The particles and the observed masses as of 1962 are 



shown in the Fig 2.1 



In July 1962, both Gell-mann and Ne'eman went to the 11 th International Conference on 
High-Energy Physics at CERN. Ne'eman learned from Sulamith and Gerson Goldhaber (a hus- 
band and wife pair) that K + (us,Is = +1/2, Y = 1) vs N (ddu,/3 = — 1/2, Y = 1) scattering 
did not lead to a resonance at (^3 = 0, Y = 2) with a mass range near 1000 MeV as one would 
expect if the pattern followed was that of the 27 |29j. Gell-mann had also learned of the Gold- 
habers' negative result, now known as the Goldhaber gap. During the discussion after the next 
talk at the conference, Gell-mann and Ne'eman both planned on announcing their observation; 
Gell-mann was called on first, and announced that one should discover the f2~ with a mass near 
1680 GeV. The was first seen in 1964 at Brookhaven [30 with a mass of Mq- = 1686 db 12 
MeV. Amazingly the spin of was not experimentally verified until 2006 [31 . 

The broken SU(3)f flavor symmetry successfully predicted and explained the masses of 
the Q~ and many other baryons and mesons. The general formula for the masses in terms of 
the broken symmetry was developed both by Gell-mann and by Okubo and is known as the 
Gell-mann- Okubo mass formula [32J. 

Despite this success (and many others) the quark-model theoretical description used here is 
still approximate; the strong force is only being parameterized. The quarks here classify the 
hadrons, but the masses of the quarks do not add up to the mass of the hadron. Quark masses 
that are thought of as being about 1/3 the mass of the baryon are called "constituent quark 
masses". These are not the quark masses that are most often referenced in particle physics. A 
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consistent definition of quark mass and explanation of how these quarks masses relate to hadron 
masses will wait for the discovery of color SU(3) C . In SU(3) C the masses of the quarks are 
defined in terms of Chiral Perturbation Theory (xPT) and are called "current" quark masses. 
We will be using current quark masses as a basis for the arguments in Chapter [3j 

Another important precedent set here is that of arranging the observed particles in a rep- 
resentation of a larger group that is broken. This is the idea behind both the Standard Model 
and the Grand Unified Theories to be discussed later. The particle content and the forces are 
arranged into representations of a group. When the group is broken, the particles distinguish 
themselves as those we more readily observe. 

2.1.4 Charm Quark Mass Prediction and Fine Tuning of Radiative Correc- 
tions 

In the non-relativistic quark model, the neutral K° is a bound state of d and s and K° is a 
bound state of s and d. The Hamilontian for this system is given by 

" (*• **) ( M J M>) (£) ' < 2 ' 10 > 

A non-zero coupling 5 between the two states leads to two distinct mass eigenstates. For 5 « 
M , the two mass eigenstates have almost equal mixtures of K° and K° and are called K\ and 
K2. Experimentally the mass splitting between the two eigenstates is Mk 2 ~ = 3.5 x 1CP 6 
eV. 

During the 1960s, a combination of chiral four-Fermi interactions with the approximate chiral 
symmetry SU(3)l x SU(3)r of the three lightest quarks was proving successful at describing 
many hadronic phenomena; however, it predicted 5 was non zero. Let's see why. The effective 
weak interaction Lagrangian that was successful in atomic decays and Kaon decays |33j [31] was 
given by 

C=^-J»4 where J" = Jfa + Jfa (2.11) 

and 

Jfa = ^i(l- 7 5 ) e + I 7^(l- 7 5 )^ (2.12) 
Jgyj = u 7 ^(l-7 5 )dcos0 c + n 7 ^(l-7 5 )ssin# c . (2.13) 

and 6 C is the Cabibbo angle. Using a spontaneously broken SU(3)l x SU(3)r |35 j one can 
calculate the loops connecting us to us. These loops are responsible for the Kaon mass splitting 
6 In this Hamilontian, we're neglecting the CP violating features. 
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and give 

Mk 2 ~M Ki = ^* C ° s2e 4 lt nec f 2 K M K G F (G F A 2 ) (2.14) 

where A is the cut off on the divergent loop and fx ~ 180 GeV [36] [37] [35] [33] . Using this 
relation, the cut-off cannot be above 2 GeV without exceeding the observed Kaon mass splitting. 
There are higher order contributions, each with higher powers of GpA 2 , that may cancel to give 
the correct answer. Indeed Glashow, Ilioupoulos, and Maiani (GIM) [37 j observe in a footnote 
"Of course there is no reason that one cannot exclude a priori the possibility of a cancelation in 
the sum of the relevant perturbation expansion in the limit A — > oo" . The need for an extremely 
low cut off was a problem of fine tuning and naturalness with respect to radiative corrections. 

This solution to the unnaturally low scale of A suggested by the mass splitting Mk 2 — Mk x 
holds foreshadowing for supersymmetry Q GIM proposed a new broken quark symmetry SU(A) 
which required the existence of the charm quark, and also discuss its manifestation in terms of 
a massive intermediate vector boson W + , W~ . They introduce a unitary matrix U, which will 
later be known as the CKM matrix. We place group the quarks into groups: d l = d, s, . . . and 
u % = u, c, . . .; then the matrix Uij links up-quark i to down-quark j to the charged current 

(J( H) r = «Y^(l-7 5 K^. (2-15) 

The coupling of the Kaon to the neutral current formed by the pair of W + with W~ . In the 
limit of exact <S£/(4) (all quark masses equal) the coupling of the Kaon to the neutral current 
is proportional to A oc (J2j =uc UdjUj s ) = 5d s = 0. The coupling 5 2 between K° and K° is 
proportional to A 2 . This means that in the limit of SU (4) quark symmetry, there would be no 
coupling to enable K\ and Ki mass splitting. 

However, the observed K\ K2 mass splitting is non-zero, and SU(4) is not an exact symmetry; 
it is broken by the quark masses. The mass splitting is dominated by the mass of the new quark 
m c . GIM placed a limit on m c S 3 GeV. One might think of the proposed SU(4) symmetry 
becoming approximately valid above scale A ~ m c . The new physics (in this case the charm 
quark) was therefore predicted and found to lie at the scale where fine tuning could be avoided. 

2.1.5 The Standard Model: Broken Gauge Symmetry: W ± , Z° Bosons Mass 

The Standard Model begins with a gauge symmetry SU (3) c x SU (2) l x U(1)y that gets sponta- 
neously broken to SU(3) C x U{1)em- For a more complete pedagogical introduction, the author 
has found useful Refs. [381 EH] SO] and the PDG [HI Ch 10]. 

7 I have learned from Alan Barr that Maiani also sees the GIM mechanism as a precedent in favor of super- 
symmetry. 
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Table 2.1: Standard Model field's transformation properties. * Indicates the field does not 
transform covariantly but rather transform as a connection^. 



The field content of the Standard Model with its modern extensions for our future use is 
given in Table |2.1| The i indexes the three generations. The fermions are all represented as 
2-component Weyl left-handed transforming states^} The gauge bosons do not transform covari- 
antly under the gauge group, rather they transform as connections^} In 1967 [12], Weinberg set 
the stage with a theory of leptons which consisted of the left-handed leptons which form SU{2) 
doublets L' 1 , the right-handed charged leptons which form SU{2) singlets (e c ) 4 , and the Higgs 
fields H which is an SU{2) doublet. Oscar Greenberg first proposed 3-color internal charges of 
SU(3) C in 1964 [43j^J It was not until Gross and Wilczek and Politzer discovered asymptotic 
freedom in the mid 1970s that SU(3) C was being taken seriously as a theory of the strong nu- 
clear force [H] [16] . The hypercharge in the SM is not the same hypercharge as in the previous 
subsection p*] The U(1)y hypercharge assignments Y are designed to satisfy Q = a 3 / 2 + Y/2 
where Q is the electric charge operator and a 3 /2 is with respect to the SU(2)i doublets (if it 
has a charge under SU{2)l otherwise it is 0). 

The gauge fields have the standard Lagrangian 



C 



W,B,G 



(2.16) 



where B 
0, 



ig'B,,. Likewise W, 



jlD^^Dy] when is the covariant derivative with respect to U(l): = 
[D^Du] with D^ = d^- igfW? and = j-[D„D v ] with 



D/j, = df! — igs^G^ with a/2 and A/2 the generators of SU{2) and SU(3) gauge symmetry 



The (2,0) is the projection of Eq(|2.7[) onto left handed states with P L = 1/2(1 — 7 ). In a Weyl basis, 



Eq(2.7l is block diagonal so one can use just the two components that survive Pl- If e c transforms as a (2,0) 
then ia 2 (e c )* transforms as (0,2). 

9 The gauge fields transform as connections under gauge transformations: WL = UW^U -1 + (i/ g)Ud / _ l U~ 1 . 
10 Greenberg, like the author, also has ties to the US Air Force. Greenberg served as a Lieutenant in the USAF 
from 1957 to 1959. A discussion of the history of SU(3) C can be found in Ref. |44] , 



In the SM, the d and s (both left-handed or right-handed) have the same hypercharge, but in flavor 5)7(3) 



/ 



(2.1.31 they have different hypercharge. 
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respectively. The \j g factors in these definitions give the conventional normalization of the field 
strength tensors. At this stage, SU{2) gauge symmetry prevents terms in the Lagrangian like 
|M^W M W which would give W a mass. 

The leptons and quarks acquire mass through the Yukawa sector given by 



C Y = Y$L\i(T 2 H*)(e c y + YfjQ^i^H^idy + Y%<?H(u c y + h.c. (2.17) 



where we have suppressed all but the flavor indices Because H and ia2H* both transform 



as a 2 (as opposed to a 2), these can be used to couple a single Higgs field to both up-like and 



down-like quarks and leptons. If neutrinos have a Dirac mass then Eq(2.17) will also have a 



term Y^L l H{v c y . If the right-handed neutrinos have a Majorana mass, there will also be a 
teimM Rjk (v c y(v c ) k . 

With these preliminaries, we can describe how the W boson's mass was predicted. The story 
relies on the Higgs mechanism that enables a theory with an exact gauge symmetry to give mass 
to gauge bosons in such a way that the gauge symmetry is preserved, although hidden. The 
Higgs sector Lagrangian is 

C H = D^D^H - fi 2 H^H - X(H^H) 2 (2.18) 

where the covariant derivative coupling to H is given by = 8^ — ig^W® — i^B^. The gauge 
symmetry is spontaneously broken if fj? < 0, in which case H develops a vacuum expectation 
value (VEV) which by choice of gauge can be chosen to be along (H) = (0, v/y2). The gauge 
boson's receive an effective mass due to the coupling between the vacuum state of H and the 
fluctuations of W and B: 

<£h) = ((WX + (WX) + y (9'B, + g(W%) 2 . (2.19) 

From this expression, one can deduce W'h and W 2 have a mass of Myy = g 2 v 2 /4, and the linear 
combination = (g 2 + g' 2 )~ x l 2 {g' B^ + g(W 3 )^) has a mass M| = (g 2 + g' 2 )v 2 /A. The massless 
photon An is given by the orthogonal combination An = (g 2 + g' 2 )~ 1 / 2 (gB fJ- — g'{W 3 )^). The 
weak mixing angle 0w is given by smOy/ = g'(g 2 + g' 2 ) -1 / 2 , and electric charge coupling e is 
given by e = g sin 9w ■ 

Before or Z° were observed, the value of v and sin 2 9\y could be extracted from the 
rate of neutral current weak processes, and through left right asymmetries in weak process 
experiments |47| . At tree-level, for momentum much less than M^y, the four- Fermi interaction 



12 Because both H and L transform as 2 their contraction to form an invariant is done like L°H b e a f, where t a \, 
is the antisymmetric tensor. 
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Figure 2.2: One Loop Contributions to the Higgs mass self energy, (a) Fermion loop, (b) Higgs 
loop. 

can be compared to predict Myy = V2e 2 / (8G f sin 2 9\y) . By 1983 when the W boson was first 
observed, the predicted mass including quantum corrections was given by M\y = 82 ±2.4 GeV to 
be compared with the UA1 Collaboration's measurement M\y = 81 ± 5 GeV [IB] . A few months 
later the Z° boson was observed Details of the W boson's experimental mass determination 
will be discussed in Section 4.4 as it is relevant for this thesis's contributions in Chapters [5] - [8] 



The Standard Model described in here has one serious difficulty regarding fine tuning of the 
Higgs mass radiative corrections. The Higgs field has two contributions to the mass self-energy 



shown in Fig. 2.2 The physical mass mu is given by 



mlK<H . E ^W1 A2+ _^ A2 (M0) 



f=u,d,e 

where m Qj u is the bare Higgs mass, A is the cut off energy scale, Y u Y d Y e are the three Yukawa 
coupling matrices, third term is from the Higgs loop. Assuming the physical Higgs mass mu is 
near the electo-weak scale (~ 100 GeV)[^]and A is at the Plank scale means that these terms need 
to cancel to some 34 orders of magnitude. If we make the cut off at A ~ 1 TeV the corrections 



from either Higgs-fermion loop 14 are equal to a physical Higgs mass of run = 100 GeV. Once 
again there is no reason that we cannot exclude a priori the possibility of a cancelation between 
terms of this magnitude. However, arguing that such a cancelation is unnatural has successfully 



predicted the charm quark mass in Sec 2.1.4 Arguing such fine-tuning is unnatural in the Higgs 



sector suggests new physics should be present below around 1 TeV. 

2.2 Dark Matter: Evidence for New, Massive, Invisible Particle 
States 

Astronomical observations also point to a something missing in the Standard Model. As- 
tronomers see evidence for 'dark matter' in galaxy rotation curves, gravitational lensing, the 



13 We use 100 GeV as a general electroweak scale. The current fits to Higgs radiative corrections suggest 
mu = 95 ± 35 GeV |8, ChlO], but the direct search limits require mn ^115 GeV. 

14 We use the fermion loop because the Higgs loop depends linearly on the currently unknown quartic Higgs self 
coupling A. 
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cosmic microwave background (CMB), colliding galaxy clusters, large-scale structure, and high 
red-shift supernova. Some recent detailed reviews of astrophysical evidence for dark-matter can 
be found in Bertone et.al[50 and Baer and Tata [51]. We discuss the evidence for dark matter 
below. 

• Galactic Rotation Curves 

Observations of the velocity of stars normal to the galactic plane of the Milky Way led Jan Oort 
|52| in 1932 to observe the need for 'dark matter'. The 1970s saw the beginning of using doppler 
shifts of galactic centers versus the periphery. Using a few ideal assumptions, we can estimate 
how the velocity should trend with changing r. The 21 cm HI line allows rotation curves of 
galaxies to be measured far outside the visible disk of a galaxy [53] To get a general estimate 
of how to interpret the observations, we assume spherical symmetry and circular orbits. If the 
mass responsible for the motion of a star at a radius r from the galactic center is dominated by 
a central mass M contained in R <C r, then the tangential velocity is vt ~ \J Gn M/r where 
Gat is Newton's constant. If instead the star's motion is governed by mass distributed with a 
uniform density p then vt = r \J 4 Gn /?7r/3. Using these two extremes, we can find when the 
mass of the galaxy has been mostly bounded by seeking the distance r where the velocity begins 
to fall like y/r. The observations |53| l54"j show that for no r does vt fall as y/l/r. Typically 
vt begins to rise linearly with r and then for some r > R Q it stabilizes at vt = const.. In our 
galaxy the constant vt is approximately 220km/sec. The flat vt vs r outside the optical disk of 
the galaxy implies p(r) oc r~ 2 for the dark matter. The density profile of the dark matter near 
the center is still in dispute. The rotation curves of the galaxies is one of the most compelling 
pieces of evidence that a non-absorbing, non-luminous, source of gravitation permeates galaxies 
and extends far beyond the visible disk of galaxies. 

• Galaxy clusters 

In 1937 Zwicky [55] published research showing that galaxies within clusters of galaxies were 
moving with velocities such that they could not have been bound together without large amounts 
of 'dark matter'. Because dark matter is found near galaxies, a modified theory of newtonian 
gravity (MOND) [56] has been shown to agree with the galactic rotation curves of galaxies. 
A recent observation of two clusters of galaxies that passed through each other, known as 
the bullet cluster, shows that the dark matter is separated from the luminous matter of the 
galaxies[57 a . The dark matter is measured through the gravitational lensing effect on the shape 
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of the background galaxies. The visible 'pressureless' stars and galaxies pass right through each 
other during the collision of the two galaxy clusters. The dark matter lags behind the visible 
matter. This separation indicates that the excess gravitation observed in galaxies and galaxy 
clusters cannot be accounted for by a modification of the gravitation of the visible sector, but 
requires a separate 'dark' gravitational source not equal to the luminous matter that can lag 
behind and induce the observed gravitational lensing. The bullet cluster observation is a more 
direct, empirical observation of dark matter. 

The dark matter seen in the rotation curves and galaxy clusters could not be due to H 
gas (the most likely baryonic candidate) because it would have been observed in the 21 cm 
HI observations and is bounded to compose no more than 1% of the mass of the system [58]. 
Baryonic dark matter through MAssive Compact Halo Objects (MACHOs) is bounded to be 
< 40% of the total dark matter |SJ Ch22] 

• Not Baryonic: The Anisotropy of the Cosmic Microwave Background and 
Nucleosythesis 

The anisotropy of the cosmic microwave background radiation provides an strong constraint 
on both the total energy density of the universe, but also on the dark-matter and the baryonic 
part of the dark matter. Here is a short explanation of why. 

The anisotropy provides a snapshot of the density fluctuations in the primordial plasma at 
the time the universe cools enough so that the protons finally 'recombine' with the electrons to 
form neutral hydrogen. The initial density fluctuations imparted by inflation follow a simple 
power law. The universe's equation-of-state shifts from radiation dominated to matter domi- 
nated before the time of recombination. Before this transition the matter-photon-baryon system 
oscillates between gravitational collapse and radiation pressure. After matter domination, the 
dark matter stops oscillating and only collapses, but the baryons and photons remain linked and 
oscillating. The peaks are caused by the number of times the baryon-photon fluid collapses and 
rebounds before recombination. 

The first peak of the CMB anisotropy power spectrum requires the total energy density of 
the universe to be very close to the critical density $7 = 1. This first peak is the scale where the 
photon-baryon plasma is collapsing following a dark-matter gravitational well, reaches its peak 
density, but doesn't have time to rebound under pressure before recombination. This attaches 
a length scale to the angular scale observed in the CMB anisotropy in the sky. The red-shift 
attaches a distance scale to the long side of an isosceles triangle. We can determine the geometry 
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from the angles subtended at vertex of this isosceles triangle with all three sides of fixed geodesic 
length. The angle is larger for such a triangle on a sphere than on a flat surface. By comparing 
the red-shift, angular scale, and length scale allows one to measure the total spatial curvature 
of space-time [59] to be nearly the critical density $7 = 1. 

We exclude baryonic matter as being the dark matter both from direct searches of H and 
MACHOs and from cosmological measurement of Qm vs £Im is the fraction of the critical 
density that needs to obey a cold-matter equation of state and is the fraction of the critical 
density that is composed of baryons. Q\ is the fraction of the critical density with a vacuum-like 
equation of state (dark energy). The anisotropy of the Cosmic Microwave Background (CMB) 
provides a constraint the reduced baryon density {h 2 ^) and reduced matter density (h 2 Q),). 

As the dark-matter collapses, the baryon-photon oscillates within these wells. The baryons 
however have more inertia than the photons and therefore 'baryon-load' the oscillation. The 
relative heights of the destructive and constructive interference positions indicate how much 
gravitational collapse due to dark matter occur during the time of a baryon-photon oscillation 
cycle. The relative heights of different peaks measure the baryon loading. Therefore, the third 
peak provides data about the dark matter density as opposed to the photon-baryon density. 

Using these concepts, the WMAP results give estimates of the baryonic dark-matter: Q b h 2 = 
0.0224±0.0009 and the dark-matter plus baryonic matter Sl M h 2 = 0.135±0.008 [50]. The CMB 
value for Of, agrees with the big-bang nucleousythesis (BBN) value 0.018 < ^h 2 < 0.023. The 
dark-matter abundance is found to be Qm = 0.25 using h = 0.73±0.03 |8, Ch21]. Combining the 
CMB results with supernova observations with baryon acoustic oscillations (BAO) compared to 
the galaxy distribution on large scales all lead to a three-way intersection on Qm vs S1a that give 
O m = 0.245 ±0.028, ft A = 0.757 ±0.021, n b = 0.042 ±0.002 [60]. This is called the concordance 
model or the ACDM model. 

Because all these tests confirm << Qm, we see that the cosmological observations confirm 
the non-baryonic dark matter observed in galactic rotation curves. 

• Dark Matter Direct Detection and Experimental Anomalies 

First we assume a local halo density of the dark matter of 0.3 GeV/ cm 3 as extracted from 
Qm and galaxy simulations. For a given mass of dark matter particle, one can now find a number 
density. The motion of the earth through this background number density creates a flux. Cryo- 
genic experiments shielded deep underground with high purity crystals search for interactions 
of this dark matter flux with nucleons. Using the lack of a signal they place bounds on the cross 
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Figure 2.3: Direct dark matter searches showing the limits from the experiments CDMS, 
XENON10, and the DAMA signal. Also shown are sample dark-matter cross sections and 
masses predicted for several supersymmetric models and universal extra dimension models. Fig- 
ures generated using 

http://dmtools.berkeley.edu/limitplots/ . 



section as a function of the mass of the dark matter particle which can be seen in Fig. 2.3 This 
figure shows direct dark-matter limits from the experiments CDMS (2008) [BT], XEONON10 
(2007) [62], and the DAMA anomaly (1999) [63 [64]. Also shown are sample dark-matter cross 
sections and masses predicted for several supersymmetric models |65j [66] and universal extra 
dimension models [67]. This shows that the direct searches have not excluded supersymmetry 
as a viable source for the dark matter. 

Although direct searches have no confirmed positive results, there are two anomalies of 
which we should be aware. The first is an excess in gamma rays observed in the EGRET 
experiment [68 j which points to dark matter particle annihilation with the dark-matter particles 
mass between 50 GeV < Mjv < 100 GeV. The uncertainties in the background gamma ray 
sources could also explain this excess. The second anomaly is an annual variation observed in 
the DAMA experiment |JH] . The annual variation could reflect the annual variation of the dark 
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matter flux as the earth orbits the sun. The faster and slower flux of dark-matter particles 
triggering the process would manifest as an annual variation. Unfortunately, CDMS did not 
confirm DAMA's initial results [69] {61]. This year the DAMA/LIBRA experiment released new 
results which claim to confirm their earlier result [63] with an 8d confidence. Until another 
group's experiment confirms the DAMA result, the claims will be approached cautiously. 

• Dark matter Candidates 

There are many models that provide possible dark-matter candidates. To name a few possi- 
bilities that have been considered we list: imparity conserving supersymmetry, Universal Extra 
Dimensions (UED), axions, degenerate massive neutrinos [70], stable black holes. In Chapter 
[3] we focus on Supersymmetry as the model framework for explaining the dark-matter, but our 
results in Chapters [6] - [8] apply to any model where new particle states are pair produced and 
decay to two semi-stable dark matter particles that escape the detector unnoticed. 

The pair-produced nature of dark-matter particles is relatively common trait of models with 
dark matter. For example the UED models |71[ I72j 173] can also pair produced dark- matter 
particles at a collider. The lightest Kaluza-Klein particle (LKP) is also a dark-matter candidate. 
In fact UED and SUSY have very similar hadron collider signatures [73J . 

2.3 Supersymmetry: Predicting New Particle States 

Supersymmetry will be the theoretical frame-work for new-particle states on which this thesis 
focuses. Supersymmetry has proven a very popular and powerful framework with many suc- 
cesses and unexpected benefits: Supersymmetry provides a natural dark-matter candidate. As 
a fledgling speculative theory, supersymmetry showed a preference for a heavy top-quark mass. 
In a close analogy with the GIM mechanism, Supersymmetry's minimal phenomenological im- 
plementation eliminates a fine-tuning problem associated with the Higgs boson in the Standard 
Model. Supersymmetry illuminates a coupling constant symmetry (Grand Unification) among 
the three non-graviational forces at an energy scale around 2 x 10 16 GeV. Supersymmetry is 



the only extension of Poincare symmetry discussed in 2.1.1 allowed with a graded algebra [75] . 
It successfully eliminates the tachyon in String Theory through the GSO projection [76]. Last, 
it is a candidate explanation for the 3 a deviation of the muons magnetic moment known as the 
anomaly [77] . 

These successes are exciting because they follow many of the precedents and clues described 
earlier in this chapter that have successfully predicted new-particle states in the past. SUSY, 
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short for supersymmetry, relates the couplings of the Standard-Model fermions and bosons to the 
couplings of new bosons and fermions known as superpartners. SUSY is based on an extension of 



the Poincare-symmetry (Sec |2.1.T ). In the limit of exact Supersymmetry, the Higgs self-energy 



problem (Sec |2.1.5 ) vanishes. In analogy to the GIM mechanism (Sec |2.1.4), the masses of the 



new SUSY particles reflect the breaking of supersymmetry. There are many theories providing 
the origin of supersymmetry breaking which are beyond the scope of this thesis. The belief that 
nature is described by Supersymmetry follows from how SUSY connects to the past successes 
in predicting new-particle states and their masses from symmetries, broken-symmetries, and 
fine-tuning arguments. 

Excellent introductions and detailed textbooks on Supersymmetry exist, and there is no 
point to reproduce these textbooks here. Srednicki provides a very comprehensible yet compact 
introduction to supersymmetry via superspace at the end of Ref [39]. Reviews of SUSY that 
have proved useful in developing this thesis are Refs [TH] EH [HQl Ell E2 [83] . In this introduction, 
we wish to highlight a few simple parallels to past successes in new-particle state predictions. 
We review supersymmetric radiative corrections to the coupling constants and to the Yukawa 
couplings. These radiative corrections can be summarized in the renormalization group equations 
(RGE). The RGE have surprising predictions for Grand Unification of the non-gravitational 
forces and also show evidence of a quark lepton mass relations suggested by Georgi and Jarlskog 
which we refer to as mass matrix unification. The mass matrix unification predictions can be 
affected by potentially large corrections enhanced by the ratio of the vacuum expectation value 
of the two supersymmetric Higgs particles v u /V(i = tan/3. Together, these renormalization group 
equations and the large tan (3 enhanced corrections provide the basis for the new contributions 
this thesis presents in Chapter [3} 

2.3.1 Supersymmetry, Radiative Effects, and the Top Quark Mass 



Supersymmetry extends the Poincare symmetry used in 2.1.1 The super- Poincare algebra 
involves a graded Lie algebra which has generators that anticommute as well as generators that 
commute. One way of understanding Supersymmetry is to extend four space-time coordinates to 
include two complex Grassmann coordinates which transform as two-component spinors. This 
extended space is called superspace. In the same way that generates space-time translations, 



we introduce one supercharga| 15 | Q that generates translations in 9 and 6. The graded super- 



3 We will only consider theories with one supercharge. 
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Poincare algebra includes the Poincare supplemented by 



Ql,P*\ =[Q a ,P"] 

{Qa,Q b } = {QtQ\} 

[Q a , MH 

~Ql,M^' 



2 rr M P 



(2.21) 
(2.22) 
(2.23) 
(2.24) 
(2.25) 



where M are the boost generators satisfying Eq(2.3), {, } indicate anticommutation, a, b, a, b are 



spinor indices, = (1, a) and <r M = (1, — a). Eq(2.21 ) indicates htat the charge Q is conserved 



by space-time translations. Eq(2.22 ) indicates that no more than states of two different spins can 



be connected by the action of a supercharge. Eq(2.23 2.24) indicate that Q and transform as 



left and right handed spinors respectively. Eq(2.25) indicates that two supercharge generators 
can generate a space-time translation. 

In relativistic QFT, we begin with fields <j)(x) that are functions of space-time coordinates x^ 
and require the Lagrangian to be invariant under a representation of the Poincare group acting 
on the space-time and the fields. To form supersymmetric theories, we begin with superfields 
<J>(x, 9, 9) that are functions of space-time coordinates x^ and the anticommuting Grassmann 
coordinates 9 a and 9 a and require the Lagrangian to be invariant under a representation of the 
super-Poincare symmetry acting on the superspace and the superfields. The procedure described 
here is very tedious: defining a representation of the super-Poincare algebra, and formulating 
a Lagrangian that is invariant under actions of the group involves many iterations of trial and 
error. Several short-cuts have been discovered to form supersymmetric theories very quickly. 
These shortcuts involve studying properties of superfields. 

Supersymmetric theories can be expressed as ordinary relativistic QFT by expressing the 
superfield $(x, 9, 9) in terms of space-time fields like <j>(x) and ip(x). The superfields, which we 
denote with hats, can be expanded as a Taylor series in 9 and 9. Because 9\ = 0? = 0, the 
superfield expansions consist of a finite number of space-time fields (independent of 9 or 9) ; some 
of which transform as scalars, and some transform as spinors, vectors, or higher level tensors. 
A supermultiplet is the set of fields of different spin interconnected because they are part of the 
same superfield. If supersymmetry were not broken, then these fields of different spin would be 
indistinguishable. The fields of a supermultiplet share the same quantum numbers (including 
mass) except spin. 
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Different members of a supermultiplet are connected by the action of the supercharge op- 
erator Q or roughly speaking Q\ boson) = |fermion) and Q| fermion) = | boson). Because 
Qa = (Qh) 2 = a supermultiplet only has fields of two different spins. Also because it is a 
symmetry transformation, the fields of different spin within a supermultiplet need to have equal 
numbers of degrees of freedom [^} A simple type of superfield is a chiral superfield. A chiral 
superfield consists of a complex scalar field and a chiral fermion field each with 2 degrees of 
freedom. Another type of superfield is the vector superfield V = V* which consists of a vector 
field and a Weyl fermion field. Magically vector superfields have natural gauge transformations, 
the spin-1 fields transform as connections under gauge transformation whereas the superpartner 
spin-1/2 field transforms covariantly in the adjoint of the gauge transformation^] 

The shortcuts to form Lagrangians invariant under supersymmetry transformations are based 
on three observations: (1) the products of several superfields is again a superfield (2) the term 
in an expansion of a superfield (or product of superfields) proportional to Q a Q\,e ah is invariant 
under supersymmetry transformations up to a total derivative (called an F-term), and (3) the 
term in an expansion of a superfield (or product of superfields) proportional to 9 a 6i,e ab 9a0j ) e ab is 
invariant under supersymmetry transformations up to a total derivative (called a D-tevm). 

These observations have made constructing theories invariant under supersymmetry a rela- 
tively painless procedure: the D-term of ^e~ v ^ provides supersymmetricly invariant kinetic 
terms. A superpotential W governs the Yukawa interactions among chiral superfields. The F- 
term of W gives the supersymmetricly invariant interaction Lagrangian. The i^-term of W can 
be found with the following shortcuts: If we take all the chiral superfields in W to be enumer- 
ated by i in Vt'j then the interactions follow from two simple calculations: The scalar potential 
is given by C D= —Q2j \dW/d^j\ 2 ) where is a the superfield and all the superfields inside 
() are replaced with their scalar part of their chiral supermultiplet. The fermion interactions 
with the scalars are given by C D —^2 i j(d 2 W/d^id^j)TpiTpj. where the superfields in () are 
replaced with the scalar part of the chiral supermultiplet and and ipj are the 2-component 
Weyl fermions that are part of their chiral supermultiplet i$fj. Superpotential terms must be 
gauge invariant just as one would expect for terms in the Lagrangian and must be holomor- 
phic function of the superfielda/^j Another way to express the interaction Lagrangian is by 



16 There are auxiliary fields in supermultiplets that, while not dynamical, preserve the degrees of freedom when 
virtual states go off mass shell. 

17 Although this seems very unsymmetric, Wess and BaggerJST] show a supersymmetric differential geometry 
with tetrads that illuminate the magic of how the spin-f fields transform as connections but the superpartners 
transform covariantly in the adjoint representation. 

18 By holomorphic we mean the superpotential can only be formed from unconjugated superfields 4> and not 
the conjugate of superfields like 9*. 
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Table 2.2: Minimial Supersymmetric Standard Model (MSSM) field's transformation properties. 
Fields grouped together are part of the same supermultiplet. * Indicated fields transform as a 
connection as opposed to covariantljff^. 

£w = / d6 2 W + f d9 2 W> where the integrals pick out the F terms of the superpotential W. 

In addition to supersymmetry preserving terms, we also need to add 'soft' terms which 
parameterize the breaking of supersymmetry. 'Soft' refers to only SUSY breaking terms which 
do not spoil the fine-tuning solution discussed below. 

To form the minimal supersymmetric version of the Standard Model, known as the MSSM, 
we need to identify the Standard Model fields with supermultiplets. The resulting list of fields 



is given in Table 2.2. Supersymmetry cannot be preserved in the Yukawa sector with just 
one Higgs field because the superpotential which will lead to the F-terms in the theory must be 
holomorphic so we cannot include both H and H* superfields in the same superpotential; instead 
the MSSM has two Higgs fields H u = and Hd = (h d ,h^). The neutral component of 

each Higgs will acquire a vacuum expectation value (VEV): = v u and {h° d ) = Vd- The 
parameter tan/3 = v u /vd is ratio of the VEV of the two Higgs fields. 

The structure of the remaining terms can be understood by studying a field like a right- 
handed up quark. The (u c ) transforms as a 3 under SU(3) so its superpartner must be also 
transform as a 3 but have spin or 1. No Standard Model candidate exists that fits either 
option. A spin-1 superpartner is excluded because the fermion component of a vector superfield 
transforms as a connection in the adjoint of the gauge group; not the fundamental representation 
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Figure 2.4: One Loop Contributions to the supersymmetric Higgs H u mass self energy, (a) top 
loop, (b) stop loop. 

like a quark. If (u c ) is part of a chiral superfield, then there is an undiscovered spin-0 partner. 
Thus the (u c ) is part of a chiral supermultiplet with a scalar partner called a right-handed squark 
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The remaining supersymmetric partner states being predicted in Table 2.2 can be deduced 
following similar arguments. The superpartners are named after their SM counterparts with 
preceding 's' indicating it is a scalar superpartner of a fermion or affixing 'ino' to the end indi- 
cating it is a fermionic partner of a boson: for example selectron, smuon, stop-quark, Higgsino, 
photino, gluino, etc. 

Supersymmetry solves the fine-tuning problem of the Higgs self energy. The superpotential 
describing the Yukawa sector of the MSSM is given by 

W = Y^VH^ey + YjQ'Haidy + Y^Q'H^uy + fj,H u H d (2.26) 

where the fields with hats like L, H Ut d, Q, etc. are all superfields. In the limit of exact supersym- 
metry the fine-tuning problem is eliminated because the resulting potential and interactions lead 



to a cancelation of the quadratic divergences between the fermion and scalar loops in Fig 2.4 
The self-energy of the neutral Higgs H u = (h^,h°) is now approximately given by 

9 . „ Tr(T"(T u )t) a9 Tr(y"(y u )t) „ , ^ 

where |// D | 2 is the modulus squared of the bare parameter in the superpotential which must 
be positive, the second term comes from the fermion loop and therefore has a minus sign, 
and the third term comes from the scalar loop. Both divergent loops follow from the MSSM 
superpotential: the first loop term follows from the fermion coupling d 2 W/dtdt c = yth° u tt c and 
the second loop term follows from the scalar potential |3VF/9t c | 2 = y||/i°| 2 |t| 2 and |<9Wy<9i| 2 = 
y||/i°| 2 |t c | 2 where we assume the top-quark dominates the process. Exact supersymmetry ensures 
these two quadratically divergent loops cancel. However two issues remain and share a common 
solution: supersymmetry is not exact, and the Higgs mass squared must go negative to trigger 



The right-handed refers to which fermion it is a partner with. The field is a Lorentz scalar. 
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spontaneous symmetry breaking (SSB). In the effective theory well above the scale where all 
superpartners are energetically accessible the cancelation dominates. The fine-tuning arguments 



in Sec. 2.1.5 suggest that if fi a ~ 100 GeV 20 this cancelation should dominate above about 1 
TeV. In an effective theory between the scale of the top-quark mass and the stop-quark mass 
only the fermion loop (Fig |2.4| a) will contribute significantly. In this energy-scale region we 
neglect the scalar loop (Fig |2.4| b) . With only the fermion loop contributing significantly and if 
yt is large enough then the fermion loop will overpower /i 2 then the mass squared m 2 D can be 
driven negative. 

In this way the need for SSB without fine tuning in the MSSM prefers a large top-quark mass 
and the existence of heavier stop scalar states. Assuming fj, ~ 100 GeV predicts t and t c below 
around 1 TeVp"] This phenomena for triggering SSB is known as radiative electroweak symmetry 
breaking (REWSB). We have shown a very coarse approach to understanding the major features; 
details can be found in a recent review jH5j or any of the supersymmetry textbooks listed above. 

The detailed REWSB |85] technology was developed in the early 1980s, and predicted Mt ~ 
55 — 200 GeV [86] [87 ; a prediction far out of the general expectation of M t ~ 15 — 40 GeV of 
the early 1980s and closer to the observed value of Mt = 170.9 ±1.9 GeV. If supersymmetric 
particles are observed at the LHC, the large top-quark mass may be looked at as the first 
confirmed prediction of supersymmetry. 

The MSSM provides another independent reason to prefer a large top quark mass. The 
top Yukawa coupling's radiative corrections are dominated by the difference between terms 
proportional to g 2 and y 2 . If the ratio of these two terms is fixed, then yt will remain fixed [92 . 
For the standard model this gives a top quark mass around 120 GeV. However for the MSSM, 
assuming 03 = g\j^ = 0.12 then one finds a top quark mass around 180 GeV assuming a 
moderate tan (5 [93] . Our observations of the top quark mass very near this fixed point again 
points to supersymmetry as a theory describing physics at scalers above the electroweak scales. 

2.3.2 Dark Matter and SUSY 

Supersymmetry is broken in nature. The masses of the superpartners reflect this breaking. Let's 
assume the superpartner masses are at a scale that avoids bounds set from current searches yet 
still solve the fine-tuning problem. Even in this case there are still problems that need to be 



20 Again we choose 100 GeV as a generic electroweak scale. 

21 As fine tuning is an aesthetic argument, there is a wide range of opinions on the tolerable amount of fine 
tuning acceptable. There is also a wide range of values for fi that are tolerable. 

22 Raby [88] Glashow [89] and others [90] |91] all made top-quark mass predictions in the range of 15 — 40 GeV. 
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resolved 23 There are many couplings allowed by the charge assignments displayed in Table 2.2 
that would immediately lead to unobserved phenomena. For example the superpotential could 
contain superfield interactions u c u c d c or QLd c or LLe c or kLH u where k is a mass scale. Each 
of these interactions is invariant under the SU(3) C x SU{2)l x U{1)y charges listed in Table 



2.2 These couplings, if allowed with order 1 coefficients, would violate the universality of the 
four-Fermi decay, lead to rapid proton decay, lepton and baryon number violation, etc. 

These couplings can be avoided in several ways. We can require a global baryon number 
or lepton number U(l) on the superpotential; in the Standard Model these were accidental 
symmetries. However, successful baryogenesis requires baryon number violation so imposing it 
directly is only an approximation. Another option is to impose an additional discrete symmetry 
on the Lagrangian; a common choice is i?-parity 

R = (_l)2j+3B+L ^^g) 

where j is the spin of the particle. This gives the Standard Model particles R = 1 and the 
superpartners R = — 1. Each interaction of the MSSM Lagrangian conserves R-parity. The 
specific choice of how to remove these interactions is more relevant for GUT model building. 
The different choices lead to different predictions for proton decay lifetime. 

The i?-Parity, which is needed to effectively avoid these unobserved interactions at tree level, 
has the unexpected benefit of also making stable the lightest supersymmetric particle (LSP). A 
stable massive particle that is non-baryonic is exactly what is needed to provide the dark matter 



observed in the galactic rotation curves of Sec. 2.2 



2.3.3 Renormalization Group and the Discovered Supersymmetry Symme- 
tries 

Unification of SU(3) C x SU(2)l x U(1)y coupling constants 

Up until now, the divergent loops have been treated with a UV cut off. Renormalization of non- 
abelian gauge theories is more easily done using dimensional regularization where the dimensions 
of space-time are taken to be d = 4 — 2e. The dimensionless coupling constants g pick up a 
dimensionfull coefficient g fi~ e where /U is an arbitrary energy scale. The divergent term in 
loops diagrams is now proportional to 1/e, and the counter-terms can be chosen to cancel these 
divergent parts of these results. By comparison with observable quantities, all the parameters 
in the theory are measured assuming a choice of /U. The couplings with one choice of /x can 
be related to an alternative choice of \i by means of a set of differential equation known as the 



There is also a flavor changing neutral current (FCNC) problem not discussed here. 
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Running of the Coupling Constants 
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Figure 2.5: Gauge couplings for the three non-gravitational forces as as a function of energy 
scale for the (left) Standard Model and (right) MSSM. 

renormalization group equations (RGE). The choice of // is similar to the choice of the zero of 
potential energy; in principle any choice will do, but in practice some choices are easier than 
others. Weinberg shows how this arbitrary scale can be related to the typical energy scale of a 
process [33]. 

The renormalization group at one loop for the coupling constant g of an SU (N) gauge theory 
coupled to fermions and scalars is 
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(2.29) 



where np is the number of 2-component fermions, ns is the number of complex scalars, S(Rf,s) 
is the Dynkin index for the representation of the fermions or scalars respectively. Applying this 
to our gauge groups in the MSSM C(G) = N for SU(N) and C(G) = for [7(1). For the 
fundamental SU(N) we have S(R) = 1/2; for the adjoint S(R) = N; for U(l) we have S(R) 
equal to the ^2,{Y/2) 2 over all the scalars (which accounts for the number of scalars). When 
the SU{3) C x (577(2) £ x U(l)y is embedded in a larger group like SU(5), the U(l)y coupling g' 
is rescaled to the normalization appropriate to the SU(5) generator that becomes hypercharge. 
This rescaling causes us to work with g\ = a/5 / 3g' pi 

Applying this general formula to both the Standard Model (SM) and the MSSM leads to 



d 




dlogfi 
n G 4/3 + 1/10 
-22/3 + n G 4/3 + n H 1/6 
-11 + no 4/3 



9i 



MSSM 



(2.30) 



h= n G 2 + n H 3/10 
b 2 = -6 + n G 2 + n H l/2 
b 3 = -9 + n G 2 



where n G = 3 is the number of generations and nn is the number of Higgs doublets (jih = 1 



in SM and nu = 2 in MSSM). A miracle is shown in Fig. 2.5 The year 1981 saw a flurry of 



Eq(2.29l is for g' and one must substitute the definition of gi to arrive at Eq(2.30l 
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papers from Dimopoulos, Ibanez, Georgi, Raby, Ross, Sakai and Wilczek who were detailing 
the consequences of this miracle [94J[95 [96 [97J. There is one degree of freedom in terms of 
where to place an effective supersymmetry scale Ms where Standard Model RG running turns 
into MSSM RG running. At one-loop order, unification requires Ms ~ Mz', at two-loop order 
200 GeV < M s < 1 TeV 



2 r , 



The MSSM was not designed for this purpose, but the particle spectrum gives this result 



effortlessl} 26 A symmetry among the coupling of the three forces is discovered through the RG 
equations. If the coupling unify, they may all originate from a common grand-unified force that 
is spontaneously broken at \x ~ 2 x 10 16 GeV. 

There are several possibilities for SUSY GUT theories: SU(5) [TOO] or 50(10) or SU(4) x 
SU{2)l x SU(2)r |101] to list just a few. SU(5), although the minimal version is now excluded 
experimentally, is the prototypical example with which we work. A classic review of GUT 
theories can be found in [102J 

Georgi-Jarlskog Factors 

It is truly miraculous that the three coupling constants unify (to within current experimental 
errors) with two-loop running when adjusted to an SU(5) grand unified gauge group and when 
the SUSY scale is placed in a region where the fine-tuning arguments suggest new-particle states 
should exist. Let's now follow the unification of forces arguments to the next level. Above the 
unification scale, there is no longer a distinction between SU(3) C and t/(l)y. If color and 
hypercharge are indistinguishable, what distinguishes an electron from a down-quark? The 
Yukawa couplings, which give rise to the quark and lepton masses, are also functions of the scale 
H and RG equations relate the low-energy values to their values at the GUT scale. Appendix |A| 
gives details of the RG procedure used in this thesis to take the measured low-energy parameters 
and use the RG equations to relate them to the predictions at the GUT scale. Do the mass 
parameters also unify? 

With much more crude estimates for the quark masses, strong force, and without the knowl- 
edge of the value of the top quark mass Georgi and Jarlskog (GJ) [TT] noticed that at the GUT 



25 This range comes from a recent study [5S] which assumes as{Mz) = 0.122. If we assumes as(Mz) = 0.119, 
then we find 2 TeV < Ms < 6 TeV. Current PDG [H Ch 10] SM global fits give a s = 0.1216 ± 0.0017. 

26 Very close coupling constant unification can also occur in non-supersymmetric models. The Standard Model 
with six Higgs doublets is one such example [55], but the unification occurs at too low a scale ~ 10 14 GeV. 
GUT-scale gauge-bosons lead to proton decay. Such a low scale proton decay at a rate in contradiction with 
current experimental bounds. 
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scale the masses satisfied the approximate relations 
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rrtft ~ 3m s 3m e 



m d . 



(2.31) 



This is a very non-trivial result. The masses of the quarks and charged leptons span more than 
5 orders of magnitude. The factor 3 is coincident ally equal to the number of colors. At the scale 
of the Z°-boson's mass \i = Mz the ratios look like ~ 2m s , and 1.6m T ~ nib so the factor 
of three is quite miraculous. Using this surprising observation, GJ constructed a model where 
this relation followed from an SU(5) theory with the second generation coupled to a Higgs in a 
different representation. 

In the SU(5) model the fermions are arranged into a 5 (tpf) and a 10 (ip a bi) where a, b, c, .. 
are the SU(5) indexes and i,j, .. are the family indexes. The particle assignments are 
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(2.32) 
(2.33) 



where c indicates the conjugate field. There are also a 5 Higgs fields (Hd) a and a 5 Higgs 



field (H u ) a . The key to getting the mass relations hypothesized in Eq(2.31) is coupling only 



the second generation to a 45 Higgs 28 The VEVs of the Higgs fields are H u = (0, 0, 0, 0, v u ), 
Hd = (0, 0, 0, 0, Vd) and (^45,^)^5 = v^(5^ — 45 a4 5b4)- The coupling to matter that gives mass 
to the down-like states is 



w Yd = Y d5ij r bl iPi(H d )b + Y d45iJ r bi 4(H dA5 y c 

and that give mass to the up-like states 

W Yu = Y 5uij r bl ^ cdj (H u ) e € abcde . 



(2.34) 



(2.35) 



Georgi and Jarlskog do not concern themselves with relating the neutrino masses to the up- 
quark masses, so we will focus on the predictions for the down-like masses. The six masses of 
both the down-like quarks and the charged leptons may now be satisfied by arranging for 
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(2.36) 



27 To the best of my knowledge, the b = r relations were first noticed by Buras et.al 103 . 

28 In tensor notation the 45 representation is given by (Hd : 4,s)ab where ab are antisymmetric and the five traces 



(Hd, 45)06 = are removed. 
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and fitting the three parameters A, B, and C. The fitting will create the hierarchy B » 



C » A. Coupling these Yukawa matricies to the Eq(2.34) gives the factor of — 3C for the (2,2) 
entry of the leptons mass matrix relative to the (2, 2) entry of the down-like quark mass matrix. 
Because B » C >> A, the equality of the (3,3) entry leads to m^/rar ~ 1. The (2,2) entry 
dominates the mass of the second generation so m^/m s ~ 3. The determinant of the resulting 
mass matrix is independent of C so at the GUT scale the product of the charged lepton masses 
is predicted to equal the product of the down-like quark masses. 

These results have been generalized to other GUT models like the Pati-Salam model FH 
Family symmetries have been used to arrange the general structure shown here |104j . The 
continued validity of the Georgi-Jarlskog mass relations is one of the novel contributions of the 
thesis presented in Chapter [3j 

2.3.4 tan/5 Enhanced Threshold Effects 

The Appelquist Carazzone |105| decoupling theorem indicates that particle states heavier than 
the energy scales being considered can be integrated out and decoupled from the low-energy 
effective theory. A good review of working with effective theories and decoupling relations is 
found in Pich |106j . The parameters we measure are in some cases in an effective theory of 
SU(S) C x U(1)em\ in other cases we measure the parameters with global fits to the Standard 
Model. At the energy scale of the sparticles, we need to match onto the MSSM effective theory. 
Finally at the GUT scale we need to match onto the GUT effective theory. 

As a general rule, matching conditions are needed to maintain the order of accuracy of the 
results. If we are using one-loop RG running, we can use trivial matching conditions at the 
interface of the two effective theories. If we are using two-loop RG running, we should use one- 
loop matching conditions at the boundaries. This is to maintain the expected order of accuracy 
and precision of the results. There is an important exception to this general rule relevant to 
SUSY theories with large tan (3. 

At tree level the VEV of H u gives mass to the up-like states (t,c,u) and the neutrinos. At 
tree-level the VEV of gives mass to the down-like states (b, s, d, r, [i, e). However 'soft' 
interactions which break supersymmetry allow the VEV of H u to couple to down-like Yukawa 
couplings through loop diagrams. Two such soft terms are the trilinear couplings C D ytA t H u tt c 
where At is a mass parameter and the gluino mass C D ggM^ where M3 is the gluino's soft mass 
parameter. 



29 In Pati-Salam the 45's VEV is such that one has a factor of 3 and not —3 for the charged leptons vs the 
down-like quark Yukawa coupling. 
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The matching conditions for two effective theories are deduced by expressing a common 
observable in terms of the two effective theories. For example the Pole mass M b of the bottom 



quark ^ u would be expressed as 




bL > J— bR % jg y b ~% >r 

m 6 « ! + — is — I — £ — + ■■■ 

" H 

where H takes on the vacuum expectation value, yb(n) and ^(/u) are the QFT parameters which 
depends on an unphysical choice of scale /i. The same observable expressed in MSSM involves 
new diagrams 



M b ^ ~[ + 

Hd 
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The parameters M (^), 9§ M {lA are not equal to the parameters y^ s (fj,) , g^f (//) . By 

expressing y^ ISSM = + <5y fe and likewise for g% lsSM = #f M + <5<?3, we find many common 

graphs which cancel. We are left with an expression for 5y b equal to the graphs not common 

between the two effective theories. 

Hu Hu 
i i 

- ! w ~ ! y t A * 

-(Sy b )v d = b ,^- s b + f,^ v ? +... (2.37) 

b L ,' M 3 ) b R b L [v- 1 b R 

*" J ~ ■ ' M Hft * ~ * ~ ^ 

03 9 9 93 y t h h y b 

The two graphs in this correction are proportional to the VEV of (h%) = v u . However the y b 
Yukawa coupling is the ratio of m b to (h° d ) = Vd- This makes the correction 5y b /yb due to the 
two loops shown proportional to v u /vd = tan/3. If tan/3 were small < 5, then the loop result 
times the tan ft would remain small and the effect would be only relevant at two-loop accuracy. 
However when tan ft > 10 then the factor of tan ft makes the contribution an order of magnitude 
bigger and the effect can be of the same size as the one-loop running itself. 



If it existed as a free state. 
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These tan (3 enhanced SUSY threshold corrections can have a large effect on the GUT-scale 
parameters. More precise observations of the low-energy parameters have driven the Georgi- 
Jarlskog mass relations out of quantitative agreement. However there is a class of tan /3 enhanced 
corrections that can bring the relations back into quantitative agreement. Chapter [3] of this the- 
sis makes predictions for properties of the SUSY mass spectrum by updating the GUT-scale 
parameters to the new low-energy observations and considering properties of tan [3 enhanced 
SUSY threshold corrections needed to maintain the quantitative agreement of the Georgi Jarl- 
skog mass relations. 

Chapter Summary 

In this chapter we have introduced the ingredients of the Standard Model and its supersymmetric 
extension and given examples of how symmetries, broken symmetries, and fine-tuning arguments 
have successfully predicted the mass of the positron, the Q~ , the charm quark, and the and 
Z° bosons. We have introduced astrophysical evidence that indicate a significant fraction of 
the mass-energy density of the universe is in a particle type yet to be discovered. We have 
introduced supersymmetry as a plausible framework for solving the fine-tuning of the Higgs self 
energy, for explaining the top-quarks large mass, and for providing a dark-matter particle. In 
addition we have discussed how SUSY predicts gauge coupling unification, and a framework 
for mass-matrix unification. Last we have introduced potentially large corrections to the RG 
running of the mass matrices. 
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Chapter 3 



Predictions from Unification and 
Fermion Mass Structure 

Chapter Overview 

Grand Unified Theories predict relationships between the GUT-scale quark and lepton masses. 
Using new data in the context of the MSSM, we update the values and uncertainties of the 
masses and mixing angles for the three generations at the GUT scale. We also update fits to 
hierarchical patterns in the GUT-scale Yukawa matrices. The new data shows that not all the 
classic GUT-scale mass relationships remain in quantitative agreement at small to moderate 
tan (3. However, at large tan /3, these discrepancies can be eliminated by finite, tan /3-enhanced, 
radiative, threshold corrections if the gluino mass has the opposite sign to the wino mass. This 
chapter is based on work first published by the author and his supervisor in Ref. |10j . 

Explaining the origin of fermion masses and mixings remains one of the most important goals 
in our attempts to go beyond the Standard Model. In this, one very promising possibility is 
that there is an underlying stage of unification relating the couplings responsible for the fermion 
masses. However we are hindered by the fact that the measured masses and mixings do not 
directly give the structure of the underlying Lagrangian both because the data is insufficient 
unambiguously to reconstruct the full fermion mass matrices and because radiative corrections 
can obscure the underlying structure. In this chapter we will address both these points in the 
context of the MSSM. 

We first present an analysis of the measured mass and mixing angles continued to the GUT 
scale. The analysis updates Ref [T] using the precise measurements of fermion masses and mixing 
angles from the b- factories and the updated top-quark mass from CDF and DO. The resulting 
data at the GUT scale allows us to look for underlying patterns which may suggest a unified 
origin. We also explore the sensitivity of these patterns to tan /3-enhanced, radiative threshold 
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Figure 3.1: Updates to the top-quark mass, strong coupling constant, and bottom-quark mass 
are responsible for the quantitative stress of the classic GUT relation for yb/y T - 
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Figure 3.2: Updates to the strange-quark mass and V c b are responsible for the quantitative 
stress of the Georgi- Jarlskog mass relations and the need to update values from Ref . [T] . 
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corrections. 

We next proceed to extract the underlying Yukawa coupling matrices for the quarks and 
leptons. There are two difficulties in this. The first is that the data cannot, without some 
assumptions, determine all elements of these matrices. The second is that the Yukawa coupling 
matrices are basis dependent. We choose to work in a basis in which the mass matrices are 
hierarchical in structure with the off-diagonal elements small relative to the appropriate combi- 



nations of on-diagonal matrix elements. Appendix [Bj Eq(B.l) defines this basis more precisely. 
This is the basis we think is most likely to display the structure of the underlying theory, for 
example that of a spontaneously broken family symmetry in which the hierarchical structure is 
ordered by the (small) order parameter breaking the symmetry. With this structure to leading 
order the observed masses and mixing angles determine the mass matrix elements on and above 
the diagonal, and our analysis determines these entries, again allowing for significant tan/3 en- 
hanced radiative corrections. The resulting form of the mass matrices provides the "data" for 
developing models of fermion masses such as those based on a broken family symmetry. 

3.1 Supersymmetric Thresholds and GUT-Scale Mass Relations 



The data set used is summarized in Table 3.1 Since the fit of reference [lj (RRRV) to the 



Yukawa texture was done, the measurement of the Standard-Model parameters has improved 



considerably. Figs 3.1 and 3.2 highlight a few of the changes in the data since: The top-quark 
mass has gone from M t = 174. 3±5 GeV to M t = 170.9±1.9 GeV. In 2000 the Particle Data Book 
reported m&(m;,) = 4.2 ± 0.2 GeV |109j which has improved to m^m^) = 4.2 ± 0.07 GeV today. 
In addition each higher order QCD correction pushes down the value of m\,{Mz) at the scale of 
the Z bosons mass. In 1998 m b (M z ) = 3.0 ±0.2 GeV [1113 and today it is m b {M z ) = 2.87 ±0.06 
GeV |lllj . The most significant shift in the data relevant to the RRRV fit is a downward revision 
to the strange-quark mass at the scale fj,L = 2 GeV from m s {ni,) ~ 120±50 MeV |109] to today's 
value m s {iJLL) = 103±20 MeV. We also know the CKM unitarity triangle parameters better today 
than six years ago. For example, in 2000 the Particle Data book reported sin 2(3 = 0.79 ± 0.4 



|109] which is improved to sin 2(3 = 0.69 ± 0.032 in 2006 [41 . Figures 3.1 and 3.2 show these 
updates visually. The sin 2(3 value is about 1.2 a off from a global fit to all the CKM data [112 , 
our fits generally lock onto the global-fit data and exhibit a 1 ex tension for sin 2(3. Together, the 
improved CKM matrix observations add stronger constraints to the textures compared to data 
from several years ago. 
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Low-Energy Parameter 



Value(Uncertainty in last digit(s)) 



Notes and Reference 



m u (fi L ) /m d (^L) 
m s (fi L )/md()iL) 
m u (fi L ) + m d (fi L ) 

m s~( m d + m u) 2 /^ 

m 3 (/LJt) 

m u {pL L ) 
m d (p L ) 
m c (m c ) 
m b (m b ) 
M t 

(M e , M/j,, M T ) 

A Wolfenstein parameter 
p Wolfenstein parameter 
A Wolfenstein parameter 
r\ Wolfenstein parameter 

| Vckm\ 

sin 2/3 from CKM 
Jarlskog Invariant 



0.45(15) 
19.5(1.5) 
[8.8(3.0), 7.6(1.6)] MeV 

22.8(4) 
[103(20), 95(20)] MeV 

3(1) MeV 

6.0(1.5) MeV 

1.24(09) GeV 

4.20(07) GeV 

170.9 (1.9) GeV 
(0.511(15), 105.6(3.1), 1777(53) ) MeV 

0.818(17) 
0.221(64) 
0.2272(10) 
0.340(45) 

^0.97383(24) 0.2272(10) 0.00396(09) N 

0.2271(10) 0.97296(24) 0.04221(80) 
v 0.00814(64) 0.04161(78) 0.999100(34), 
0.687(32) 
3.08(18) x 10" 5 



PDB Estimation [41] 

PDB Estimation [41] 

PDB, Quark Masses, pg 15 [IT]. 

( Non-lattice, Lattice ) 

Martemyanov and Sopov |107] 

PDB, Quark Masses, pg 15 [41]. 

[Non-lattice, lattice] 

PDB, Quark Masses, pg 15 gj. 

Non-lattice. 

PDB, Quark Masses, pg 15 [41]. 
Non-lattice. 

PDB, Quark Masses, pg 16 [IT]. 
Non-lattice. 

PDB, Quark Masses, pg 16,19 

[IT] . Non-lattice. 

CDF & DO [I08j Pole Mass 

3% uncertainty from neglecting 

Y e thresholds. 

PDB Ch 11 Eq. 11.25 [UJ 

PDB Ch 11 Eq. 11.25 [IT] 

PDB Ch 11 Eq. 11.25 [UJ 

PDB Ch 11 Eq. 11.25 [UJ 

PDB Ch 11 Eq. 11.26 [UJ 

PDB Ch 11 Eq. 11.19 [UJ 
PDB Ch 11 Eq. 11.26 [UJ 



Uncertainty expanded. [UJ 
PDB Sec 10.6 [IT] 



VHig S s(M Z ) 

( q-^(Mz), a B (M z ), sin 2 6 W {M Z 



246.221(20) GeV 
( 127.904(19), 0.1216(17), 0.23122(15)) 



Table 3.1: Low-energy observables. Masses in lower-case m are MS running masses. Capital M 
indicates pole mass. The light quark's (u,d,s) mass are specified at a scale fj,L = 2 GeV. Vckm 
are the Standard Model's best fit values. 



We first consider the determination of the fundamental mass parameters at the GUT scale 
in order simply to compare to GUT predictions. The starting point for the light-quark masses 



at low scale is given by the % fit to the data of Table 3.1 



m u {ix L ) = 2.7 ± 0.5 MeV m d (fi L ) = 5.3 ± 0.5 MeV m s (/i L ) = 103 ± 12 MeV. 



(3.1) 



Using these as input we determine the values of the mass parameters at the GUT scale for 
various choices of tan/3 but not including possible tan/3 enhanced threshold corrections. We do 
this using numerical solutions to the RG equations. The one-loop and two-loop RG equations 
for the gauge couplings and the Yukawa couplings in the Standard Model and in the MSSM 
that we use in this study come from a number of sources [110J [113 [40j[ll4] and are detailed in 



Appendix [AJ The results are given in the first five columns of Table 3.2 These can readily be 
compared to expectations in various Grand Unified models. The classic prediction of SU (5) with 
third generation down-quark and charged-lepton masses given by the coupling B 5/.10/.5j^]is 

1 5/, 10/ refer to the SU(5) representations making up a family of quarks and leptons while 5h is a five 
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mb(Mx)/m T (Mx) = 1 [103J. This ratio is given in Table 3.2 where it may be seen that the 
value agrees at a special low tan (3 value but for large tan (5 it is some 25% smaller than the 
GUT prediction^] A similar relation between the strange quark and the muon is untenable 
and to describe the masses consistently in SU{5) Georgi and Jarlskog proposed that the 
second generation masses should come instead from the coupling C 5j.10j.45j/ leading instead 



to the relation 3m s (Mx) /m^(Mx) = 1- As may be seen from Table 3.2 in all cases this ratio is 
approximately 0.69(8). The prediction of Georgi and Jarlskog for the lightest generation masses 
follows from the relation Det(M d ) / Det(M l ) = 1. This results from the form of their mass 
matrix which is given bjj^] 

M d = I A C I , M l = I A -3C J (3.2) 





in which there is a (1, 1) texture zerdjand the determinant is given by the product of the (3, 3), 
(1,2) and (2, 1) elements. If the (1,2) and (2, 1) elements are also given by 5y.l0y.5# couplings 
they will be the same in the down-quark and charged-lepton mass matrices giving rise to the 



equality of the determinants. The form of Eq(3.2) may be arranged by imposing additional 



continuous or discrete symmetries. One may see from Table 3.2 that the actual value of the 
ratio of the determinants is quite far from unity disagreeing with the Georgi Jarlskog relation. 

In summary the latest data on fermion masses, while qualitatively in agreement with the 
simple GUT relations, has significant quantitative discrepancies. However the analysis has 
not, so far, included the SUSY threshold corrections which substantially affect the GUT mass 
relations at large tan/? [116 . A catalog of the full SUSY threshold corrections is given in 
|117j . The particular finite SUSY thresholds discussed in this letter do not decouple as the 
superpartners become massive. We follow the approximation described in Blazek, Raby, and 
Pokorski (BRP) for threshold corrections to the CKM elements and down-like mass eigenstates 
|118j . The finite threshold corrections to Y e and Y u and are generally about 3% or smaller 

5Y U , 5Y d < 0.03 (3.3) 

and will be neglected in our study. The logarithmic threshold corrections are approximated by 
using the Standard- Model RG equations from Mz to an effective SUSY scale Ms- 



dimensional representation of Higgs scalars. 

2 We'd like to thank Ilja Dorsner for pointing out that the tan/3 dependence of mb/m T (Mx) is more flat than 
in previous studies (e.g. ref. [115) ') . This change is mostly due to the higher effective SUSY scale Ms, the higher 
value of a 3 (Mz) found in global standard model fits, and smaller top-quark mass Mt- 

3 The remaining mass matrix elements may be non-zero provided they do not contribute significantly to the 
determinant. 

4 Below we discuss an independent reason for having a (1, 1) texture zero. 
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Parameters 


Input SUSY Parameters 
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Parameters 


Corresponding GUT-Scale Parameters with Propagated Uncertainty 


y\M x ) 
y b (M x ) 
y T (M x ) 

{m u /m c ){M x ) 
{m d /m s ){M x ) 
(m e /m^)(M x ) 
(m c /m t )(M x ) 
(m s /m b )(M x ) 
(m M /m T )(M x ) 

A{M X ) 

X(M X ) 

p(M x ) 

fj(M x ) 
J(M X ) x 10" 5 


0.0113+™? 02 

0.0114(3) 

0.0027(6) 

0.051(7) 

0.0048(2) 

00 09+ 001 
u.uuuy_ ooo6 

0.014(4) 

0.059(2) 

56 +0 - 34 

0.227(1) 

0.22(6) 

0.33(4) 

1 4+2-2 


0.48(2) 
0.051(2) 
0.070(3) 
0.0027(6) 
0.051(7) 
0.0048(2) 
0.0025(2) 
0.019(2) 
0.059(2) 
0.77(2) 
0.227(1) 
0.22(6) 
0.33(4) 
2.6(4) 


0.49(2) 
0.23(1) 
0.32(2) 
0.0027(6) 
0.051(7) 
0.0048(2) 
0.0024(2) 
0.017(2) 
0.054(2) 
0.75(2) 
0.227(1) 
0.22(6) 
0.33(4) 
2.5(4) 


0.51(3) 
0.37(2) 
0.51(4) 
0.0027(6) 
0.051(7) 
0.0048(2) 
0.0023(2) 
0.016(2) 
0.050(2) 
0.72(2) 
0.227(1) 
0.22(6) 
0.33(4) 
2.3(4) 


0.51(2) 
0.34(3) 
0.34(2) 
0.0026(6) 
0.051(7) 
0.0048(2) 
0.0023(2) 
0.018(2) 
0.054(2) 
0.73(3) 
0.227(1) 
0.22(6) 
0.33(4) 
2.3(4) 


0.51(2) 
0.34(3) 
0.34(2) 
0.0026(6) 
0.051(7) 
0.0048(2) 
0.0023(2) 
0.010(2) 
0.054(2) 
0.46(3) 
0.227(1) 
0.22(6) 
0.33(4) 
1.0(2) 


Parameters 


Comparison with GUT Mass Ratios 


{m b /m T ){M x ) 
(3m s /m M )(M x ) 
(m d /3m e ){M x ) 


0.82(7) 


0.73(3) 
0.69(8) 
0.83(7) 


0.73(3) 
0.69(8) 
0.83(7) 


0.73(4) 
0.69(8) 
0.83(7) 


1.00(4) 

0.9(1) 

1.05(8) 


1.00(4) 
0.6(1) 
0.68(6) 


^7+ - 08 

u - ' -0.26 


0.42(7) 


0.42(7) 


0.42(7) 


0.92(14) 


0.39(7) 



Table 3.2: The mass parameters continued to the GUT-scale M x for various values of tan/3 
and threshold corrections Jtbd- These are calculated with the 2-loop gauge coupling and 2-loop 
Yukawa coupling RG equations assuming an effective SUSY scale Ms = 500 GeV. 



The finite, tan /3-enhanced Y d SUSY threshold corrections are dominated by the a sbottom- 
gluino loop, a stop-higgsino loop, and a stop-chargino loop. Integrating out the SUSY particles 
at a scale Mg leaves the matching condition at that scale for the Standard-Model Yukawa 
couplings: 

5m sch Y uSM = sin/3 Y u (3.4) 
5m sch Y dSM = cos p Uf (l + T d + V^ KM T u V C km) Y d ag U d R (3.5) 
Y eSM = cos/? Y e . (3.6) 

All the parameters on the right-hand side take on their MSSM values in the DR scheme. The 
factor 5m sc h converts the quark running masses from MS to DR scheme. Details about this 
scheme conversion are listed in Appendix |A.3[ The (5 corresponds to the ratio of the two Higgs 
VEVs v u /vd = tan/3. The U matrices decompose the MSSM Yukawa couplings at the scale 
M s : Y u = U^Y&JTfr and Y d = U^Y d U%. The matrices Y£ and Y d are diagonal and 
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correspond to the mass eigenstates divided by the appropriate VEV at the scale Ms- The CKM 
matrix is given by Vckm = U^llf 1 . The left-hand side involves the Standard- Model Yukawa 
couplings. The matrices T u and T d encode the SUSY threshold corrections. 

If the squarks are diagonalized in flavor space by the same rotations that diagonalize the 
quarks, the matrices T u and T d are diagonal: T d = diag(7d, 7^, 75), T u = diag(7„, j u , 74). In 
general the squarks are not diagonalized by the same rotations as the quarks but provided the 
relative mixing angles are reasonably small the corrections to flavour conserving masses, which 
are our primary concern here, will be second order in these mixing angles. We will assume T u 
and T d are diagonal in what follows. 

Approximations for T u and T d based on the mass insertion approximation are found in 
[113] [120] [121] : 

2 a tan {3 _ , o 2 2\ 2 tciii f3 11 . . 

11 ~ VtJ1 T^^Kv" 1 *^ ) ~ y *32^^T (3 - 7) 

9 tclll ^3 r /222\ /\ 

7 U -g 2 M 2 ii^^h{rn Xl ,m X2 ,m il ) ~ (3.8) 

8 2 tan ^ 2 2 2^ 4 2 tan/3/xM 3 

76 * 3 ff3 16^ M3/xl3 H' m 5 2 ' M3 ) ~ 3 53 i6^^T (3 - 9) 

b 

ld * 3 ff3 16^ M3 ^ /3(m ^' m dV M3 ) ~ 3 53 16^ (3 - 10) 



where ^3 is given by 



2 2 2, a^logg + feVlogg + cVlogg 

,6 ,C ) - (a 2_ 6 2 )(6 2_ c 2 )(a 2_ c2) • 

In these expressions q refers to superpartner of q. x J indicate chargino mass eigenstates. \i 
is the coefficient to the H u H d interaction in the superpotential. Mi, M 2 , M3 are the gaugino 
soft breaking terms. A t refers to the soft top-quark trilinear coupling. The mass insertion 
approximation breaks down if there is large mixing between the mass eigenstates of the stop or 
the sbottom. The right-most expressions in Eqs fl3.7|3.9|3~T0l ) assume the relevant squark mass 



eigenstates are nearly degenerate and heavier than M3 and fi. These expressions ( eqs 3.7 - 3.10) 
provide an approximate mapping from a supersymmetric spectra to the 7« parameters through 
which we parameterize the threshold corrections; however, with the exception of Column A of 



Table 3.4 we do not specify a SUSY spectra but directly parameterize the thresholds corrections 
through 7j. 

The separation between 7^ and 7^ is set by the lack of degeneracy of the down-like squarks. 
If the squark masses for the first two generations are not degenerate, then there will be a 
corresponding separation between the (1,1) and (2,2) entries of T d and T u . If the sparticle 
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spectra is designed to have a large A 1 and a light stop, 74 can be enhanced and dominate over 
76. Because the charm Yukawa coupling is so small, the scharm-higgsino loop is negligible, 
and 7 U follows from a chargino squark loop and is also generally small with values around 0.02 
because of the smaller §2 coupling. In our work, we approximate ~ r" x ~ 0. The only 
substantial correction to the first and second generations is given by 7^ [116J. 

As described in BRP, the threshold corrections leave \ V US \ and | V^b/V^l unchanged to a good 
approximation. Threshold corrections in T u do affect the V u b and V c b at the scale Ms giving 

ySM _ yMSSM ySM _ yMSSM 

yMSSM — yMSSM — ^ '* ' u > ' {O.iz.) 

The threshold corrections for the down-quark masses are given approximately by 

m d ±2 m° d (l +Jd + lu)~ 1 
m s ^2 m°(l +7^ + 7n) _1 
m b ^ m°(l +76 + 7*)" 1 

where the superscript denotes the mass without threshold corrections. Not shown are the 
nonlinear effects which arise through the RG equations when the bottom Yukawa coupling 
is changed by threshold effects. These are properly included in our final results obtained by 
numerically solving the RG equations. 

Due to our assumption that the squark masses for the first two generations are degenerate, 
the combination of the GUT relations given by (detM 1 / det M d ) (3m s /m^) 2 (mb/m T ) = 1 is 
unaffected up to nonlinear effects. Thus we cannot simultaneously fit all three GUT relations 
through the threshold corrections. A best fit requires the threshold effects given by 

7fc + 7t -0.22 ±0.02 (3.13) 
ld + lu ~ -0.21 ±0.02. (3.14) 



giving the results shown in the penultimate column of Table 3.2 just consistent with the GUT 
predictions. The question is whether these threshold effects are of a reasonable magnitude and, if 
so, what are the implications for the SUSY spectra which determine the 7i? From Eqs( |3.9|3.10) , 
at tan j3 = 38 we have 



2 



-0.5, — % ~ 1.0 



2 <-/.u, n 

b d 



The current observation of the muon's {g — 2)„ is 3.4 a |122j away from the Standard-Model 
prediction. If SUSY is to explain the observed deviation, one needs tan (5 > 8 [123 and [XM2 > 
|77j . With this sign we must have 11M3 negative and the d, s squarks only lightly split from the b 
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squarks. M3 negative is characteristic of anomaly mediated SUSY breaking[124] and is discussed 
in [125] [126] [121] [127] . Although we have deduced M 3 < from the approximate Eqs( |3.9|3.10] ), 



the correlation persists in the near exact expression found in Eq(23) of Ref |118] . Adjusting to 
different squark splitting can occur in various schemes |128j . However the squark splitting can 
readily be adjusted without spoiling the fit because, up to nonlinear effects, the solution only 



requires the constraints implied by Eq(|3.13|), so we may make 7^ > 7^ and hence make m~ < m ~ 



by allowing for a small positive value for jt. In this case A* must be positive. 

It is of interest also to consider the threshold effects in the case that is positive. This 



is illustrated in the last column of Table 3.2 in which we have reversed the sign of 7^, consistent 
with positive 11M3 , and chosen 7^ ~ 7^ as is expected for similar down squark masses. The 
value of 74 is chosen to keep the equality between nib and m T . One may see that the other GUT 
relations are not satisfied, being driven further away by the threshold corrections. Reducing 
the magnitude of 7;, and jd reduces the discrepancy somewhat but still limited by the deviation 



found in the no-threshold case (the fourth column of Table 3.2). 

At tan (3 near 50 the non-linear effects are large and b — r unification requires 7b + 74 ~ —0.1 
to —0.15. In this case it is possible to have t — b — t unification of the Yukawa couplings. For 
fi > 0, M 3 > 0, the "Just-so" Split-Higgs solution of references [1291 US [131] [[32] can achieve 
this while satisfying both b — > s 7 and (g — 2)^ constraints but only with large 7^ and 74 and a 
large cancellation in jb+7t- In this case, as in the example given above, the threshold corrections 
drive the masses further from the mass relations for the first and second generations because 
//M 3 > 0. It is possble to have t — b — r unification with /iM 3 < 0, satisfying the b — > s 7 
and {g — 2)^ constraints in which the GUT predictions for the first and second generation of 
quarks is acceptable. Examples include Non-Universal Gaugino Mediation [133 and AMSB; 
both have some very heavy sparticle masses ( > 4 TeV) |121j . Minimal AMSB with a light 
sparticle spectra( < 1 TeV), while satisfying (g — 2) M and & -» s 7 constraints, requires tan/3 
less than about 30 [77] . 

3.2 Updated fits to Yukawa matrices 

We turn now to the second part of our study in which we update previous fits to the Yukawa 
matrices responsible for quark and lepton masses. As discussed above we choose to work in a 
basis in which the mass matrices are hierarchical with the off-diagonal elements small relative 



to the appropriate combinations of on-diagonal matrix elements defined in Eq(B.l). This is 



the basis we think is most likely to display the structure of the underlying theory, for example 
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Parameter 


2001 RRRV 


Fit AO 


Fit BO 


Fit Al 


Fit Bl 


Fit A2 


Fit B2 


tan (5 


Small 


1.3 


1.3 


38 


38 


38 


38 


a' 


Oil) 














-2.0 


-2.0 


e u 


0.05 


0.030(1) 


0.030(1) 


0.0491(16) 


0.0491(15) 


0.0493(16) 


0.0493(14) 


£d 


0.15(1) 


0.117(4) 


0.117(4) 


0.134(7) 


0.134(7) 


0.132(7) 


0.132(7) 


|6'| 


1.0 


1.75(20) 


1.75(21) 


1.05(12) 


1.05(13) 


1.04(12) 


1.04(13) 


arg(6 ) 


90° 


+ 93(16)° 


- 93(13)° 


+ 91(16)° 


-91(13)° 


+ 93(16)° 


- 93(13)° 


a 


1.31(14) 


2.05(14) 


2.05(14) 


2.16(23) 


2.16(24) 


1.92(21) 


1.92(22) 


b 


1.50(10) 


1.92(14) 


1.92(15) 


1.66(13) 


1.66(13) 


1.70(13) 


1.70(13) 


c 


0.40(2) 


0.85(13) 


2.30(20) 


0.78(15) 


2.12(36) 


0.83(17) 


2.19(38) 


arg(c) 


-24(3)° 


-39(18)° 


-61(14)° 


-43(14)° 


-59(13)° 


-37(25)° 


-60(13)° 



Table 3.3: Results of a x 2 fit °f eqs( 3.15|3.16 ) to the data in Table 3.2 in the absence of threshold 
corrections. We set a' as indicated and set d = d! = d = and / = /' = 1 at fixed values. 



that of a spontaneously broken family symmetry, in which the hierarchical structure is ordered 
by the (small) order parameter breaking the symmetry. With this structure to leading order 
in the ratio of light to heavy quarks the observed masses and mixing angles determine the 
mass matrix elements on and above the diagonal provided the elements below the diagonal are 
not anomalously large. This is the case for matrices that are nearly symmetrical or for nearly 
Hermitian as is the case in models based on an 50(10) GUT. 

For convenience we fit to symmetric Yukawa coupling matrices but, as stressed above, this 
is not a critical assumption as the data is insensitive to the off-diagonal elements below the 
diagonal and the quality of the fit is not changed if, for example, we use Hermitian forms. For 
comparison Appendix [B] gives the observables in terms of Yukawa matrix entries following a 
general hierarchical texture. We parameterize a set of general, symmetric Yukawa matrices as: 



Y U (M X 
Y d (M 




x, 



b'el 


c'el 






a' el 






ce 3 d \ 






a 4 





(3.15) 
(3.16) 



Although not shown, we always choose lepton Yukawa couplings at Mx consistent with the 
low-energy lepton masses. Notice that the / coefficient and are redundant (likewise in Y u ). 
We include / to be able to discuss the phase of the (2,2) term. We write all the entries in terms 
of e so that our coefficients will be 0(1). We will always select our best e parameters such that 

1/1 = 1. 

RRRV noted that all solutions, to leading order in the small expansion parameters, only 
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depend on two phases 4>i and <j>2 given by 



- 4>'f) - (4>b - </>/) 

i>c - <t>a) ~ Ob - 4>f) 



(3.17) 
(3.18) 



where <j) x is the phase of parameter x. For this reason it is sufficient to consider only b' and c as 
complex with all other parameters real. 

As mentioned above the data favours a texture zero in the (1, 1) position. With a symmetric 
form for the mass matrix for the first two families, this leads to the phenomenologically successful 
Gatto Sartori Tonin [134J relation 



V US {M X ) 



be d -\b'\e l ^'e u 



(— )o- J^W* 



m. 



(3.19) 



This relation gives an excellent fit to V us with cfii ~ ±90°, and to preserve it we take d, d! to 
be zero in our fits. As discussed above, in SU(5) this texture zero leads to the GUT relation 
Det(M d )/Det(M l ) = 1 which, with threshold corrections, is in good agreement with experiment. 
In the case that c is small it was shown in RRRV that <j)\ is to a good approximation the CP 
violating phase 5 in the Wolfenstein parameterization. A non-zero c is necessary to avoid the 
relation V u bl^cb = \/ m u/ r nc an d with the improvement in the data, it is now necessary to have 
c larger than was found in RRRV As a result the contribution to CP violation coming from 
4>2 is at least 30%. The sign ambiguity in (p\ gives rise to an ambiguity in c with the positive 



sign corresponding to the larger value of c seen in Tables 3.3 and 3.4 



Table 3.3 shows results from a x 2 fit of Eqs( 3.15|3.16 ) to the data in Table 3.2 in the absence 
of threshold corrections. The error, indicated by the term in brackets, represent the widest 
axis of the la error ellipse in parameter space. The fits labeled 'A' have phases such that we 
have the smaller magnitude solution of |c|, and fits labeled 'B' have phases such that we have 
the larger magnitude solution of |c|. As discussed above, it is not possible unambiguously to 
determine the relative contributions of the off-diagonal elements of the up and down Yukawa 
matrices to the mixing angles. In the fit A2 and B2 we illustrate the uncertainty associated 



with this ambiguity, allowing for 0(1) coefficients a' . In all the examples in Table 3.3 the mass 



ratios, and Wolfenstein parameters are essentially the same as in Table 3.2 



The effects of the large tan/3 threshold corrections are shown in Table 3.4 The threshold 
corrections depend on the details of the SUSY spectrum, and we have displayed the effects 



As shown in ref. [135| . it is possible, in a basis with large off-diagonal entries, to have an Hermitian pattern 
with the (1,1) and (1,3) zero provided one carefully orchestrates cancelations among Y u and Y d parameters. We 
find this approach requires a strange-quark mass near its upper limit. 
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Parameter 


A 


B 


c 


B2 


C2 




tan /3 




30 


38 


38 


38 


38 




7b 




0.20 


-0.22 


+0.22 


-0.22 


+0.22 




7* 




-0.03 





-0.44 





-0.44 




Id 




0.20 


-0.21 


+0.21 


-0.21 


+0.21 




a' 













-2 


-2 


w\ 

\ u 1 
arg(6') 
a 
b 

c 

arg(c) 

VU 

2 33 
vd 
x 33 
ye 
1 33 


0.0495(17) 
0.131(7) 
1.04(12) 
90(12)° 
2.17(24) 
1 69(4 3 s ! 
0.80(16) 
-41(18)° 
0.48(2) 
0.15(1) 
0.23(1) 


0.0483(16) 
0.128(7) 
1.07(12) 
91(12)° 
2.27(26) 
1 73(4 3^ 
0.86(17) 
-42(19)° 
0.51(2) 
0.34(3) 
0.34(2) 


0.0483(18) 
0.102(9) 
1.07(11) 
93(12)° 
2.30(42) 
9 91 f4 R) 
1.09(33) 
-41(14)° 
0.51(2) 
0.34(3) 
0.34(2) 


0.0485(17) 

0. 127(7) 
1.05(12) 
95(12)° 
2.03(24) 
1 74(4 0"! 

1. J 4:1 1U 1 

0.81(17) 
-53(10)° 
0.51(2) 
0.34(3) 
0.34(2) 


0.0485(18) 
0.101(9) 
1.06(10) 
95(12)° 
1.89(35) 

1.10(35) 
-41(12)° 
0.51(2) 
0.34(3) 
0.34(2) 


(m 6 /m T )(M x ) 

(3m./m M )(Afx) 
K/3m e )(Mx) 


0.67(4) 
0.60(3) 
0.71(7) 


1.00(4) 
0.9(1) 
1.04(8) 


1.00(4) 
0.6(1) 
0.68(6) 


1.00(4) 
0.9(1) 
1.04(8) 


1.00(4) 
0.6(1) 
0.68(6) 




dot y"(M x ) 
dot Y C (M X ) 




0.3(1) 


0.92(14) 


0.4(1) 


0.92(14) 


0.4(1) 



Table 3.4: A x 2 fit of Eqs( 3.15|3.16 ) including the SUSY threshold effects parameterized by the 
specified 7j. 

corresponding to a variety of choices for this spectrum. Column A corresponds to a "standard" 
SUGRA fit - the benchmark Snowmass Points and Slopes (SPS) spectra lb of ref([3]). Because 
the spectra SPS lb has large stop and sbottom squark mixing angles, the approximations given 



in Eqns(3.7 3.10) break down, and the value for the correction 7j in Column A need to be 
calculated with the more complete expressions in BRP |118] . In the column A fit and the next 
two fits in columns B and C, we set a' and d to zero. Column B corresponds to the fit given in 



the penultimate column of Table 3.2 which agrees very well with the simple GUT predictions. 
It is characterized by the "anomaly-like" spectrum with M3 negative. Column C examines the 
M3 positive case while maintaining the GUT prediction for the third generation mi, = m T . 
It corresponds to the "Just-so" Split-Higgs solution. In the fits A, B and C the value of the 
parameter a is significantly larger than that found in RRRV. This causes problems for models 
based on non-Abelian family symmetries, and it is of interest to try to reduce a by allowing a', b' 
and c' to vary while remaining 0(1) parameters. Doing this for the fits B and C leads to the fits 
B2 and C2 given in Table |3.4| where it may be seen that the extent to which a can be reduced 
is quite limited. Adjusting to this is a challenge for the broken family-symmetry models. 

Although we have included the finite corrections to match the MSSM theory onto the Stan- 
dard Model at an effective SUSY scale Ms = 500 GeV, we have not included finite corrections 
from matching onto a specific GUT model. Precise threshold corrections cannot be rigorously 
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calculated without a specific GUT model. Here we only estimate the order of magnitude of 



corrections to the mass relations in Table 3.2 from matching the MSSM values onto a GUT 



model at the GUT scale. The tan/3 enhanced corrections in Eq(3.7 3.10) arise from soft SUSY 
breaking interactions and are suppressed by factors of Msusy /Mqut m the high-scale match- 
ing. Allowing for 0(1) splitting of the mass ratios of the heavy states, one obtains corrections 



to y /y T (likewise for the lighter generations) of 0(j£^) from the X and Y gauge bosons and 

v 2 

0{ , 4 A 2 ) from colored Higgs states. Because we have a different Higgs representations for differ- 



■ (4tt) 2 

ent generations, these threshold correction will be different for correcting the Srus/m^ relation 



than the mb/m T relation. These factors can be enhanced in the case there are multiple Higgs 
representation. For an SU(5) SUSY GUT these corrections are of the order of 2%. Planck 
scale suppressed operators can also induce corrections to both the unification scale |136| and 
may have significant effects on the masses of the lighter generations [137J. In the case that the 
Yukawa texture is given by a broken family symmetry in terms of an expansion parameter e, 
one expects model dependent corrections of order e which may be significant. 

Chapter Summary 

In summary, in the light of the significant improvement in the measurement of fermion mass 
parameters, we have analyzed the possibility that the fermion mass structure results from an 
underlying supersymmetric GUT at a very high-scale mirroring the unification found for the 
gauge couplings. Use of the RG equations to continue the mass parameters to the GUT scale 
shows that, although qualitatively in agreement with the GUT predictions coming from simple 
Higgs structures, there is a small quantitative discrepancy. We have shown that these discrep- 
ancies may be eliminated by finite radiative threshold corrections involving the supersymmetric 
partners of the Standard-Model states. The required magnitude of these corrections is what is 
expected at large tan /3, and the form needed corresponds to a supersymmetric spectrum in which 
the gluino mass is negative with the opposite sign to the wino mass. We have also performed 
a fit to the recent data to extract the underlying Yukawa coupling matrices for the quarks and 
leptons. This is done in the basis in which the mass matrices are hierarchical in structure with 
the off-diagonal elements small relative to the appropriate combinations of on-diagonal matrix 
elements, the basis most likely to be relevant if the fermion mass structure is due to a spon- 
taneously broken family symmetry. We have explored the effect of SUSY threshold corrections 
for a variety of SUSY spectra. The resulting structure has significant differences from previous 
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fits, and we hope will provide the "data" for developing models of fermion masses such as those 
based on a broken family symmetry. 

Since this work was first published, its conclusions have been confirmed by studies of other 
research groups [138 . The updated fits to the Yukawa textures and viability of the Georgi- 
Jarlskog relations have been used in numerous string theory and family symmetry models. 
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Chapter 4 



Mass Determination Toolbox at 
Hadron Colliders 

Chapter Overview 

In the previous chapter, we presented arguments that predict the sign of the gluino mass relative 
to the wino mass and relationships between these masses needed to satisfy classic Georgi- Jarlskog 
mass relationships at the GUT scale. 

This chapter and the remaining chapters discuss mostly model-independent experimental 
techniques to determine the mass of a pair-produced dark-matter particle. As a test case, we 
take the dark-matter to be the lightest supersymmetric particle (LSP) which we assume is the 
neutralino xt\ Determining the mass of the dark-matter particle is a necessary prerequisite to 
determining the entire mass spectrum of the new particle states and to determining the sign 
of the gluino's mass (given that supersymmetry is the correct model) predicted in the previous 
chapter. 

This chapter reviews challenges of hadron-collider mass-determination and mostly model- 
independent techniques to address these challenges existent in the literature. In the subsequent 
chapters of this thesis, we will improve on these techniques and develop new techniques that can 
perform precision measurements of the dark-matter's mass with only a few assumptions about 
the underlying model. 

4.1 Mass Determination Challenges at Hadron Colliders 

Measuring all the masses and phases associated with the predictions in Chapter [3] or the pre- 
dictions of any of the many competing models depends on the successful resolution of several 
challenges present at any hadron collider. First to avoid selection bias and because of the large 

1 The four neutralinos are superpositions of the two Higgsinos, the bino, and the neutral wino (hZ, h^, B, W 3 ). 
They are numbered from 1 having the smallest mass to 4 having the largest mass. 
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number of possible models, we would prefer kinematic model-independent techniques to mea- 
sure parameters like the mass instead of model-dependent techniques. Kinematic techniques are 
made complicated because of the possibility of producing dark matter in the collider which would 
lead to new sources (beyond neutrinos) of missing transverse momentum. Kinematic techniques 
are also made complex because the reference frame and center-of-mass energy of the parton 
collision is only known statistically^} Last our particle detector has fundamental limitations on 
the shortest track that can be observed leading to combinatoric ambiguities identifying decay 
products from new particle states that are created and decay at the "effective" vertex of the 
collision. 

4.1.1 Kinematic versus Model-Dependent Mass Determination 

Any new physics discovered will come with unknown couplings, mixing angles, spins, interac- 
tions, and masses. One approach is to assume a model and fit the model parameters to the 
observed data. However, we have hundreds of distinct anticipated theories possessing tens of 
new particle states with differing spins and couplings each with tens of parameters. Global fits 
to approximately 10 9 events recorded per year are an enormous amount of work, and assume 
that the 'correct model' has been anticipated. 

We would like to measure the particle properties with minimal model assumptions. The 
term 'model-independent' technique is misleading because we are always assuming some model. 
When we say model-independent, we mean that we assume broadly applicable assumptions 
such as special relativity. We also try to ensure that our techniques are largely independent of 
a priori unknown coupling coefficients and model parameters. For this reason our desire for 
model-independence constrains our toolbox to kinematic properties of particle production and 
decay. 

If we are able to determine the mass of the dark-matter particle, other properties follow 
more easily. Knowing the dark-matter particle's mass, the remaining mass spectrum follows 
from kinematic edges. If we know the dark-matter's mass and gluino's mass, then the sign of 
the gluino's mass predicted in Chapter |3]can be determined from the distribution of the invariant 
mass of jet pairs from the decay of gluino pair production [139 . Another avenue to the sign 
of the gluino's mass requires measuring the masses of the gluino and the two stop squarks 
and the measuring the decay width of gluino to these states |140| Appendix B]. In addition, 
2 To the extent we understand the uncertainties in the measured parton distribution functions. 
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the measurement of spin correlations, a study that can contribute to spin determination, has 
historically relied on knowing the masses to reconstruct the event kinematics [141J. 

In short, mass determination is a key part to identifying the underlying theory which lies 
beyond the Standard Model. Newly discovered particles could be Kaluza-Klein (KK) states from 
extra-dimensions, supersymmetric partners of known states, technicolor hadrons, or something 
else that we have not anticipated. Models predict relationships between parameters. Super- 
symmetry relates the couplings of current fermions to the couplings of new bosons, and the 
supersymmetric particle masses reflect the origin of supersymmetry breaking. Masses of KK 
states tell us about the size of the extra dimensions. If these parameters, such as the masses and 
spins, can be measured without assuming the model, then these observations exclude or confirm 
specific models. In general, mass determination of new particle states is central to discerning 
what lies beyond the Standard Model. 

4.1.2 Dark matter: particle pairs carrying missing transverse momentum 

If dark matter is produced at a hadron collider, the likely signature will be missing transverse 
momentum. The lightest supersymmetric particle (LSP) or lightest Kaluza-Klein particle (LKP) 
is expected to be neutral, stable, and invisible to the particle detectors. The astrophysical 
dark matter appears to be stable. Whatever symmetry makes dark matter stable and / or 
distinguishes superpartners from their Standard-Model cousins will likely also require that dark 
matter particles be pair produced. Therefore, events with supersymmetric new particle states 
are expected to end in two LSP's leaving the detector unnoticed. Events with Kaluza-Klien 
particles are expected to end with two LKP's leaving the detector unnoticed. The presence of 
two invisible particles prevents complete reconstruction of the kinematics of any single event 
and leads to missing transverse momentum in the events. 

4.1.3 Reference frame unknown due to colliding hadrons 

At a hadron collider, because the partons colliding each carry an unknown fraction of the 
hadron's momentum we do not know the reference frame of the collision. Although the individual 
parton momentum is unknown, the statistical distribution of the parton's momentum can be 
measured from deep inelastic scattering |142j . The measured parton distribution functions 
(PDFs) f^(x,Q) give the distribution of parton i within the hadron H carrying x fraction of 
the hadron's momentum when probed by a space-like momentum probe g M with Q 2 = —q 2 . Q 
is also the factorization scale at which one cuts off collinear divergences. The dominant parton 
types i are u, d, s, their antiparticles u, d, s, and gluons g. 
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The events produced in the collider follow from a convolution of the cross section over the par- 
ton distribution functions . If the two colliding protons have 4-momentum p\ = ^slhc{ 1 > 0, 0, 1) /2 
and p2 = \/slhc(1,0, 0, — 1)/2, then for example the u and u quarks colliding would have 4- 
momentum x u p\ and x^P2- The spatial momentum of the collision along the beam axis is given 
by (x u — Xu)y/SLHC /2, and the center-of-mass energy of the parton collision is yj x u x^slhc- 
Because x u and xu are only known statistically, any individual collision has an unknown center- 
of-mass energy and an unknown momentum along the beam axis. 

LHC processes are calculated by convolutions over the PDFs as in 

a = J dx u dx n J%(X U , Q)fl(%u, Q)v(Xu X n S L CH, Q)uu^final (4.1) 

where a(s, Q) U u—> final is the total cross section for a given parton collision center-of-mass squared 
s, factorization scale Q, and ^/slhc is the center-of-mass collision energy of the protons colliding 
at the LHC. The region of integration of x u and x^ is based on the kinematically allowable 
regions. The factorization scale Q is chosen to minimize the size of the logarithms associated 
with regulating collinear divergences in the cross section. See Ref |143] for a recent review on 
choosing the factorization scale. 

Calculating these input cross sections even at tree level is an arduous process that has 
been largely automated. Performing these convolutions over experimentally determined PDFs 
adds even further difficulties. Monte- Carlo generators typically perform these calculations. The 
calculations in subsequent chapters made use of Mad Graph and MadEvent [144J, HERWIG[17 , 
CompHEP [US]. 

Despite this automation, the author has found it useful to reproduce these tools so that we 
can deduce and control what aspects of the observed events are caused by what assumptions. 
For this purpose, parton distributions can be downloaded^] from [146 . Analytic cross sections 
for processes can be produced by CompHEP [145 . 

These review studies have led to an appreciation of uncertainties in the parton distribu- 
tion functions, uncertainties about the correct Q to use^J and uncertainties in the process of 
hadronization of the outgoing partons. All these uncertainties contribute to significant un- 
certainties in the background of which new physics must be discovered, and to uncertainties 
underling model-dependent, cross-section-dependent determinations of particle masses. 

3 One should also look for the Les Houches Accord PDF Interface (LHAPDFs) interfaces to interface PDFs to 
different applications (http://projects.hepforge.org/lhapdf/) . 

4 Ref 1 140] suggests using Q w 2Msusy or whatever energy scale is relevant to the particles being produced. 
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4.1.4 Detector Limitations 

Particle detectors have finite capabilities which limit what we can learn. We will discuss the 
effects of finite energy resolution later in this thesis. Here I focus on the lack of the traditional 
tracks one sees in event representations. 

The early history of particle physics shows beautiful bubble chamber tracks of particles. 
Studies of Pions and Kaons relied on tracks left in bubble chambers or modern digitized detectors. 
In the case of the K$ the lifetime is r = 0.9 x 10" 10 seconds or cr = 2.7 cm, or equivalently 
the decay width is T = 7 x 10 -6 eV. A particle with a width of T = 1 eV has a lifetime of 
r = 6 x 10~ 16 seconds with cr = 0.2 /mi. The more massive known states of , Z° bosons and 
the top quark have widths of T\y = 2.141 GeV, V z = 2.50 GeV, and Tt ~ 1.5 GeV. These decay 
widths give these states tracks with cr ~ 0.1 fm or about 1/10 the size of an atomic nuclei. 

Most supersymmetric states, depending on the model, have decay widths between around 1 
MeV and 2 GeV^J These observations, while trivial, make us realize that new physics discoveries 
will most likely need to deduce properties of the new states from their decay products only. The 
fact that all events effectively occur at the origin creates combinatoric problems because we 
cannot know the order in which the visible particles came off of a decay chain or even to which 
decay chain they belong. 

4.2 Invariant Mass Edge Techniques 

We now turn to existent tools to overcome the challenges inherent in mass determination at 
hadron colliders in the presence of missing transverse momentum. 

4.2.1 Invariant Mass Edges 

Lorentz invariant distributions are optimal for making a distribution independent of the unknown 



frame of reference. Fig |4.1| shows a two-body decay and a three-body decay. The decay products 
from Y in both cases involve a dark-matter particle N and two visible states with 4-momenta 
a± and a.% we will abbreviate a = a\ + ai- The Lorentz invariant mass of these two visible 
states is defined as 

m\ 2 = a 2 = (ai + a 2 ) 2 . (4.2) 

Assuming no spin correlation between a± and Q2, the distribution for the two-body decay looks 
more like a right-triangle whereas the three-body decay case looks softer. Spin correlations 



5 This is not a universal statement. Some viable models have long lived charged states (cr >> 1 cm) [1471 114£?l 
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Figure 4.1: (Left:) Two body decay and it's associated mu distribution. (Right:) Three body 
decay and it's associated mu distribution. 

between these states may change the shape and might be used to measure spin in some cases 
[150] . The distribution shape can also be affected by competing processes in the three-body 
decay. For example, the degree of interference with the slepton can push the peak to a smaller 
value[151j. 

The end-point of the distribution gives information about the masses of the new particle 
states. In the two-body decay case the endpoint gives 

, _ {Ml-M\){M\-Ml) 



maxm 



12 



(4.3) 



and in the three-body decay case 



max mi2 = My — Mpf. 



(4.4) 



These edges in invariant mass combinations provide information about mass differences or mass 
squared differences but not about the mass-scale itself. Measuring these edges requires some 
model of the events creating the distribution in order to simulate the effect of the detectors energy 
resolution. For this reason, even end-point techniques, although mostly model-independent, 
still require some estimate of the distribution near the end point to get an accurate end-point 
measurement. Radiative corrections also play a role in shifting this endpoint slightly [152] . 
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These techniques have been studied to select a model if one assumes a particular set of 
starting models [153 . They have also been used with Bayesian techniques to measure the mass 
difference between slepton states |154j . 

4.2.2 Constraints from Cascade Decays 

If there are many new particle states produced at the collider, then we can use edges from 
cascade decays between these new states to provide constraints between the masses. Given 
some luck regarding the relative masses, these constraints may be inverted to obtain the mass 



scale. Fig 4.2 shows a cascade decay from Z to Y to X ending in N and visible particle momenta 
a>i, a 2 , «3- 

There are four unknown masses and potentially four linearly-independent endpoints. In 



addition to Eq(4.3), we also have [155, 153 



maxm§ 2 = (M§ - M$)(M$ - M x )/M$, 



(4.5) 



and 



maxm 



123 



(M 2 -M 2 ){M 2 -M 2 N )/M 2 
(M 2 -M 2 X )(M 2 X -M 2 N )/M 2 X 



iSM 2 < M N M Z 
iEM N M z < M 2 X 



(M|M| - M$M%)(M$ - M\) /(M^M X ) ffiM x M z < M N M^ 



(4.6) 



^M Z -M N ) 2 
The fourth endpoint is 



otherwise. 



maxm? 3 = (Mf - M$){M 2 X - M 2 N )/M 2 X . 



(4.7) 



Depending on whether we can distinguish visible particles 1, 2 and 3 from each other, some 
of these endpoints may be obscured by combinatorics problems. If the visible particles and 
the masses are such that these four can be disentangled from combinatoric problems, and if 
the masses are such that these four end-points provide independent constraints, and if there 
is sufficient statistics to measure the endpoints well-enough, then we can solve for the mass 
spectrum. 

This technique has been studied in |156| I155| 1153] and subsequently by many others. In 
some cases, there is more than one solution to the equations. This degeneracy can be lifted by 
using the shape of these distributions |157] I158[ I159j . Using the Supersymmetric benchmark 
point SPS la, this technique has been shown to determine the LSP mass to ±3.4 GeV with 
about 500 thousand events from 300 fb _1 . They are able to determine the mass difference to 



vm^o-aUo = 0.2 GeV and a Mr 



-IvUo 



0.3 GeV [157 . 
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Figure 4.2: Cascade decay from Z to Y to X ending in ./V and visible particle momenta ai, 02, 
as. 
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Figure 4.3: Events with the new state Y is pair produced and in which each Y decays through 
a two-body decay to a massive new state X and a visible state 1 and then where X subsequently 
decays to a massive state N invisible to the detector and visible particles 2. All previous decay 
products are grouped into the upstream transverse momentum, k. 

4.3 Mass Shell Techniques 



What if there are not as many new-states accessible at the collider, or what if the masses arrange 
themselves to not enable a solution to the four invariant mass edges listed above? There is also 
a series of approaches called Mass Shell Techniques (MST^ where assumption are used about 
the topology and on-shell conditions to solve for the unknown masses. 

An MST was used in the early 1990s by Goldstein, Dalitz (GD) [161 , and Kondo, Chika- 
matsu, Kim (KCK)|162 as a suggested method to measure the mass of the top quark. The 



top-quark pair production and decay has the topology shown in Fig 4.3 Here Y = t or t 



X 



[1) = b, (2) = e and N = v e . In each event, one measures ai2, Pi,2, and the 



6 The title MST is suggested in Ref. [160 
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missing transverse momentum = (p + q)t- The resulting mass shell and missing momentum 
equations are 

M 2 N = p 2 =q 2 , (4.8) 

M\ = (p + a 2 ) 2 = (q + (3 2 ) 2 , (4.9) 

My = {p + a 1 + a 2 ) 2 = {q + Pi + fc) 2 , (4.10) 

f t = (p + «)r- (4.11) 

If we assume a mass for My, then we have 8 unknowns (the four components of each p and q) 
and 8 equations. The equations can be reduced to the intersection of two ellipses which have 
between zero and four distinct solutions. In GD and KCK, they supplemented these solutions 
with some data from the cross section and used Bayes' theorem to create a likelihood function 
for the unknown mass of the top quark. The approach has a systematic bias [163 [164 that 
must be removed by modeling [165^] 

A modern reinvention of this MST is found in Cheng, Gunion, Han, Marandella and McElrath 
(CHGMM) [166 where they assume a symmetric topology with an on-shell intermediate state 



such as in Fig 4.3 This reinvention assumes Y = x 2 , X = I and N = x°. CHGMM do a test 
of each event's compatibility with all values of (My,Mx,Mn). Using this approach, they find 
that with 300 fb _1 of events that they can determine the LSP mass to ±12 GeV. Unlike GD 
and KCK, CHGMM make no reference to Bayesian techniques. 

Another MST assumes a longer symmetric decay chain and assumes the masses in two events 
must be equal |167] . This results in enough equations to equal the number of unknowns to solve 
for the masses directly. For SPS la after 300 fb -1 , they reach ±2.8 GeV using 700 events 
satisfying their cuts but have a 2.5 GeV systematic bias that needs modeling to remove. 

Finally there is a suggestion for hybrid MST where one combines on-shell conditions with 
the information from the many edges in cascade decays [168] . The M 2 c variable introduced in 
Chapters [6] is a simple example of such a hybrid technique that predated their suggestion. 

4.4 Transverse Mass Techniques 
4.4.1 M T and Measuring M w 

The invariant mass distributions are Lorentz invariant. In events with one invisible particle, like 
the neutrino, leaving the detector unnoticed, we only know the neutrino's transverse momentum 
inferred from the missing transverse momentum. If a parent particle X decays into a visible 



7 Details of this bias are discussed more in Section 8.2.1 
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particle with observed 4-momentum a\ and an invisible particle p, then what progress can 
be made in mass determination? There is a class of techniques which, although not Lorentz 
invariant, is invariant with respect to longitudinal boosts along the beam line. The 4-momentum 
Pn = (Po,Px,Py,Pz) is cast in new coordinates 

-Kp) = |ln^ (4.12) 
2 Po-Pz 



E T (p) = ^ml + pl+pl (4.13) 
p M = (E T coshr],p x ,p y ,E T smhri) (4.14) 



where m 2 = p 2 . The mass of X in the decay X{a\ + p) —* + l(ai) in Fig. 



4.1 



is 



M|- = (p + ai) 2 = M^ + M 1 2 + 2p-«i (4.15) 
= M 2 N + Ml + 2E T (p)E T ( ai ) cosh(?7(p) - 77(01)) - 2p T ■ (a x ) T . (4.16) 

We observe the 4-momentum a± in the detector. From the missing transverse momentum, we 
can deduce the transverse components p x and p z . The components p and p z remain unknown. 
If we know the mass of N, we can bound the unknown mass of X from below by minimizing 
with respect to the unknown rapidity r]{p). The minimum 

M\ > M T (p, ai) = M 2 N + Ml + 2E T (p)E T ( ai ) - 2p T ■ (ai) T (4.17) 

is at r]{p) = r]{a\) and is guaranteed to be less than the true mass of X. This lower bound is 
called the transverse mass Mt- 

This technique is used to measure the boson's mass with the identification X = W^, N = 
v e , and (1) = e^ 1 . The first observation of the was reported in 1983 at UA1 Collaboration at 
Super Proton Synchrotron (SPS) collider at CERN[48 . They used M T to bound M w > 73 GeV 
at the 90% confidence level. Assuming the Standard Model and fitting the data to the model 
with Mw as a free parameter gave the UA1 collaboration the measurement Mw = 81 ± 5 GeV, 
which agreed with the prediction described by Llewellyn Smith and Wheater [37]. Mt is still 
used even for the more recent 2002 DO-Collaboration model-independent measurement of My/ = 



80.483 ±0.084 GeV [169J. 

4.4.2 The Stransverse Mass Variable M T2 

If a hadron collider pair produces dark-matter particles, then there are two sources of miss- 
ing transverse momentum, each with the same unknown mass. The 'stransverse mass' 771^2 
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introduced by Lester and Summers p21 [13] adapts the transverse mass mr to this task. Mt2 
techniques have become widely used with Refs. |12| [T3] having more than 45 citations. 

The stransverse mass is used to determine the mass difference between a parent particle and 
a dark-matter candidate particle given an assumed mass for the dark-matter candidate based on 



a topology similar to Figs 4.4 or 4.3 The variable vtlt2 accepts three inputs: xn (an assumed 
mass of the two particles carrying away missing transverse momenta), a and (3 (the visible 
momenta of each branch), and fx = (p j i~q)t (the missing transverse momenta two- vector). We 
can define rriT2 m terms of the transverse mass of each branch where we minimize the maximum 
of the two transverse masses over the unknown split between p and q of the overall missing 
transverse momenta: 

M% 2 ( X N,a,p,P T ) = min [max {M$.(a,p), M|(A q)}] . (4.18) 

p T +q T =f> T 

In this expression xn is the assumed mass of N, a and are the four momenta of the vis- 
ible particles in the two branches, the transverse mass is given by M^(a,p) = rr? a + x% + 
2{Et{p)Et{oi) —pr • olt) and the transverse energy Et{oi) = \jVr + X% 1S determined from the 



transverse momentum of p and the assumed mass of p. The My here is identical to Eq(4.17) 
with Mtv replaced by the assumed mass xn- An analytic formula for the case with no transverse 
upstream momentum Ut can be found in the appendix of |170j . For each event, the quantity 
Mt2(xn = mNi cti + «2 ; Pi + /?2, Pt) gives the smallest mass for the parent particle compatible 
with the event's kinematics. Under ideal assumptions, the mass of the parent particle Y is 
given by the end-point of the distribution of this m^2 parameter over a large number of events 



like Figs. |4.4| or 4.3. Because a priori we do not know M^v, we need some other mechanism to 
determine Mat. We use x to distinguish assumed values of the masses (xy, Xx, Xn) from the 
true values for the masses (my, mx, rn^). Because of this dependence on the unknown mass, 
we should think of max tot2 as providing a relationship or constraint between the mass of Y 
and the mass of N. This forms a surface in the (xy, Xx, Xn) space on which the true mass will 
lie. We express this relationship as Xy(xn) 

In addition to invariance with respect to longitudinal boosts, Mt2, m the limit where Ut = 0, 
is also invariant with respect to back-to-back boosts of the parent particle pair in the transverse 
plane |171j . 

Mt2 is an ideal tool for top-mass determination at the LHC |172j . Mt2 would apply to events 
where both pair-produced top quarks decay to a b quark and a and where in both branches 



8 In principle this surface would be considered a function of Xy{xx,Xn), but rriT2 makes no reference to the 
mass of X and the resulting constraints are therefore independent of any assumed value of the mass of X. 
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Figure 4.4: Events with the new state Y is pair produced and in which each Y decays through 
a three-body decay to a massive state N invisible to the detector and visible particles 1, 2, 3, 
and 4. All previous decay products are grouped into the upstream transverse momentum, k. 

the decays leptonically to I and V[. To use Mt2 to find the mass of the parent particle, a 
value must be assumed for xn = Mm- For a v we approximate xn = 0. However, even in models 
of new physics where new invisible particles are nearly massless (like the gravitino studied in 
|173j ). we would rather not just assume the mass of the lightest supersymmetric particle (LSP), 
which is needed as an input to the traditional analysis, without measuring it in some model 
independent way. 

4.4.3 Max M T2 Kink Technique 

Mt2 assumes a mass xn for N. Is there any way that one can determine from Mt2 alone? 
At the same time as the author was completing work on M<ic described in Chapter [6j Cho, 
Choi Kim, Park (CCKP) [174) 1171] published an observation on how the mass could be deduced 
from Mt2 alone. If Y is pair produced and each undergoes a three-body decajflto N as shown 



in Fig. 4.4, then a 'kink' in the max m^PS US], will occur at a position which indicates the 



invisible particle mass. This 'kink' is corroborated in Refs. |176| 1177] . Fig 4.5 shows the 'kink' 
in max m,T2 in & n idealized simulation where My = 800 GeV and Mn = 100 GeV. In the CCKP 
example, they study a mAMSB senario where a gluino decays via a three-body decay to the 
LSP and quarks. In this case, they determined the LSP's mass to ±1.7 GeV. The difficulty in 
mass determination with this technique scales as (M+ /M_) 2 = (My + Mjy) 2 / (My — M^) 2 . For 
CCKP's example M + /M_ = 1.3; we will consider examples where M+/M_ w 3. The (M + /M_) 2 
scaling behavior follows by propagating error in M_ to position of the intersection of the two 
curves that form the kink. 



9 The presence of a > 3-body decay is a sufficient but not necessary condition. Two-body decays can also 
display kinks [1751 1176| provided the decaying particles have sufficiently large transverse boosts. 
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Figure 4.5: The maxM^Cx) shows a kink at the true Mjy and My. For this simulation, 
my = 200 GeV, and m N = 99 GeV. 

This kink is quantified by the constrained mass variable mic [15 j [16J that is the subject of 
later chapters of this thesis. We will be applying it to a case about 5 or 6 times more difficult 
because M + /M_ ~ 3 for the neutralino studies on which we choose to focus. 

Chapter Summary 

In this chapter, we have introduced the need for making mass measurements in a model inde- 
pendent manner of collider-produced dark matter. We have outlined many approaches based 
on edges of invariant mass distributions, mass shell techniques, and transverse mass techniques. 
Each technique has distinct domains of validity and requirements. In the remaining chapters 
of this thesis, we introduce and test new model-independent techniques to determine mass of 
dark-matter particles. Which approach turns out to be best is likely to depend on what scenario 
nature hands us, since the various techniques involve different assumptions. Having different 
approaches also offers the advantage of providing a system of redundant checks. 
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Chapter 5 

Using with Cascade decays 



Chapter Overview 

Recently, a paper by Tovey fTTB] introduced a new variable, Mct, with the powerful concept 
of extracting mass information from intermediate stages of a symmetric decay chain. In this 
chapter, we compare Mt2 with Mct- The variable Mct is new but shares many similarities 
and differences with Mt2- We briefly define Mqt and explain when it gives identical results to 
Mt2 and when it gives different results. We comment on benefits of each variable in its intended 
applications. We find that for massless visible particles Mct equals Mt2 m a particular limit, 
but Mt2 has better properties when there is initial state radiation (ISR) or upstream transverse 
momentum (utmJ) We argue that M?2 is a more powerful tool for extracting mass information 
from intermediate stages of a symmetric decay chain. This chapter is based on work first 
published by the author in Ref. [14J. 



5.1 Mt2 and Mct in this context 

Both Mt2 and Mct assume a pair-produced new-particle state followed by each branch decaying 
symmetrically to visible states and dark-matter candidates which escape detection and appear 



as missing transverse momentum. Fig 4.3 is the simplest example on-which we can meaningfully 
compare the two kinematic quantities. The figure shows two partons colliding and producing 
some observed initial state radiation (ISR) or upstream transverse momentum (UTM) with four 
momenta k and an on-shell, pair-produced new state Y. On each branch, Y decays to on-shell 
states X and v i with masses Mx and m Vl , and X then decays to on-shell states iV and V2 with 
masses Mjv and m V2 . The four-momenta of v±, V2 and iV are respectively ai, 02 and p on one 



1 ISR and UTM both indicate kr 7^ 0. We take ISR to mean the case when kr is small compared to the energy 
scales involved in the collision. 
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branch and fli, P2 and q in the other branch. The missing transverse momenta fx is given by 
the transverse part of p + q. 

Tovey |178] recently defined a new variable Mct which has many similarities to Mx2- The 
variable is defined as 

M&riauPi) = (£T(ai) + £ T (/3i)) 2 -(aiT-/3i T ) 2 . (5.1) 

Tovey's goal is to identify another constraint between masses in the decay chain. He observes 
that in the rest frame of Y the momentum of the back-to-back decay products X and v\ is given 
by 

(uf M m \\* (My - K + M X ) 2 )(M* - (m Vl - M x ) 2 ) 
{h{M Y ,Mx,m Vl )) = ^2 (5.2) 

where fc* is the two-body excess-momentum parameter (2BEMP) [^] In the absence of transverse 
ISR (&t = 0) and if the visible particles are effectively massless (m Vl = 0), Tovey observes that 
maxMcr(ai,|3i) is given by 2k*; this provides an equation of constraint between My and Mx- 
Tovey observes that if we could do this analysis at various stages along the symmetric decay 
chain, then all the masses in principle could be determined. 

The big advantage of Mct is m its computational simplicity. It is a simple one line formula 
to evaluate. Also, Mct is intended to only be calculated once per event instead of at a variety 
of choices of the hypothetical LSP mass \- I n contrast, Mt2 is a more computationally intensive 
parameter to compute; but this is aided by the use of a common repository of community tested 
C++ libraries found at |179j . 

How are these two variables similar? Both Mct and Mt2, m the limit of kx = 0, are 
invariant under back-to-back boosts of the parent particles momenta |171] . The variable Mct 
equals Mt2 m the special case where \ = 0, and the visible particles are massless (a\ = fi\ = 0), 
and there is no transverse ISR or UTM {kx = 0) 

Mcr(ctl,Pl) = M T2 (x = 0,a 1 ,p 1 ,f' T = (p + q + a2 + p2)T) if <4 = Pi = °- ( 5 - 3 ) 
= 2(a 1T -Pi T + |«it| IAtD- (5-4) 

The Mct side of the equation is straight forward. The Mt2 side of the expression can be derived 
analytically using the formula for Mt2 given in [T70j; we also show a short proof in Appendix 
[C| Eq(5.3) uses a Mt2 m a unconventional way; we group the observed momenta of the second 



decay products into the missing transverse momenta. In this limit, both share an endpoint of 



2 Tovey refers to this as the 2-body mass parameter Mi. We feel calling this a mass is a bit misleading so we 
are suggesting 2BEMP. 
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2k* = (Mf- - M 2 X )IM Y . To the best of our knowledge, this endpoint was first pointed out by 
CCKP pH] |5 We find it surprising that a physical relationship between the masses follows 



from Mt2 evaluated at a non physical x- I n the presence of ISR or UTM, Eq(5.3) is no longer an 



equality. Furthermore in the presence of ISR or UTM, the end point of the distribution given by 



either side of Eq(5.3) exceeds 2k*. In both cases, we will need to solve a combinatoric problem 
of matching visible particles to their decay order and branch of the event which we leave for 
future research. 

In the case where the visible particle v\ is massive, the two parameters give different end- 
points 

maxM CT (ai,A) = ^ + and (5.5) 

aly aly 



maxM T2 {x = 0,a 1 ,/3i,f> T = ( P + q+a2 + l3 2 )T) = J m 2 Vl + 2(k* + K ^Jk^ + m^J 



(5.6) 



where k* is given by Eq(5.2). Unfortunately, there is no new information about the masses in 



these two endpoints. If we solve Eq(5.5) for Mx and substitute this into Eq(5.6) and (5.2), all 
dependence on My is eliminated. 

5.2 Application to symmetric decay chains 

Tovey's idea of analyzing the different steps in a symmetric decay chain to extract the masses is 
powerful. Up to now, we have been analyzing both variables in terms of the first decay products 
of Y. This restriction is because Mct requires no transverse ISR to give a meaningful endpoint. 
If we were to try and use 02 and P2 to find a relationship between Mx and Mx , then we would 
need to consider the transverse UTM to be {k + a± + Pi)t, which is unlikely to be zero. 

We suggest Mt2 is a better variable with which to implement Tovey's idea of analyzing the 
different steps in a symmetric decay chain because of its ISR properties. With and without ISR, 
Mt2 's endpoint gives the correct mass of the parent particle when we assume the correct value of 
the missing-energy-particle's mass Q For this reason, max M?2 gives a meaningful relationship 
between masses (My,Mx,Mx) for all three symmetric pairings of the visible particles across 
the two branches. A relationship between My and Mx is given by 

Xy(xx) = maxMr 2 (xx,ai,/3i,fV = (p + q + a 2 + Pi)t)- (5.7) 
3 The endpoint given by CCKP is violated for non-zero ISR at xn < Afjv and \n > Afjv- 

4 In principle we could plot the man Mt2{xx , 011, Pi, Ft = ( Q 2 + P2 + P + ?)t) vs xx as a function of transverse 
ISR and the value of xx at which the end point is constant would give the correct value of Mx; at which point 
the distributions end point would give the correct My. In practice we probably will not have enough statistics of 
ISR events. 
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X (GeV) 



Figure 5.1: Shows constraints from max Mt2 used with different combinations as described 
in eqs( 5.7|5.8|5.9 ) and the max mi2 described in Eq(5.11). Intersection is at the true mass 
(97 GeV, 144 GeV, 181 GeV) shown by sphere. Events include ISR but otherwise ideal condi- 
tions: no background, resolution, or combinatoric error. 



A relationship between Mx and Mx can be found by computing 



Xx(xn) = m&xM T2 {xN,a2,P2,f > T = (p + q)t) 



(5.8) 



where we have grouped a± + Pi with the k as a part of the ISR. A relationship between My and 
Mx can be found by using M?2 hi the traditional manner giving 



Xy(xn) = maxM T2 (x7V,ai + a 2 ,Pi + fh,pT = (p + q)t)- 



(5.9) 



Lastly, we can form a distribution from the invariant mass of the visible particles on each branch 



I'D 



12 



(a?i + a 2 ) 2 or mf 2 = (Pi + P2) 2 ■ The endpoint of this distribution gives a relationship 



between My, Mx, and Mn given by 



2 (M^-M> x )(Mi-M> N ) 

maxm 19 = i — 7, — —. 

12 M\ 



(5.10) 



Solving this expression for My gives the relationship 

Xx((maxmf 2 ) + x\ ~ X%) 



Xy(xn,Xx) 



x\ X% 



(5.11) 



Fig 5.1 shows the constraints from Eqs(5.7 5.8 5.9 5.11) in an ideal simulation using (My = 
181 GeV, M x = 144 GeV, M N = 97 GeV), 1000 events, and massless visible particles, and ISR 
added with an exponential distribution with a mean of 50 GeV. These four surfaces in principle 



intersect at a single point (My, Mx, Mx) given by the sphere in the figure 5.1 Unfortunately. 
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all these surfaces intersect the correct masses at a shallow angles so we have a sizable uncertainty 
along the direction of the sum of the masses and tight constraints in the perpendicular directions. 
In other words, the mass differences are well-determined but not the mass scale. From here one 
could use a shape fitting technique like that described in Chapters [6] - |8j to find a constraint on 
the sum of the masses. Tovey's suggestion for extracting information from these intermediate 
stages of a symmetric cascade chain clearly provides more constraints to isolate the true mass 



than one would find from only using the one constraint of Eq(5.9) as described in |174j . However, 
Tovey's suggestion is more feasible using the Mt2 rather than Mct because the constraint 
surfaces derived from Mt2 intersect the true masses even with UTM. 

Chapter summary 

In summary, we have compared and contrasted Mqt with Mt2- The variable Mct is a special 



case of Mt2 given by Eq(5.3) when ISR can be neglected and when the visible particles are 
massless. In this case, the end-point of this distribution gives 2k*, twice the two-body excess 
momentum parameter (2BEMP). If m Vl 7^ 0, the two distributions have different endpoints but 
no new information about the masses. In the presence of ISR the two functions are not equal; 
both have endpoints that exceed 2k*. Because of it's better properties in the presence of UTM 
or ISR, Mt2 is a better variable for the task of extracting information for each step in the decay 
chain. Extracting this information requires solving combinatoric problems which are beyond the 
scope of this chapter. 
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Chapter 6 

The Variable M^c'- Direct Pair 
Production 



Chapter Overview 

In this chapter, we propose an improved method for hadron-collider mass determination of new 
states that decay to a massive, long-lived state like the LSP in the MSSM. We focus on pair- 
produced new states which undergo three-body decay to a pair of visible particles and the new 
invisible long-lived state. Our approach is to construct a kinematic quantity which enforces 
all known physical constraints on the system. The distribution of this quantity calculated for 
the observed events has an endpoint that determines the mass of the new states. However we 
find it much more efficient to determine the masses by fitting to the entire distribution and not 
just the end point. We consider the application of the method at the LHC for various models 
under the ideal assumption of effectively direct production with minimal ISR and demonstrate 
that the method can determine the masses within about 6 GeV using only 250 events. This 
implies the method is viable even for relatively rare processes at the LHC such as neutralino 
pair production. 

This chapter, which is based on work first published by the author and his supervisor in 
Ref. |15j . concentrates on mass determination involving the production of only two new states. 
Our particular concern is to use the available information as effectively as possible to reduce 
the number of events needed to make an accurate determination of My and Mat. The main 
new ingredient of the method proposed is that it does not rely solely on the events close to the 
kinematic boundary but makes use of all the events. Our method constrains the unobserved 
energy and momentum such that all the kinematical constraints of the process are satisfied in- 



cluding the mass difference, Eq(4.4), which can be accurately measured from the 11 spectrum 



discussed in Sec 14.2.11 This increases the information that events far from the kinematic bound- 
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ary can provide about My and significantly reduces the number of events needed to obtain a 
good measurement of the overall mass scale. We develop the method for the case where Y is 
directly pair produced in the parton collision with minimal ISR and where each Y decays via a 
three-body decay to a on-shell final states N + l + + l~. Its generalization to other processes is 
straightforward and considered later in this thesis^} 



The chapter is structured as follows: In Section 6.1 we introduce the M2C distribution 



whose endpoint gives My, and whose distribution can be fitted away from the endpoint to 



determine My and before we have enough events to saturate the endpoint. Section 6.2 
discusses the relationship between our distribution and the kink in Mt2(x) °f CCKP and how 
this relationship can be used to calculate M2C in a computationally efficient manner. We then 
in 



6.3 discuss symmetries and dependencies of the M2C distribution. Section 6.4 estimates the 
performance for a few SUSY models where we include approximate detector resolution effects 
and where we expect backgrounds to be minimal but where we assume (Jct ~ 0). Finally we 
summarize the chapter's findings. 

6.1 An improved distribution from which to determine My 



We consider the event topology shown in Fig 4.4 The new state Y is pair produced. Each 



branch undergoes a three-body decay to the state N with 4-momentum p (g) and two visible 
particles 1 + 2 (3 + 4) with 4-momentum a (/?). The invariant mass 77142 (w.34) of the particles 
1 + 2 (3 + 4) will have an upper edge from which we can well-determine M_ . Other visible 
particles not involved can be grouped into V with 4-momentum k. In general we refer to any 
process that creates non-zero kx as upstream transverse momentum (UTM) . One type of UTM 
is initial state radiation (ISR), which tends to be small compared to the mass scales involved in 
SUSY processes. Another type of UTM would be decays of heavier particles earlier in the decay 
chain. In the analysis presented in this chapter, we tested the concepts against both k = and 
k < 20 GeV commiserate with what we might expect for ISR. 

We adapt the concept from Mt2 of minimizing the transverse mass over the unknown mo- 
menta to allow for the incorporation of all the available information about the masses. To do 
this we form a new variable M2C which we define as the minimum mass of the second to lightest 



1 We note that the on-shell intermediate case studied by CGHMM is also improved by including the relationship 
measured by the edge in the 11 distribution on each event's analysis. The Y decay channel with an on-shell 
intermediate state X has an edge in the U invariant mass distribution which provides a good determination of 
the relationship maxmf; = (My — M X ){M X — M%)/M x . This relationship forms a surface in Mjv ,Mx ,My space 
that only intersects the allowed points of CGHMM's Fig 3 near the actual masses. We will investigate this case 
in Chapter |SJ 
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new state in the event My constrained to be compatible with the observed 4-momenta of Y"'s 
visible decay products with the observed missing transverse energy, with the four-momenta of 
Y and N being on shell, and with the constraint that M_ = My — Mjy is given by the value 
determined by the end point of the m\i distribution. The minimization is performed over the 
eight relevant unknown parameters which may be taken as the 4-momenta p and q of the particle 
N. We neglect any contributions from unobserved initial state radiation (ISR). Thus we have 

Mf c = min(p + a) 2 (6.1) 

p,q 

subject to the 5 constraints 

(p + a) 2 = (q + P) 2 , (6.2) 

p 2 = q 2 , (6.3) 

Pt = (p + q)T (6-4) 
ViP + a) 2 - VW = M. (6.5) 

where fij< is the missing transverse momentum and (p + q)x are the transverse components of 
p+q. Although we can implement the minimization numerically or by using Lagrange multipliers, 
we find the most computationally efficient approach is to modify the Mt2 analytic solution from 
Lester and Barr |170] . Details regarding implementing M2C and the relation of Mic to Mt2 



and the approach of CCKP are in Sec. 6.2 



Errors in the determined masses propagated from the error in the mass difference in the limit 
of k = are given by 



SM Y = — 1-^ 6M N = -—(l + -^) (6.6) 



where JM_ is the error in the determination of the mass difference M_ . To isolate this source 
of error from those introduced by low statistics, we assume we know the correct M_, and we 



should consider the error described in Eq(6.6) as a separate uncertainty from that reported in 



our initial performance estimates in the Section 6.4 



Because the true p, q are in the domain over which we are minimizing, M2C will always 
satisfy M2C < My- The equality is reached for events with either m\2 or 771,34 smaller than M_, 
with. p z /p Q = a z /a Q , and q z /q = Pz/floi and with the transverse components of a parallel to the 
transverse components of j3. When m\2 = m?>\ = M_ the event only gives information about 
about the mass difference. 

The events that approximately saturate the bound have the added benefit that they are 
approximately reconstructed (p and q are known) . We present a proof of this in Appendix [E] If 
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Y is produced near the end of a longer cascade decay, then this reconstruction would allow us 
to determine the masses of all the parent states in the event. The reconstruction of several such 
events may also aid in spin-correlation studies |141] . 



In order to determine the distribution of Mic for the process shown in Fig 4.4 we computed 



it for a set of events generated using the theoretical cross section and assuming perfect detector 



resolution and no background. Figure 6.1 shows the resulting distribution for three cases: My = 
200 GeV, My = 150 GeV and My = 100 GeV each with M_ = 50 GeV. Each distribution was 
built from 30000 events. Note that the minimum My for an event is M_. The three examples 
each have endpoints that give the mass scale, and we are able to distinguish between different 
My for a given M_. The end-point for My = 100 GeV is clear, and the endpoint for My = 150 
GeV and My = 200 GeV becomes more difficult to observe. The shape of the distribution 



exhibits a surprising symmetry discussed in Sec 6.3 



We can also see that as M + /M_ becomes large, the My determination will be hindered by 
the small statistics available near the endpoint or backgrounds. To alleviate this, we should 
instead fit to the entire distribution. It is clear that events away from the endpoint also contain 
information about the masses. For this reason we propose to fit the entire distribution of M2C 
and compare it to the 'ideal' distribution that corresponds to a given value of the masses. As 
we shall discuss this allows the determination of My with a significant reduction in the number 
of events needed. This is the most important new aspect of the method proposed here. 

6.2 Using Mt2 to Find M<ic and the max Mti Kink 

The calculation of M<ic is greatly facilitated by understanding its relation to Mt2- The variable 
Mt2j which was introduced in by Lester and Summers [T2], is equivalent to 

M| 2 ( X ) =min(p + a) 2 (6.7) 
subject to the 5 constraints 

(p + a) 2 = (q + P) 2 , (6.8) 

p 2 = q 2 (6.9) 

f T = {p + q)T (6.10) 
P 2 = X 2 - (6.H) 

As is suggested in the simplified example of |175j . the minimization over the longitudinal frame 
of reference and center-of-mass energy is equivalent to assuming p and a have equal rapidity 
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Figure 6.1: The distribution of 30000 events in 5 GeV bins with perfect resolution and no 
background. The three curves represent My = 200 GeV (dot-dashed), My = 150 GeV (dotted) 
and My = 100 GeV (solid) each with M_ = 50 GeV. Each distribution cuts off at the correct 

My. 



and q and (3 have equal rapidity. Implementing this Eq(6.7) reduces to the traditional definition 
of the Cambridge transverse mass [^j 



By comparing Mt2(x) as defined above to M2C defined in Eq(6.1), we can see that they 



are very similar with the exception that the constraint Eq(6.5) is replaced by the constraint 



Eq(6.11). M2C can be found by scanning Mt2{x) f° r the- x value that such that the constraint 



in Eq(6.5) is also satisfied. At the value of x that satisfies Mt2(x) ~X = both Eq(6.5) and 



Eq(6.11) are satisfied. 



We can see the M2C and M?2 relationship visually. Each event provides a curve Mt2(x)i 
Fig |6.2| shows curves for four events with Mjv = 50 GeV and My = 100 GeV. For all events 
Mt2(x) is a continuous function of x- CCKP point that out that at x > Mat and at k = 
the maximum Mt2(x) approaches x + M— At x < Mat and at k = the maximum Mt2{x) 
approaches 

maxM T2 (x) = k* + \J (k*) 2 + x 2 (6.12) 

where k* = (M Y -Mfj)/2My is the 2BEMP. This maximum occurs for events with a 2 = (3 2 = 0. 
Putting this together, if Mt2^X = 0) > M__, as is true for two of the four events depicted in 



"We have only tested and verified this equivalence numerically for events satisfying Mt2(x — 0) > M_ 
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Figure 6.2: The Mt2(x) curves for four events with = 50 GeV and My = 100 GeV. 
Only the events whose curves starts off at Mt2(0) > M_ intersect the straight line given by 
Mt2(x) ~ X = M-. The Mt2 at the intersection is M2C for that event. 



Fig. 6.2 then because the event is bounded above by Eq(6.12) and is continuous, then it must 



cross x + M— at x < Mat- At this intersection there is a solution to Mt2(x) = X + where 
Af T 2(x) = min My I constraints = M 2 c- Equivalently 

M 2C = M T2 at x where M T2 (x) = X + M_ if M T2 (x = 0) > M_ (6.13) 
= M_ otherwise. (6-14) 

Assuming A; = 0, the maximum x of such an intersection occurs for x = which is why the 
endpoint of M2C occurs at the correct My and why this corresponds to the kink of CCKP. 
Because Barr and Lester have an analytic solution to Mt2 m Ref. j!70| in the case k = 0, this is 
computationally very efficient as a definition. We will study the intersection of these two curves 
when kx 7^ in Chapter [7j 

6.3 Symmetries and Dependencies of the M2C Distribution 

All transverse mass variables Mt and Mt2 are invariant under longitudinal boosts. Mt2 has 
an additional invariance. CCKP [171j prove that if kx = then M?2 is invariant under back to 
back boosts of the parent particle pair in the plane perpendicular to the beam direction. This 
means that for any event of the topology in Fig. |4.4| with kx = and barring spin correlation 
effects, the M?2 distribution will be the same for a fixed center-of-mass energy as it is for a 
mixed set of collision energies. We now verify this argument numerically. 



In order to numerically determine the distribution of M2C for the processes shown in Fig 4.4 



it is necessary to generate a large sample of "ideal" events corresponding to the physical process 
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shown in the figure. We assume kr = where we expect the back-to-back boost invariance to 
lead to a symmetry of the distribution under changes in the y/s of the collision. We assume each 
branch decays via an off-shell Z°-boson as this is what could be calculated quickly and captures 
the essential elements to provide an initial estimate of our approach's utility. 

Even under these assumptions without knowing about the invariance of Mt2 under back- 
to-back boosts, we might expect that the shape of the distribution depends sensitively on the 
parton distribution and many aspects of the differential cross section and differential decay rates. 
This is not the case; in the case of direct pair production the shape of the distribution depends 
sensitively only on two properties: 

(i) the shape of the mi2 (or equivalently 171,34) distributions. In the examples studied here 
for illustration, we calculate the m\2 distribution assuming it is generated by a particular su- 
persymmetric extension of the Standard Model, but in practice we should use the measured 
distribution which is accessible to accurate determination. The particular shape of m\2 does 
not greatly affect the ability to determine the mass of N and Y so long as we can still find the 
endpoint to determine My — Mn and use the observed raj; distribution to model the shape of 
the M2C distribution. 

(ii) the angular dependence of the iV's momenta in the rest frame of Y . In the preliminary 
analysis presented here we assume that in the rest frame of x£j Xi' s momentum is distributed 
uniformly over the 47T steradian directions. While this assumption is not universally true it 
applies in many cases and hence is a good starting point for analyzing the efficacy of the method. 

Under what conditions is the uniform distribution true? Note that the \2 S s Pi n ^ s the ° nr y 
property of \% that can break the rotational symmetry of the decay products. For \2 S s Pi n to 
affect the angular distribution there must be a correlation of the spin with the momentum which 
requires a parity violating coupling. Consider first the Z° contribution. Since we are integrating 
over the lepton momenta difference^) the parity violating term in the cross section coming 
from the lepton-Z° vertex vanishes and a non-zero correlation requires that the parity violating 
coupling be associated with the neutralino vertex. The Z°-boson neutralino vertex vanishes 
as the Z° interaction is proportional to X^l^l^Xi^li or X^^Xi^u depending on the relative 
sign of and eigenvalues. However if the decay has a significant contribution from an 
intermediate slepton there are parity violating couplings and there will be spin correlations. In 
this case there will be angular correlations but it is straightforward to modify the method to 

take account of correlation^ 

3 M2C only depends on the sum of the two OSSF lepton momenta that follow from a decay of Y. 
4 Studying and exploiting the neutralino spin correlations is discussed further in Refs [1801 118f1 1182] , 
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Figure 6.3: Effect of this maximally spin correlated process on the Mic distribution. Modeled 
masses are My = 200 GeV and Mn = 150 GeV. The solid black distribution is the uncorrelated 
case and red dotted distribution is maximally spin correlated. 



How big an effect could spin correlations have on the shape of the Mic distribution? To 



demonstrate, we modeled a maximally spin-correlated direct production process. Fig 6.3 show 
the spin-correlated process that we consider, and the Mic distributions from this process com- 
pared to the Mic distribution from the same topology and masses but without spin correlations. 
The modeled case has perfect energy resolution, m Vl = GeV, My = 200 GeV, and M^ = 150 
GeV. Our maximally spin-correlated process involves pair production of Y through a pseu- 
doscalar A. The fermion Y in both branches decays to a complex scalar N and visible fermion 
v\ through a purely chiral interaction. The production of the pseudoscalar ensures that the 
Y and Y are in a configuration \f 7 l 1 (\ fl) + | ||)). The Y then decays with N preferentially 
aligned with the spin. The Y decays with N* preferentially aligned against the spin. This causes 
the two sources of missing transverse momentum to be preferentially parallel and pushes Mic 
closer to the endpoint. For this reason the spin-correlated distribution (red dotted distribution) 



is above the uncorrelated distribution (solid thick distribution) in Fig. 6.3 For the remainder 



of the chapter we assume no such spin correlations are present as is true for neutralino pair 
production and decay to leptons through a Z° boson. 

How likely will we have no spin correlations in supersymmetric LHC events? We showed 
earlier that Z°-boson dominated three-body decays lack spin correlations. Even in the case that 
the slepton contribution is significant the correlations may still be largely absent. Because we 
are worried about a distribution, the spin correlation is only of concern to our assumption if a 



mechanism aligns the spin's of the Xi s m the two branches. Table 6.2 shows that most of the 
XiS that we expect follow from decay chains involving a squark, which being a scalar should 
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make uncorrelated the spin of the x% i n the two branches. We would then average over the spin 
states of X2 an d recover the uniform angular distribution of x°'s momentum in x*2 s res t frame. 

Once we have fixed the dependencies (i) and (ii) above, the shape of the distribution is 
essentially independent of the remaining parameters. We calculate the "ideal" distributions for 
M2C assuming that k = and that in the rest frame of Y there is an equal likelihood of N going 
in any of the An steradian directions. The observable invariant a 2 is determined according to 
the differential decay probability of X2 t° e+ e ~ an d Xi through a Z°-boson mediated three- 
body decay. Analytic expressions for cross sections were obtained from the Mathematica output 
options in CompHEP |145] To illustrate the symmetry with respect to changes in y/s and the 



angle of Y's production, we show in Fig 6.4 two cases: 

(1) The case that the collision energy and frame of reference and angle of the produced Y 
with respect to the beam axis are distributed according to the calculated cross section for the 



process considered in Section 6.4 in which x.2 decays via Z° exchange to the three-body state 
l + + l~ + xi-, convoluted with realistic parton distribution functions. 

(2) The case that the angle of the produced Y with respect to the beam axis is arbitrarily 
fixed at 9 = 0.2 radians, the azimuthal angle <j> fixed at radians, and the center-of-mass energy 
set to y/s = 500 GeV. 



The left plot of Fig 6.4 shows the two distributions intentionally shifted by 0.001 to allow us 
to barely distinguish the two curves. On the right side of Fig |6.4| we show the difference of the 
two distributions with the 2 a error bars within which we expect 95% of the bins to overlap 
if the distributions are identical. 

Analytically, this symmetry is because both the Mt2 and M2C are invariant under back to 
back boosts when kx = 0. This symmetry implies that the distributions (1) and (2) above should 



be the same which has been verified numerically in Fig 6.4 In this way we have identified both 
the mathematical origin of the symmetry and numerically verified the symmetry. Thus the M2C 
distribution has a symmetry with respect to changes in y/s when = 0. 

In addition to tests with kx = 0, we also tested that k < 20 GeV does not change the shape 
of the distribution to within our numerical uncertainties. We constructed events with My = 150 
GeV, M N = 100 GeV, \fk? uniformly distributed between 2 and 20 GeV, k/k Q = 0.98, and with 
uniform angular distribution. We found the M2C distribution agreed with the distribution shown 



in Fig. 6.4 within the expected error bars after 10000 events. Scaling this down to the masses 
studied in model PI we trust these results remain unaffected for kx 5j 20 GeV. Introduction of 
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Figure 6.4: Demonstrates the distribution is independent of the COM energy, angle with which 
the pair is produced with respect to the beam axis, and the frame of reference. 

cuts on missing traverse energy and and very large UTM {kx <i M_) changes the shape of the 
distribution. These effects will be studied in Chapter [7j 

Inclusion of backgrounds also changes the shape. Backgrounds that we can anticipate or 
measure, like di-rs or leptons from other neutralino decays observed with different edges can be 
modeled and included in the ideal shapes used to perform the mass parameter estimation. A 
more complete study of the background effects also follows in Chapter [7j 

6.4 Application of the method: SUSY model examples 

To illustrate the power of fitting the full Mic distribution, we now turn to an initial estimate 
of our ability to measure My in a few specific supersymmetry scenarios. Our purpose here is to 
show that fitting the M<ic distribution can determine My and Mn with very few events. We 
include detector resolution effects and assume k = (equivalent to assuming direct production), 
but neglect backgrounds until Chapter [7} We calculate M2C for the case where the analytic Mt2 
solution of Barr and Lester can be used to speed up the calculations as described in Sec |6.2| We 



make the same modeling assumptions described in Sec. 6.3 

Although fitting the M2C distribution could equally well be applied to the gluino mass studied 
in CCKP, we explore its applications to pair-produced Xt We select SUSY models where \% 
decays via a three-body decay to / + + l~~ + Xi- The four momenta a = pi+ +p\- for the leptons 
in the top branch, and the four momenta (3 = pi+ + pi- for the leptons in the bottom branch. 

The production and decay cross section estimates in this section are calculated using Mad- 
Graph/MadEvent |144j and using SUSY mass spectra inputs from SuSpect |183j . We simulate 
the typical LHC detector lepton energy resolution [9j 1184] by scaling the a and j3 four vectors 
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by a scalar normally distributed about 1 with the width of 

5a 0.1 0.003 

— = ^^^= + ———+0.007. 6.15 
a ^a (GeV) a D (GeV) 

The missing transverse momentum is assumed to be whatever is missing to conserve the trans- 
verse momentum after the smearing of the leptons momenta. We do not account for the greater 
uncertainty in missing momentum from hadrons or from muons which do not deposit all their 
energy in the calorimeter and whose energy resolution is therefore correlated to the missing 
momentum. Including such effects is considered in Chapter [7j These finite resolution effects are 
simulated in both the determination of the ideal distribution and in the small sample of events 
that is fit to the ideal distribution to determine My and M^. We do not expect expanded energy 
resolutions to greatly affect the results because the resolution effects are included in both the 
simulated events and in the creation of the ideal curves which are then fit to the low statistics 
events to estimate the mass. 

We consider models where the three-body decay channel for Xq, wm dominate. These models 
must satisfy M-%° — < Mz and must have all slepton masses greater than the M^o. The 
models considered are shown in Table 16.11 The Min-Content model assumes that there are 
no other SUSY particles accessible at the LHC other than x.2 an d X° an d we place and 
Mjfo at the boundary of the PDG Live exclusion limit [41 . SPS 6, PI, and 7 are models 
taken from references [3], [2], and |185j respectively. Each has the x.2 decay channel to leptons 
via a three-body decay kinematically accessible. We will only show simulation results for the 
masses in model PI and SPS 6 because they have the extreme values of M + /M_ with which the 
performance scales. The Min-Content model and the 7 model are included to demonstrate the 
range of the masses and production cross sections that one might expect. 

Bisset, Kersting, Li, Moortgat, Moretti, and Xie (BKLMMX) |186j have studied the four 
lepton with missing transverse momentum Standard-Model background for the LHC. They in- 
cluded contributions from jets misidentified as leptons and estimated about 190 background 
events at a C = 300 fb" 1 which is equivalent to 0.6 fb. Their background study made no refer- 
ence to the invariant mass squared of the four leptons, so we only expects a fraction of these to 
have both lepton pairs to have invariant masses less than M_ . Their analysis shows the largest 
source of backgrounds will most likely be other supersymmetric states decaying to four leptons. 
Again, we expect only a fraction of these to have both lepton pairs with invariant masses within 
the range of interest. The background study of BKLMMX is consistent with a study geared to- 
wards a 500 GeV e + e~ linear collider in Ref. [187 which predicts 0.4 fb for the standard model 
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Model 


Min Content (rcf 41 ) 


SPS 6 (ref [3]) 


PI (ref [2]) 


7 ( Rcf. 185 ) 


1 toll 1 1 1 1 ion 


v? a,nd vS 
are the only 
LHC accessible 
SUSY States 
with smallest 
allowed masses. 


Non Universal 
Gaugino Masses 
m = 150 GeV 
m i/2 = 300 GeV 
tan f3 = 10 

men fill — — 1— 

A = 

Mj = 480 GeV 

M 2 = M 3 = 300 GeV 


mSUGRA 
m = 350 GeV 
m 1/2 = 180 GeV 
tan/3 = 20 
sign(/x) = + 
A = 


Non-Universal 
Higgs Model 
m = 330 GeV 
m l/2 = 240 GeV 

j o on 

tan p = 2\) 

qi crn III) — -I— 

A = 

#2 _ _ (242 GeV) 2 
H 2 = +(373 GeV) 2 




46 GeV 


189 GeV 


69 GeV 


95 GeV 




62.4 GeV 


219 GeV 


133 GeV 


178 GeV 


M+/M- 


6.6 


13.6 


3.2 


3.3 



Table 6.1: Models with X2 decaying via a three-body decay to leptons. We only show simulation 
results for the masses in model PI and SPS 6 because they have the extreme values of M + /M_ 
with which the performance scales. 









Events with 


Model 


a n -0 Direct 




4 leptons +Ef missing 


a x° 2 x° 2 Via 9 or <? 




+ possible extra jets 
C = 300 fb^ 1 


Min Content 


2130 fb 
N/A 


0.067 
0.69 


2893 


SPS 6 


9.3 fb 
626 fb 


0.18 
0.05 


6366 


PI 


35 fb 
12343 fb 


0.025 
0.66 


2310 


7 


17 fb 
4141 fb 


0.043 
0.64 


2347 



Table 6.2: The approximate breakdown of signal events. 

contribution to 4 leptons and missing transverse momentum. The neutralino decay to r leptons 
also provide a background because the r decay to a light leptons I = e, fj, (T T ^,wjT as 0.34) 
cannot be distinguished from prompt leptons. The neutrinos associated with these light leptons 
will be new sources of missing transverse momentum and will therefore be a background to our 
analysis. The di-T events will only form a background when both opposite sign same flavor rs 
decay to the same flavor of light lepton which one expects about 6% of the time. 



Table 6.2 breaks down the LHC production cross section for pair producing two X2 111 eac h 
of these models. In the branching ratio to leptons, we only consider e and /x states as the r will 
decay into a jet and a neutrino introducing more missing transverse momentum. Direct pair 
production of \\ nas a rather modest cross section, but production via a gluino or squark has a 
considerably larger cross section but will be accompanied by additional QCD jets. We do expect 
to be able to distinguish QCD jets from r jets |188] . The events with gluinos and jets will lead 
to considerable k^. In this chapter we assume kx ^$ 20 GeV, but we take up the case of large 
kx in the following chapter. 

We now estimate how well we may be able to measure M^j and in these models under 
these simplifying ideal assumptions. Figures 6.5 and 6.6 show a y 2 fitj^lof the M2C distribution 



3 See Appendix [d| for details of how \ 2 is calculated. 
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Model PI with 250 Events 

Fit Gives M^. = 133 +/- 6 GeV 

True M Xl °=69 GeV, M^,o=133 GeV 




120 125 130 135 140 145 

M Xl ° of Theoretical Distribution (GeV) 



Figure 6.5: x 2 fit °f 250 events from model PI of Ref [2] to the theoretical distributions calculated 
for different M x ° values but fixed M x ° - M x °. The fit gives M x ° = 133 ± 6 GeV. 



Model SPS 6 with 3000 Events 
Fit Gives M^ 2 == 221 +/- 20 GeV 
True M^o = 189 GeV, M X2 ° = 219 GeV 




180 200 220 240 260 280 

M„ 2 « of Theoretical Distribution (GeV) 



Figure 6.6: x 2 fit °f 3000 events from model SPS 6 of Ref [3] to the theoretical distributions 
calculated for different M x » values but fixed M x ° - M x °. The fit gives M x ° = 221 ± 20 GeV. 
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from the observed small set of events to 'ideal' theoretical Mic distributions parameterized by 
M^o. The 'ideal' theoretical distributions are calculated for the observed value of M_ using 
different choices for M%°- A second-order interpolation is then fit to these points to estimate 
the value for M^°- The 1 a uncertainty for is taken to be the points where the x 2 increases 
from its minimum by one. 

The difficulty of the mass determination from the distribution grows with the ratio M+/M-. 



Figures |6.5| and 6.6 show the two extremes among the cases we consider. For the model PI 
M+/M- = 3.2, and for model 7 M + /M_ = 3.3. Therefore these two models can have the 
mass of M^o and M^o determined with approximately equal accuracy with equal number of 



signal events. Figure 6.5 shows that we may be able to achieve ±6 GeV resolution after about 



30 fb _1 . Model SPS 6 shown in Fig 6.6 represents a much harder case because M + /M_ = 13.6. 
In this scenario we can only achieve ±20 GeV resolution with 3000 events corresponding to 
approximately 150 fb -1 . In addition to these uncertainties, we need to also consider the error 



propagated from <5M_ in Eq(6.6). 

Chapter Summary 

We have proposed a method to extract the masses of new pair-produced states based on a 
kinematic variable, M2C*, which incorporates all the known kinematic constraints on the observed 
process and whose endpoint determines the new particle masses. However the method does not 
rely solely on the endpoint but uses the full data set, comparing the observed distribution for 
M2C with the ideal distribution that corresponds to a given mass. As a result the number of 
events needed to determine the masses is very significantly reduced so that the method may be 
employed at the LHC event for processes with electroweak production cross sections. 

This chapter is an initial feasibility study of the method for several supersymmetric models. 
We have made many idealized assumptions which amount to conditions present if the new 
particle states are directly produced. We have included the effect of detector resolution but not 
backgrounds, or cuts. Our modeling assumed that k = 0. We demonstrated that for some of 
the models studied we are able to determine the masses to within 6 GeV from only 250 events. 
This efficiency is encouraging, a study including more of the real-world complications follows in 
Chapter [7j 

The constrained mass variables we advocate here can be readily extended to other processes. 
By incorporating all the known kinematical constraints, the information away from kinematical 
end-points can, with some mild process-dependent information, be used to reduce the number 
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of events needed to get mass measurements. We shall illustrate an extension to three on-shell 
states in Chapter [8} 
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Chapter 7 



The Variable M<ic* Significant 
Upstream Transverse Momentum 

Chapter Overview 

In the previous chapter, we introduced the Mic kinematic variable which gives an event- by-event 
lower bound on the dark-matter particle's absolute mass given the mass difference between the 
dark matter candidate and its parent. The previous chapter focused on direct pair produc- 
tion with minimal initial state radiation (ISR). In this chapter, we introduce a complementary 
variable Mzc,UB which gives an event-by-event upper bound on the same absolute mass. The 
complementary variable is only relevant in the presence of large upstream transverse momen- 
tum (UTM) . Our study shows that the technique presented is as good if not better than other 
model-independent invisible-particle mass determination techniques in precision and accuracy. 

In this chapter, which is based on work first published by Barr, Ross, and the author in 
Ref |16j . we demonstrate the use of the variable M2C and M<ic,\jb in LHC conditions. The 
variables M2c(M_) and M2c,ub(M_) give an event-by-event lower-bound and upper-bound 
respectively on the mass of Y assuming the topology in Fig. |4.4| and the mass difference M_ = 
My — Mtv- To get the mass difference, we use events where Y decays into N and two visible 
states via a three-body decay in which we can easily determine the mass difference from the 
end point of the visible-states invariant- mass distribution, m\ 2 - One might also conceive of 
a situation with M^c supplementing an alternative technique that gives a tight constraint on 
the mass difference but may have multiple solutions or a weaker constraint on the mass scale 
|157| |14j. Given this mass difference and enough statistics, il^c's endpoint gives the mass of Y . 
However the main advantage of the M<ic method is that it does not rely simply on the position 
at the endpoint but it uses the additional information contained in events which lie far from the 
endpoint. As a result it gives a mass determination using significantly fewer events and is less 
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sensitive to energy resolution and other errors. 

To illustrate the method, in this chapter we study in detail the performance of the M%c 
constrained mass variable in a specific supersymmetric model. We study events where each of 
the two branches have decay chains that end with a X2 decaying to a x° an d a P a h °f opposite- 
sign same-flavor (OSSF) leptons. Thus the final states of interest contain four isolated leptons 



(made up of two OSSF pairs) and missing transverse momentum. Fig 4.4 defines the four 
momentum of the particle states with Y = X2.1 N = x°, and the OSSF pairs forming the visible 
particles 1 — 4. Any decay products early in the decay chains of either branch are grouped into k 
which we generically refer to as upstream transverse momentum (UTM). Nonzero k could be the 
result of initial state radiation (ISR) or decays of heavier particles further up the decay chain. 
Events with four leptons and missing transverse momentum have a very small Standard-Model 
background. To give a detailed illustration of the M<ic methods, we have chosen to analyze the 
benchmark point PI from [2 which corresponds to mSUGRA with m = 350 GeV, mi/2 = 180 
GeV, tan/3 = 20, sign(^) = +, A a = 0. Our SUSY particle spectrum was calculated with 
ISAJET [189J version 7.63. We stress that the analysis technique employed applies generically 
to models involving decays to a massive particle state that leaves the detector unnoticed. 

A powerful feature of the M2C distribution is that, with some mild assumptions, the shape 
away from the endpoint can be entirely determined from the unknown mass scale and quantities 
that are measured. The ideal shape fit against early data therefore provides an early mass esti- 
mate for the invisible particle. This study is meant to be a guide on how to overcome difficulties 
in establishing and fitting the shape: difficulties from combinatoric issues, from differing en- 
ergy resolutions for the leptons, hadrons, and missing transverse momentum, from backgrounds, 
and from large upstream transverse momentum (UTM) R As we shall discuss, UTM actually 
provides surprising benefits. 



The chapter is structured as follows: In Section 7.1 we review M2C and introduce the new 
observation that, in addition to an event-by-event lower bound on My, large recoil against 
UTM enables one also to obtain an event- by-event upper bound on My. We call this quantity 



M2C,UB- Section 7.2 describes the modeling and simulation employed. Section 7.3 discusses the 
implications of several effects on the shape of the distribution including the mvi (in our case 
mu) distribution, the UTM distribution, the backgrounds, combinatorics, energy resolution, 



and missing transverse momentum cuts. In Section 7.4 we put these factors together and 



1 Our references to UTM correspond to the Significant Transverse Momentum (SPT), pair production category 
in |176| where SPT indicates that the relevant pair of parent particles can be seen as recoiling against a significant 
transverse momentum. 
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estimate the performance. We conclude in Section 7.4 with a discussion about the performance 
in comparison to previous work. 



7.1 Upper Bounds on My from Recoil against UTM 

We will now review the definition of M%c as providing an event-by-event lower bound on My. 
In generalizing this framework, we find a new result that one can also obtain an upper bound 
on the mass My when the two parent particles Y recoil against some large upstream transverse 
momentum fey. 

7.1.1 Review of the Lower Bound on My 



Fig 4.4 gives the relevant topology and the momentum assignments. The visible particles 1 
and 2 and invisible particle TV" are labeled with momentum ot\ and ai (which we group into 
a = ot\ + a-z) and p, respectively P = fli + 02 and q in the other branch. We assume that 
the parent particle Y is the same in both branches so (p + a) 2 = (q + (3) 2 . Any earlier decay 
products of either branch are grouped into the upstream transverse momentum (UTM) 4-vector 
momentum, k. 

In the previous chapter we showed how to find an event-by-event lower bound on the true 
mass of Mtv and My. We assume that the mass difference M_ = My — Mn can be accurately 
measured from the invariant mass edges max 777,12 or max 711,34. For each event, the variable 
M2C is the minimum value of the mass of Y (the second lightest state) after minimizing over 
the unknown division of the missing transverse momentum between the two dark-matter 
particles N as described in Eq( 6.1]|6.5 ). 



One way of calculating M2C for an event is to use Mt2(xn) [El El US], which provides a 
lower bound on the mass of Y for an assumed mass xn of N. The true mass of Y lies along the 
line Xy{xn) = M_ + xn where we use xy to denote the possible masses of Y and to distinguish 
it from the true mass of Y denoted with My. In other words Mti provides the constraint 
Xy(xn) > Mt2(xn)- Thus we can see that for xn to be compatible with an event, we must 
have M T2 (xn) < Xy(xn) = Xn + M_. 

For a given event, if one assumes a mass xn for N, and if the inequality Mt2(xn) < Xn +M_ 
is satisfied, then there is no contradiction, and the event is compatible with this value of xn- If 
however, Mt2(xn) > M_ +xn, then we have a contradiction, and the event excludes this value 
Xn as a viable mass of N. Using this observation, M2C can be found for each event by seeking 
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the intersection between Mt2{xn) and xn + M_ [15J. Equivalently, the lower bound on My is 
given by My > M_ + Xjv where Xat is the zero of 

ff(Xiv) = M T2 {xn) - xn- M- 

with (/(xSr)<0. (7.1) 

In the case k = 0, the extreme events analyzed in CCKP [174 demonstrate that g{x) will only 
have one positive zero or no positive zeros, and the slope at a zero will always be negative. For no 
positive zeros, the lower bound is the trivial lower bound given by M_. Note that a lower bound 
on the value of My corresponds to a lower bound on the value of Mjy ■ The Appendix in Ref . |15j 



shows that at the zeros of g(xN) which satisfy Eq(7.1), the momenta satisfy Eqns(6.2 6.5). 
7.1.2 A New Upper Bound on My 

If there is large upstream transverse momentum (UTM) (fcy M_) against which the system 
recoils, then we find a new result. Using the Mt2 method to calculate Mic gives one the 
immediate ability to see that My can also have an upper bound when requiring Eqns( (T2]|6.5 ). 



This follows because for large UTM the function g(xN) may have two zeros^] which provides 
both an upper and a lower bound for My from a single event. We have also found regions of 
parameter space where g{x) has a single zero but </(Xjv) > corresponding to an upper bound 
on the true mass of Mn ( and My) and only the trivial lower bound of Mjy > 0. 

We can obtain some insight into the cases in which events with large UTM provides upper 
bounds on the mass by studying a class of extreme event with two hard jets, j a and jp against 
which Y recoils {k = j a + jp). We will describe this extreme event and solve for the regions of 
parameter space for which one can analytically see the intersection points giving a lower bound 
and/or an upper bound. The event is extremal in that Mj-2(xn), which gives a lower bound on 
My, actually gives the true value of My when one selects xn equal to the true mass M/v- 

The ideal event we consider is where a heavier state G is pair-produced on shell at threshold. 
For simplicity we assume the lab-frame is the collision frame. Assume that the Gs, initially at 
rest, decay into visible massless jets j a , jp and the two Y states with the decay product momenta 
a + p and (3 + q. Both jets have their momenta in the same transverse plane along the negative 
x-axis, and both Y's momentum are directed along the x-axis. Finally, in the rest frame of the 
two Ys, both decay such that the decay products visible states have their momentum a and j3 



2 There may be regions in parameter space where function g(x) has more than two zeros, but we have not 
encountered such cases in our simulations. 
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along the x-axis and both invisible massive states N have their two momenta along the negative 
rc-axis. In the lab frame, the four-vectors are given by 



Mg ( ,M? 
2 I M} 



j a = jp = ) {1,-1,0,0} (7.2) 



Maf M 2 
2 I M? 



<* = P = "TT 1 ( 1 - iqT ) {1)1,0,0} (" ;> >> 



P = g = M° I (Mn_My } 00 . (74 ) 

P Q 2 \\M* + Ml) \m 2 Ml) ' ' h 1 



For the event given by Eqns( 7T2]]7.4[ ), we can exactly calculate Mt2(xn) 



M 2 , , 2 , (M 2 N ~ M$) (M 2 N M 2 G - My) 
M T2 (xn) = Xn + 

(M Y - M%) ^4 Ml X %M Y + (M Y - M 2 N M 2 G f 
+ 2My ' (7 - 5j 

This is found by calculating the transverse mass for each branch while assuming xn to be the 
mass of N. The value of p x is chosen so that the transverse masses of the two branches are 
equal. Substituting this value back into the transverse mass of either branch gives Mt2{xn)- 



Fig 7.1 shows g{xN), given in Eq(7.1), for several choices of Mg for the process described 
by Eqs( |7T2f7T4"l ) with M_ = 53 GeV and M N = 67.4 GeV. Because G is the parent of Y, we 
must have M G > My. If My < M G < 2M Y /(M N + My), then M T2 (xn < M N ) is larger 
than xn + M— up until their point of intersection at xn = Mjy- In this case their point of 



intersection provides a lower bound as illustrated by the dotted line in Fig. 7.1 for the case with 



M G = 150 GeV. For 2M Y /(M N + My) < M G < ^M Y /M N there are two solutions 

XJV,Min = M N (7.6) 

(Mjv - My) (-2M Y + M N MlM Y + M 2 N M G ) 

XN ' Max ~ (M N M G + {Mg - 2M Y )M Y )(M N M G + M Y (M G + 2M Y )) ' ^ ' ' 



When M G = y M y /Mn, the function §(xn) has only one zero with the lower bound equalling 
the upper bound at Mat- The solid line in Fig. |7.l| shows this case. Between M Y jMfq < 
M G < V6M Y /^/(M N + M Y )(2M N + My) we again have two solutions but this time with 

(M N - My) (-2M Y + M N M 2 G M Y + M 2 N M G ) 

XN,Min (M N M G + (M G -2My)My)(M N M G + My(M G + 2My)) ^ ' ' 

XiV,Max = M N . (7.9) 



The dashed line in Fig. 7.1 shows this case with M G = 170 GeV. For M G greater than this, we 



have XiV.Max = M N and XiV,Min = 0. 
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Figure 7.1: Shows g(xN) for the extreme event in Eq(7.2 7.4) with M_ = 53 GeV and Mjy 



67.4 GeV. The red dotted line has Mq = 150 GeV and shows an event providing a lower bound 
on My. The blue dashed line Mq = 170 GeV and shows an event with both a lower bound and 
an upper-bound on My. The black solid line shows Mq = <J My- /Mn where the lower bound 
equals the upper bound. 



This example illustrates how Mic can provide both a lower-bound and an upper-bound on 
the true mass for those events with large UTM. The upper-bound distribution provides extra 
information that can also be used to improve early mass determination, and in what follows we 
will refer to the upper bound as Miqjjb- We now move on to discuss modeling and simulation 
of this new observation. 

7.2 Modeling and Simulation 

As a specific example of the application of the M2C method, we have chosen a supersymmetry 
model mSUGRA, m = 350 GeV, m 1/2 = 180 GeV, tan/3 = 20, sign( / u) = +, A D = Q The 
spectrum used in the simulation has M^° = 67.4 GeV and M^° = 120.0 GeV. We have employed 
two simulation packages. One is a Mathematica code that creates the 'ideal' distributions based 
only on very simple assumptions and input data. The second is HERWIG |17| \TE\ IT9] which 
simulates events based on a SUSY spectrum, MSSM cross sections, decay chains, and appropriate 
parton distribution functions. If the simple Mathematica simulator predicts 'ideal' shapes that 
agree with HERWIG generator, then we have reason to believe that all the relevant factors relating 
to the shape are identified in the simple Mathematica simulation. This is an important check 
in validating the benefits of fitting the M<ic and M2c t jjB distribution shape as a method to 
measure the mass of new invisible particles produced at hadron colliders. 



This was model PI from which we also used in [TSj. 
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7.2.1 Generation of "Ideal" Distributions 

Our 'ideal' distributions are produced from a home-grown Monte Carlo event generator written 
in Mathematica. This generater serves to ensure that we understand the origin of the distribution 
shape. It also ensures that we have control over measuring the parameters needed to determine 
the mass without knowing the full model, coupling coefficients, or parton distribution functions. 
We also use this simulation to determine on what properties the ideal distributions depends. 



The simulator is used to create events satisfying the topology shown in Fig 4.4 for a set 
of specified masses. We take as given the previously measured mass difference — M-%° = 
52.6 GeV, which we use in this chapter's simulations. We neglect finite widths of the particle 
states as most are in the sub GeV range for the model we are considering. We neglect spin 
correlations between the two branches. We perform the simulations in the center-of-mass frame 
because M.2C and M2c,ub are transverse observables and are invariant under longitudinal boosts. 
The collision energy y/s is distributed according to normalized distribution 



Is - 4 ML 

p(y/s) = 12 M| * -^—^ (7.10) 

unless otherwise specified. The x*2 ^ s produced with a uniform angular distribution, and all 
subsequent decays have uniform angular distribution in the rest frame of the parent. The UTM 
is simulated by making kx equal to the UTM with k 2 = (100 GeV) 2 (unless otherwise specified), 
and boosting the other four-vectors of the event such that the total transverse momentum is 
zero. As we will show, these simple assumptions capture the important elements of the process. 
Being relatively model independent, they provide a means of determining the mass for various 
production mechanisms. If one were to assume detailed knowledge of the production process, it 
would be possible to obtain a better mass determination by using a more complete simulation like 
HERWIG to provide the 'ideal' distributions against which one compares with the data. Here we 
concentrate on the more model independent simulation to demonstrate that it predicts the M2C 
and M2CUB distributions well-enough to perform the mass determination that we demonstrate 
in this case-study. 

7.2.2 HERWIG "Data" 

In order to obtain a more realistic estimate of the problems associated with collision data, 
we generate samples of unweighted inclusive supersymmetric particle pair production, using 
the HERWIG Monte Carlo program with LHC beam conditions. These samples produce a more 
realistic simulation of the event structure that would be obtained for the supersymmetric model 
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Figure 7.2: The M.2C and Mic,ub distributions of HERWIG events before smearing (to simulate 
detector resolution) is applied. The distributions' end-points show M^g ~ 120 GeV. The top 
thick curve shows the net distribution, the next curve down shows the contribution of only the 
signal events, and the bottom dashed curve shows the contribution of only background events. 

studied here, including the (leading order) cross sections and parton distributions. It includes 
all supersymmetric processes and so contains the relevant background processes as well as the 



particular decay chain that we wish to study. Figure 7.2 shows the M2C and M2CUB distributions 
of a sample of HERWIG generated signal and background events. 

Charged leptons {e^ and fj^) produced in the decay of heavy objects (SUSY particles and 
W and Z bosons) were selected for futher study provided they satisfied basic selection criteria 
on transverse momentum (px > 10 GeV) and pseudorapdity (\rj\ < 2.5). Leptons coming from 
hadron decays are usually contained within hadronic jets and so can be experimentally rejected 
with high efficiency using energy or track isolation criteria. This latter category of leptons was 
therefore not used in this study. The acceptance criterion used for the hadronic final state was 



\t)\ < 5. The detector energy resolution functions used are described in Section 7.3.2 

7.3 Factors for Successful Shape Fitting 

There are several factors that control or affect the shape of the M%c and M2CUB distributions. 
We divide the factors into those that affect the in-principle distribution and the factors that 
affect the observation of the distribution by the detector like energy resolution and selection 
cuts. 

The in-principle distribution of these events is influenced by the presence or absence of spin- 
correlations between the branches, the mu distribution of the visible particles, any significant 
upstream transverse momentum (UTM) against which the system is recoiling {e.g. gluinos or 
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Figure 7.3: We show the M^c and M2CUB ideal distributions for five choices of assuming 
the HERWIG generated mu and UTM distributions. 

squarks decaying further up the decay chain), and background coming from other new- physics 
processes or the Standard Model. As all these processes effectively occur at the interaction 
vertex, there are some combinatoric ambiguities. These are the factors that influence the in- 
principle distribution of events that impinges on the particle detector. 

The actual distribution recorded by the detector will depend on further factors. Some factors 
we are able to regulate - for example cuts on the missing transverse momentum. Other factors 
depend on how well we understand the detector's operation - such as the energy resolution and 
particle identification. 

Where the effect of such factors is significant, for example for the m\2, kx, and background 
distributions, our approach has been to model their effect on the ideal distributions by using 
appropriate information from the 'data', much as one would do in a real LHC experiment. For 
the present our 'data' are provided by HERWIG, rather than LHC events, but the principle is the 
same. 

7.3.1 Factors Affecting the In-principle Distribution 
• Mass Difference and Mass Scale 



The end-point of M2C and M2CUB distributions give the mass of Xt Therefore the mass 
scale, M^°, is a dominant factor in the shape of the 'ideal' distribution. This is the reason we 



can use these distributions to determine the mass scale. Fig 7.3 shows the M2C and M2CUB 
distributions for five choices of assuming the HERWIG generated mu and UTM distributions. 
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How does the shape change with mass scale? The shape is typically sharply peaked at 
M2C = M- followed by a tail that ends at the mass of M%°- The peak at M_ is due to events that 
are compatible with M^.° = 0. We say these events give the trivial constraint. Because we bin the 
data, the height of the first bin depends on the bin size. As M + /M_ = (M^° +M^°)/ (M^° — M^°) 
becomes larger, then the non-trivial events are distributed over a wider range and the endpoint 
becomes less clear. In general if all other things are equal, the larger the mass, the more events 
in the first bin and the longer and flatter the tail. 

The distribution also depends on the mass difference M_ which we assume has been deter- 
mined. We expect that experimentally one should be able to read off the mass difference from 
the ma kinematic end-point with very high precision. Gjelsten, Miller, and Osland estimate this 
edge can be measured to better than 0.08 GeV |157| 1159] using many different channels that 
lead to the same edge, and after modeling energy resolution and background. 

Errors in the mass determination propagated from the error in the mass difference in the 
limit of kx = are given approximately by 



5Myo = 1 + J 5M^o = 1 + + 7.11 



where 5M- is the error in the determination of the mass difference M_. An error in M_ will 
lead to an M<ic distribution with a shape and endpoint above or below the true mass in the 



direction indicated by Eq(7.11). 

To isolate this source of error from the uncertainty in the fit, we assume that the mass 
difference is known exactly in our stated results. In our case an uncertainty of 5M- = 0.08 GeV 
would lead to an additional 5M%-o = ±0.5 GeV to be added in quadrature to the error from 
fitting. 

• Spin Correlations Between Branches 



The potential impact of spin correlation between branches on M2C were studied in Sec 6.3 
There are also no spin correlations if the X2 parents are part of a longer decay chain which 
involves a scalar at some stage as is the case in most of the events in the model we study here. 
In the simple Mathematica simulations, we have assumed no spin dependence in the production 
of the hypothetical ideal distribution. 

• Input m\2 Distributions 
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Figure 7.4: Dependence of M^c distribution on the mu distribution. Left: The mu distri- 
butions. Right: The corresponding M2C distributions. The solid curves show the case where 
the mu distribution when the three-body decay is dominated by the Z boson channel, and 
the dashed curves show the case where the mu distribution is taken directly from the HERWIG 
simulation. 



The mu distribution affects the M^c distribution. Fig |7.4| shows two mu distributions and 
the corresponding M2C distributions with k? = (no UTM). The solid lines show the case where 
the three-body decay from X2 1,0 X\ 1S completely dominated by a Z boson. The dashed line 
shows the case where the mu distribution is extracted from the 'realistic' HERWIG simulation. 
We can see that the M2C distribution is affected most strongly in the first several non-zero bins. 
If we were to determine the mass only from the shape of these first several bins using only the Z 
contribution for the mu difference, we would estimate of the mass to be about 4 GeV below the 
true mass. This can be understood because the shape change of the mu distribution effectively 
took events out of the first bin and spread them over the larger bins simulating the effect of a 
smaller mass. 



• Input Upstream Transverse Momentum Distribution 



As we discussed in Section 7.1 if there is a large upstream transverse momentum (UTM) 
against which the two \2 S recoil) then we have both an upper and lower bound on the mass 



scale. The left frame of Fig. 7.5 shows the UTM distribution observed in the 'realistic' HERWIG 



data. The right frame of Fig. 7.5 shows the M2C and M2CUB distributions for fixed UTM (kx) 
of 0, 75, 175, 275, 375, and 575 GeV all with k 2 = (100 GeV) 2 . As we discuss under the next 
bullet, we also find the distribution is not sensitive to the value of k 2 . For hp > 275 GeV, 
these curves begin to approach a common shape. These are ideal M2C upper and lower bound 
distributions where Mjv = 70 GeV and My = 123 GeV. Notice that there is no upper-bound 
curve for the case with zero k^ UTM. The UTM makes the distribution have a sharper endpoint 
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Figure 7.5: Left: The UTM distribution observed in the HERWIG simulation. Right: Ideal 
M2C upper bound and lower bound distribution for a range of upstream transverse momentum 
(UTM) values (k T = 0,75,175,275,375,575 GeV) where M N = 70 GeV and My = 123 GeV. 

and thereby make the mass easier to determine. This is equivalent to having a sharper kink in 
max Mt2 in the presence of large UTM [176J. 

How do we determine k^ from the data? Because we demand exactly four leptons (two 
OSSF pairs), we assume all other activity, basically the hadronic activity, in the detector is 
UTM. The shape used in the 'ideal' distribution is a superposition of the different fixed UTM 



distributions, shown on right frame of Fig. 7.5 weighted by the observed UTM distribution, 
shown on the left frame of Fig. |7.5| Equivalently, we obtain the ideal distribution by selecting 
kx in the Mathematica Monte Carlo according to the observed UTM distribution. 

• Shape Largely Independent of Parton Distributions and Collision Energy 

In the limit where there is no UTM, then M<ic is invariant under back-to-back boosts of the 
parent particles; therefore, M2C is also invariant to changes in the parton distribution functions. 

How much of this invariance survives in the presence of large UTM? The answer is that it 
remains largely independent of the parton collision energy and largely independent of the mass k 2 



as shown in Fig 7.6 numerically. On the left frame, we show three distributions and in the right 



frame their difference with 2 a error bars calculated from 15000 events. The first distribution 
assumes k T = 175 GeV, k 2 = (100 GeV) 2 , ^fs distributed via Eq(|7.10[). The second distribution 



assumes k T = 175 GeV, k 2 = (2000 GeV) 2 , ^/s distributed via Eq( |7.10[ ). The third distribution 
assumes k? = 175 GeV, k 2 = (100 GeV) 2 , and a fixed collision energy of yfs = 549 GeV. 



Backgrounds 
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Figure 7.6: Figure shows that even with large UTM, the distribution is independent of k 2 
and the parton collision energy to the numerical accuracies as calculated from 15000 events. 
Shown are three distributions and their difference, 
distributed via Eq( |7 10| >. (2) k T = 175 GeV, k 2 = 
(3) k T = 175 GeV, k 2 ^ (100 GeV) 2 , ^ = 549 GeV. 



(1) k T = 175 GeV, k 2 = (100 GeV) 2 , yfs 
(2000 GeV) 2 , <Js distributed via Eq(|7.10|. 
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Figure 7.7: The invariant mass of the OSSF leptons from both branches of forming a Dalitz- 
like wedgebox analysis. The events outside the mu < 53 GeV signal rectangle provide control 
samples from which we estimate the background shape and magnitude. The dark events are 
signal, the lighter events are background. 
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Backgrounds affect the shape and, if not corrected for, could provide a systematic error in the 
estimated mass. In Section 



7.4 



we will see that the position of the minimum x 2 m a fit to M^j 



is barely affected by the background. The main effect of the background is to shift the parabola 
up, giving a worse fit. To improve the fit, we may be able to estimate the M<ic and M2C,ub 
distribution and magnitude of the background from the data itself. We first discuss the sources 
of background, and then we describe a generic technique using a Dalitz-like wedgebox analysis 
to estimate a background model which gives approximately the correct shape and magnitude of 
the background. 

One reason we study the four-lepton with missing transverse momentum channel is because 
of the very-low Standard-Model background \187\ 1186] . A previous study [186] estimates about 
120 Standard-Model four-lepton events (two OSSF pairs) for 100 ftr 1 with a f T > 20 GeV cut. 
They suggest that we can further reduce the Standard-Model background by requiring several 
hadronic jets. Because we expect very-little direct x\ P a i r production, this would have very little 
effect on the number of signal events. Also, because Z°s are a part of the intermediate states 
of these background processes, very few of these events will have mu significantly different from 
M z . 

What is the source of these Standard-Model backgrounds? About 60% is from Z-pair pro- 
duction events with no invisible decay products, in which the missing transverse momentum can 
only arise from experimental particle identification and resolution errors. This implies that a 
slightly stronger cut could further eliminate this background. Another 40% are due to t, i, Z 
production. Not explicitly discussed in their study but representing another possible source of 
backgrounds are events containing heavy baryons which decay leptonically. If we assume 6-quark 
hadrons decay to isolated leptons with a branching ratio of 0.01, then LHC t, i production will 
lead to about 10 events passing these cuts for 100 fb _1 where both OSSF leptons pairs have 
mu < M Z - 

Tau decays also provide a background for our specific process of interest. The process X2 ~* 
t + t~x° will be misidentified as e~e + or about 3% of the time (6% total). Because the 

r decays introduce new sources of missing transverse momentum (v T ), these events will distort 
the M2C calculation. This suggests that the dominant background to the X21 X2 ~ y ^ + fr+ 
hadrons will be from other SUSY processes. 

We now create a crude background model from which we estimate the magnitude and dis- 
tribution of the background using the 'true' HERWIG data as a guide. We follow the suggestion 
of Ref. [1601 186 and use a wedgebox analysis plotting the invariant mass m ee against m^u 
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to supplement our knowledge of the background events mixed in with our signal events. This 



wedgebox analysis, seen in Fig. 7.7 for our HERWIG simulation, shows patterns that tell about 
other SUSY states present. The presence of the strips along 91 GeV indicate that particle states 
are being created that decay to two leptons via an on-shell Z. The observation that the inten- 
sity changes above and below m™ = 53 GeV shows that many of the states produced have one 
branch that decays via a x\ ana - the other branch decays via an on-shell Z. The lack of events 
immediately above and to the right of the (53 GeV, 53 GeV) coordinate but below and to the 
left of (91 GeV, 91 GeV) coordinate suggest that symmetric process are not responsible for this 
background. 

We also see the density of events in the block above 53 GeV and to the right of 53 GeV suggest 
a cascade decay with an endpoint near enough to 91 GeV that it is not distinguishable from 
Mz- Following this line of thinking, we model the background with a guess of an asymmetric set 
of events where one branch has new states G, X and iV with masses such that the mu endpoint 
is 



max mf, (odd branch) = {Mg Mn> = (85 GeV) 2 (7.12) 



and the other branch is our X2 decay. The masses one chooses to satisfy this edge did not prove 
important so long as the mass differences were reasonably sized; we tried several different mass 
triplets ending with the LSP, and all gave similar answers. 

We now describe the background model used in our fits. One branch starts with a massive 
state with Mq = 160 GeV which decays to a lepton and a new state Mx = 120 GeV which 
in-turn decays to a lepton and the LSP. The second branch has our signal decay with the \2 
decaying to x° an d two leptons via a three-body decay. We added UTM consistent with that 
observed in the events. 

By matching the number of events seen outside the mu < 53 GeV region, we estimate the 
number of the events within the signal cuts that are due to backgrounds. We estimate 0.33 of 
the events with both OSSF pairs satisfying mu < 53 GeV are background events. The model 
also gives a reasonable distribution for these events. Inspecting the actual HERWIG results showed 
that actual fraction of background events was 0.4. If we let the fraction be free and minimize 
the x 2 with respect to the background fraction, we found a minimum at 0.3. 

Our background model is simplistic and does not represent the actual processes, but it does 
a good job of accounting for the magnitude and the shape of the background mixed into our 
signal distribution. Most of the HERWIG background events came from W and charginos which 
introduce extra sources of missing transverse momentum. Never the less, the shape fits very 
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accurately and the performance is discussed in Section 7.4 It is encouraging that our estimate 
of the background shape and magnitude is relatively insensitive to details of the full spectrum. 
Even ignoring the background, as we will see in Section 7.4 still leads to a minimum x 2 a t the 
correct mass. 

• Combinatoric Ambiguities 

If we assume that the full cascade effectively occurs at the primary vertex (no displaced 
vertices), then the combinatoric question is a property of the ideal distribution produced in the 
collisions. There are no combinatoric issues if the two opposite-sign same- flavor lepton pairs are 
each different flavors (e.g. the four leptons are e + , e~ , fi + and However if all four leptons 
are the same flavor, we have found that we can still identify unique branch assignments 90% 
of the time. The unique identification comes from the observation that both pairs must have 
an invariant mass raj; less than the value of the max mu edge. In 90% of the events, there is 
only one combination that satisfies this requirement. This allows one to use 95% of the four 
lepton events without ambiguity. The first 50% are identified from the two OSSF pairs being 
of different flavors and 90% of the remaining can be identified by requiring that both branches 
of OSSF lepton for an event's combinatoric pairing satisfy m^Event < max ah Events mil- The 
events which remain ambiguous have two possible assignments, both of which are included with 
a weight of 0.5 in the distribution. 

7.3.2 Factors Affecting Distribution Recorded by the Detector 

As just described, the 'ideal' in-principle distribution is created from the observed mu distri- 
bution and the observed UTM distribution. We include combinatoric effects from events with 
four leptons of like flavors. Last, we can estimate the magnitude of background events and their 
M2C and M2CUB shape. We now modify the in-principle distribution to simulate the effects of 
the particle detector to form our final 'ideal' distribution that includes all anticipated effects. 
The two main effects on the M<ic and M2CUB distributions are the energy resolution and the 
Pt cuts. 

• Shape Dependence on Energy Resolution 

Energy resolution causes the M2C and M2c,ub distributions to be smeared. Here we as- 
sume the angular resolution is negligible. For both the Mathematica Monte Carlo model and 



97 



the HERWIG events we simulate the detector's energy resolution by scaling the four vectors for 
electrons, muons, and hadrons by 

SE e 0.1 0.003 „ ,„ , 

+ -^+0.007 (7.13) 



E e \[E e E e 

■ 0.03 (7.14) 



5 Eh 0.58 0.018 

+ -^+0.025 (7.15) 



Eh \[E~h Eh 

respectively [190 |9 . The muon energy resolution is different because they are typically not 
contained by the calorimeter. A more detailed detector simulation is of course possible, but 
since we do not know the true behavior of any LHC detector until the device begins taking data, 
a more sophisticated treatment would be of limited value here. In practice the dependence of 
the ideal distribution shapes on the missing transverse momentum resolution should reflect the 
actual estimated uncertainty of the missing transverse momentum of the observed events. 

Smearing of the distributions decreases the area difference between two normalized distri- 
butions, thereby decreasing the precision with which one can determine the mass from a given 



number of signal events. This expanded uncertainty can be seen in Section 7.4 

The M2C calculations depend on the mass difference, the four-momenta of the four leptons, 
and the missing transverse momentum. As the lepton energy resolution is very tight, the missing 
transverse momentum's energy resolution is dominated by the hadronic energy resolution. We 
model the energy resolution of the UTM as a hadronic jet. This significantly increases the 
uncertainty in the missing transverse momentum because hadrons have about five times the 
energy resolution error. 

In our Mathematica model, we represent the UTM as a single four-vector k, but in reality 
it will be the sum of many four-vectors. Because we apply the energy resolution smearing to k, 
if k is small the simple Mathematica model will have a smaller missing transverse momentum 
resolution error. However, an events with almost UTM could have a large missing momentum 
energy resolution if it has a lot of hadronic jets whose transverse momentum mostly cancels. Fig 



7.5 shows that most of the time we have considerable hadronic UTM, so this effect is a minor 



correction on our results. 

• Shape Dependence on Missing Transverse Momentum Cuts 

A key distinguishing feature of these events is missing transverse momentum. To eliminate 
the large number of Standard Model events with four-lepton and with no fx-, we wm need to cut 
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Figure 7.8: The missing transverse momentum vs M 2 q values for HERWIG data. This shows that 
a Pt > 20 GeV cut would not affect the distribution for M 2 c > 65 GeV. 



on this parameter. Fig. 7.8 shows the HERWIG simulation's missing transverse momentum versus 
the M 2 c- A non-trivial M 2 q requires substantial fx- Small °f l ess than about 20 GeV only 
affects the M 2 c shape below about 65 GeV. The shape of the M 2 c < 65 GeV therefore will 
require a higher fidelity model from which to train the shapes. Instead, we just choose to not 
fit bins with M 2C < 65 GeV. 

All events near the end of M 2 c,ub distribution require significant Pj<, therefore cuts will 
not affect the part of this distribution which we fit. The number of events with no non-trivial 
upper-bounds will also be affected by fx cuts. We only fit the M 2 c,ub distribution up to about 
233 GeV. 

7.4 Estimated Performance 

Determining the mass based on the shape of the distribution enables one to use all the events 
and not just those near the end point. We fit both upper-bound and lower-bound shapes to the 
data as described in the Appendix[D] As one expects, fitting the lower-bound shape more tightly 
constrains the mass from below and fitting the upper-bound shape more tightly constrains the 
mass from above. Combining the two gives approximately even uncertainty. We calculate ideal 
distributions assuming at five values 50, 60, 70, 80, and 90 GeV. We then fit a quadratic 
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Figure 7.9: The result of x 2 fits to the data with differing assumptions for 100 fb -1 (left 
panel) and 400 fb" 1 (right panel). The thick line with filled squares shows the final result with 
all cuts, resolution error, combinatorics, and backgrounds included and estimated in the shape 
fitting. This gives us = 63.2 ±4.1 GeV with 700 events (signal or background) representing 
100 fb -1 . After 400 fb -1 this improves to = 66.0 ± 1.8 GeV. The error-free best case gives 
Mx° = 67.0 ± 0.9 GeV. The correct value is A% = 67.4 GeV. 

interpolation through the points. Our uncertainties are based on the value where x 2 increases 
by 1 from its minimum of this interpolation. This uncertainty estimate agrees with about 2 /3 
of the results falling within that range after repeated runs. Our uncertainty estimates do not 



include the error propagated from the uncertainty in the mass difference (see Eq(7.11)). 

We present results for an early LHC run, about 100 fb , and for the longest likely LHC run 
before an upg rade, about 400 fb -1 . After about 100 fb" 1 , we have 700 events ( about 400 signal 
and 300 background). After 400 fb" 1 , we have about 2700 events (about 1600 signal and 1100 
background). Only 4 events out of 1600 are from direct pair production. Most of our signal 
events follow at the end of different decay chains starting from gluinos or squarks. The upstream 
decay products produce significant UTM against which the two x% parent particles recoil. 

First for the ideal case. After 400 fb" 1 , using only signal events and no energy resolution, the 



X fits to the predicted shapes give = 67.0 ± 0.9 GeV (filled circles in Fig. 7.9). This mass 



determination can practically be read off from the endpoints seen in Fig. 7.2 the M^c endpoint 
is near 120 GeV and subtracting the mass differences gives = 120 GeV — M_ = 67 GeV. 
We now explore how well we can do with fewer events and after incorporating the effects listed 
in Section [?T3l 

How does background affect the fit? If we ignore the existence of background in our sample, 
and we fit all the events to the signal-only shapes, then we find a poor fit shown as the empty 



circle curve in Fig. 7.9. By poor fit, we mean the x 2 1S substantially larger than the 72 bins 
being compared (36 bins from each the upper-bound and lower-bound distributions). Despite 
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this worse fit, the shape fits still give a very accurate mass estimate: = 65.4 db 1.8 GeV 

after 100 ftT 1 and Mjj = 67.4 ± 0.9 GeV after 400 fb" 1 . At this stage, we still assume perfect 
energy resolution and no missing transverse momentum cut. 



Next, if we create a background model as described in Section 7.3, we are able to improve 
the x 2 fit to nearly 1 per bin; the mass estimate remains about the same, but the uncertainty 
increases by about 20%. We find a small systematic shift (smaller than the uncertainty) in our 
mass prediction as we increase the fraction of the shape due to the background model vs the 
signal model. As we increased our fraction of background, we found the mass estimate was 
shifted down from 66.5 at 0% background to 65.6 when we were at 60% background. The best 
X 2 fit occurs with 30% background; which is very close to the 33% we use from the estimate, but 
farther from the true background fraction of about 40%. With 400 fb -1 of data, the systematic 
errors are all but eliminated with the endpoint dominating the mass estimate. These fits are 
shown as the triangles with dashed lines and give = 65.1 ±2.4 GeV which after the full run 
becomes A% = 67.3 ±1.1 GeV. 



Including energy resolution as described in Section 7.3 shows a large increase in the un- 
certainty. The dashed-line with empty square markers shows the x 2 fit when we include both 
a background model and the effect of including energy resolution. These fits are shown as the 
empty squares with dashed lines and give = 63.0 ±3.6 GeV which after the full run becomes 
Mg» = 66.5 ± 1.6 GeV. 

The final shape factor that we account for are the cuts associated with the missing transverse 
momentum. After we apply cuts requiring f T > 20 GeV and fit only M^c > 65 GeV we have 
our final result shown by the thick lines with filled squares. This includes all cuts, resolution 
error, combinatorics and backgrounds. We find M^j = 63.2 ±4.1 GeV with 700 events (signal or 
background) representing 100 fb -1 , and after 400 fb -1 this improves to = 66.0 ± 1.8 GeV. 
The true mass on which the HERWIG simulation is based is = 67.4, so all the estimates are 
within about 1 a of the true mass. 



Fig 7.10 shows the ideal curve expected if = 70 GeV including all effects from energy 
resolution, background, combinatoric, and fj> cuts. The x 2 corresponds to the solid square on 



the left panel of Fig. 7.9 



The error in mass determination obtained with limited statistics can be estimated using 
Poisson statistics. In our studies we find that, as one would expect, increasing the number of 
events by a factor of four, we bring down our error by about a factor of two. This means that one 
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Figure 7.10: HERWIG data for 100 fb _1 (thick line) and the smooth ideal expectation assuming 
Mjjo = 70 GeV generated by Mathematica with all resolution, background, and combinatoric 
effects included (thin line). The \ 2 °f this curve to the HERWIG gives the solid-black square on 
the left frame of Fig. |7.9[ 

could expect ±8 GeV after about 25 fb _1 which represents 100 signal events and 75 background 
events. 

Chapter Summary 

Despite adding some of the complicating effects one would encounter with real data, we have 
discovered other factors which demonstrate that one could obtain an even better precision than 
estimated in Chapter [6j There, we used only our simple Mathematica model and neglected most 
sources of realistic uncertainty. We assumed all the events could be modeled as being direct 
production without spin-correlations. With these simplifications, we argued the mass could be 
determined to ±6 GeV using 250 signal events. 

In this chapter, we performed a case study to show that that the relevant Mic and M2C,UB 
shapes can be successfully determined from the mass difference M_ , the ma distribution obser- 
vation, the upstream transverse momentum (UTM) distribution observation. We included and 
accounted for many realistic effects: we modeled the large energy-resolution error of hadronic 
jets. We also included the effects of backgrounds, fx cuts, and combinatorics. Our signal and 
backgrounds were generated with HERWIG. We discussed how a Dalitz-like plot can estimate the 
background fraction and shape. Observed inputs were used in a simple model to determine ideal 
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distribution shapes that makes no reference to parton distribution functions, cross sections, or 
other model-dependent factors that we are not likely to know early in the process. 

Despite these extra sources of uncertainty, we found a final mass determination of ±4.1 GeV 
with about 400 signal events which is still better than the appropriately scaled result from RS 
[15] . The sources of the mass determination improvement are twofold: (1) the prediction and 
fitting of upper-bound distribution, and (2) the sharper end-point in the presence of large UTM. 
Under equivalent circumstances, the sharper endpoint is enough to give a factor of 2 improvement 
in the uncertainty over the direct production case assumed in [15]. Fitting the upper bound 
tends to improve the determination by an additional factor of y/2. This improvement is then 
used to fight the large hadronic-jet energy resolution and background uncertainty. 

Mass determination using M2C and M2CUB applies to many other processes. We have 
focused on cases where the mass difference is given by the end-point of an mu distribution 
involving a three-body decay. If there is not a three-body decay, then the mass difference may 
be found by applying other mass determination techniques like the mass shell techniques (MST) 
[166| 1167} 1168] or edges in cascade decays |156| I158[ 1159] or Mt2 at different stages in symmetric 
decay chains [14J. 

How does our method's performance compare to previous mass determination methods? 
Firstly, this technique is more robust than the max Mt2 'kink' because in fitting to the shape of 
the distribution, it does not rely entirely on identification of the events near kinematic boundary. 
One can view M2C and M2CUB as variables that event-by-event quantify the 'kink'. Other than 
the 'kink' technique, the previous techniques surveyed in Chapter [4] apply to cases where there 
is no three-body decay from which to measure the mass difference directly. However, each of 
those techniques still constrains the mass difference with great accuracy. The technique of [156 , 
157, 158, 159 which uses edges from cascade decays determines the LSP mass to ±3.4 GeV with 
about 500 thousand events from 300 fb _1 . The approach of [167] assumes a pair of symmetric 
decay chains and assumes two events have the same structure. They reach ±2.8 GeV using 700 
signal events after 300 fb _1 , but have a 2.5 GeV systematic bias that needs modeling to remove. 
By comparison, adjusting to 700 signal events we achieve ±2.9 GeV without a systematic bias 
after propagating an error of 0.08 GeV in the mass difference and with all discussed effects. 
Uncertainty calculations differ amongst groups, some use repeated trial with new sets of Monte 
Carlo data, and others use x 2 - Without a direct comparison under like circumstance, the optimal 
method is not clear; but it is clear that fitting the M<ic and M2CUB distributions can determine 
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the mass of invisible particles at least as well, if not better than the other known methods in 
both accuracy and precision. 

In summary, we have developed a mass determination technique, based on the constrained 
mass variables, which is able to determine the mass of a dark-matter particle state produced at 
the LHC in events with large missing transverse momentum. The M2C method, which bounds 
the mass from below, was supplemented by a new distribution Micjjb which bounds the mass 
from above in events with large upstream transverse momentum. A particular advantage of 
the method is that it also obtains substantial information from events away from the end point 
allowing for a significant reduction in the error. The shape of the distribution away from the 
end-point can be determined without detailed knowledge of the underlying model, and as such, 
can provide an early estimate of the mass. Once the underlying process and model generating 
the event has been identified the structure away from the end-point can be improved using, for 
example, HERWIG to produce the process dependent shape. We performed a case-study simulation 
under LHC conditions to demonstrate that mass-determination by fitting the Mic and M2CUB 
distributions survives anticipated complications. With this fitting procedure it is possible to 
get an early measurement of the mass - with just 400 signal events in our case study we found 
we would determine = 63.2 ± 4.1. The ultimate accuracy obtainable by this method is 
M^o = 66.0 ±1.8 GeV. We conclude that this technique's precision is as good as, if not better 
than, the best existing techniques. 
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Chapter 8 

The Variable M 3C : On-shell 
Interlineate States 



Chapter Overview 

The main concept of the constrained mass variable M2C PS] p2| is that after studying several 
kinematic quantities we may have well determined the mass difference between two particle 
states but not the mass itself. We then incorporate these additional constraints in the analysis 
of the events. We check each event to test the lower bounds and upper bounds on the mass 
scale that still satisfies the mass difference and the on-shell conditions for the assumed topology. 
Because the domain over which we are minimizing contains the true value for the mass, the 
end-points of the lower bounds and upper bounds distributions' give the true mass. 

The subject of this chapter is extending the constrained mass variable to the case with three 



new on-shell states as depicted in Fig 4.3 The constrained mass variable for this case will be 
called M30 We structure the chapter around a case study of the benchmark point SPS la [3]. 
In this study, the three new states are identified as Y = yQ, X = I and N = x°. The visible 
particles leaving each branch are all opposite-sign same-flavor (OSSF) leptons (/i or e). This 
allows us to identify any hadronic activity as upstream transverse momentum. 



The chapter is structured as follows: Sec. 8.1 introduces the definition of M%c and how to 



calculate it. Sec. 8.3 discuses the dependence of M%c on complications from combinatorics, large 
upstream transverse momentum (UTM), fix cuts, parton distributions, and energy resolution. 
Sec |8.4| applies M^c variables to HERWIG data from the benchmark supersymmetry spectrum 
SPS la. Finally we summarize the chapter's contributions. 
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8.1 Introducing M%c 

We will now introduce the definition of M 3 c, how to calculate it, its relationship to previous 
mass shell techniques. 

8.1.1 Definition of M 3C 

The upper bound and lower bound on the mass of the third lightest new particle state in the 
symmetric decay chain is the constrained mass variable M 3 c,lb and M^cub- This variable 



applies to the symmetric, on-shell intermediate state, topology from Fig. 4.3 which depicts two 
partons that collide and produce some observed upstream transverse momentum (UTM) with 
four momenta k and an on-shell, pair-produced new state Y. On each branch, Y decays to 
on-shell intermediate particle state X and a visible particle v\ with masses Mx and m Vl . Then 
X decays to the dark-matter particle N and visible particle V2 with masses Mn and m V2 . The 
four-momenta of vi, V2 and ./V are respectively a\, a2 and p on one branch and j3i, 02 and q in 
the other branch. The missing transverse momenta is given by the transverse part of p + q. 

We initially assume that we have measured the mass differences from other techniques. For 
an on-shell intermediate state, there is no single end-point that gives the mass difference. The 



short decay chain gives a kinematic endpoint maxmi2 described in Eq(4.3) that constrains a 



combination of the squared mass differences. Unless two of the states are nearly degenerate, the 



line with constant mass differences lies very close to the surface given by Eq(4.3). The two mass 
differences are often tightly constrained in other methods. The mass differences are constrained 
to within 0.3 GeV from studying long cascade decay chains where one combines constraints from 
several endpoints of different invariant mass combinations |157] . The concepts from Chapter [5] 
also provide another technique to determine the mass differences. After initially assuming that 
we know the mass difference, we show that our technique can also find the mass differences. 
The M3C distribution shape is a function of both the mass scale and mass differences. We can 
constrain both the mass differences and the mass scale by fitting the max m\2 edge constrains and 
the ideal M 3C {M N , AM Y n, AM X n) distribution shapes to the observed M 3C {AM Y n, AM X n)- 
To find all three parameters from this fit, we will take Mat, AM Y n, and AMxn as independent 
variables. 

For this first phase of the analysis, let's assume the mass differences are given. For each 
event, the variable M 3 c,lb is the minimum value of the mass of Y (third lightest state) after 
minimizing over the unknown division of the missing transverse energy between the two dark 



106 



End Point Gives Mass Scale 



o 



T3 

U 
N 

s 

o 



0.04 



0.03 



0.02 



0.01 - 



0.00 




100 



150 200 250 

M 3c ,lb ^3 cub (GeV) 



300 



Figure 8.1: Ideal M^clb and M^c,ub distribution for 25000 events in two cases both sharing 
AM YN = 100 GeV and AM X n = 50 GeV. The solid, thick line shows My = 200 GeV, and the 
dashed, thin line shows My = 250 GeV. 
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where AMyjv = My — Mx and AMxn = Mx — Mx- There are eight unknowns in the four 
momenta of p and q and seven equations of constraint. Likewise we define M^c,ub as the 
maximum value of My compatible with the same constraints. We discuss how to numerically 



implement this minimization and maximization in Sec. 8.2 Because the true p and q are 



within the domain over which we are minimizing ( or maximizing), the minimum (maximum) is 
guaranteed to be less than (greater than) or equal to My. 
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Figure 8.1 shows an ideal M^c,LB and M^c,ub distributions for 25000 events in two cases both 
sharing AMyn = 100 GeV and AMxn = 50 GeV. The dashed line represents the distributions 
from events with My = 250 GeV, and in the solid line represents the distributions from events 
with My = 200 GeV. One can clearly see sharp end-points in both the upper bound and lower 
bound distributions that give the value My. The upper-bound distribution is shown in red. 

We expect an event to better constrain the mass scale if we are given additional information 
about that event. In comparison to M 2 c where Y decays directly to N, the on-shell intermediate 
state has an additional state X. The extra state X and information about its mass difference 
AMxn enables M^c to make an event-by-event bound on My stronger than in the case of 
M 2 q- We will see that this stronger bound is partially offset by greater sensitivity to errors in 
momentum measurements. 

The variable M^c, hke other variables we have discussed M 2 Ci Mt 2 Mt and Mqt, is in- 
variant under longitudinal boosts of its input parameters. We can understand this because all 
the constraint equations are invariant under longitudinal boosts. The unknown p and q are 
minimized over all possible values fitting the constraints so changing frame of reference will not 
change the extrema of the Lorentz invariant quantity (p + a) 2 . 

8.2 How to calculate M%c 

To find the M3C*, we observe that if we assume masses of Y, X, and N to be {xy,Xx,Xn) Q 
with the given mass differences then there are eight constraints 



(8.9) 

( X y - AM YN + AM X n? (8.10) 
(xy-AMyn) 2 (8.11) 

and eight unknowns, p^ and q^. The spatial momenta p and q can be found as linear functions 
of the th component of p and q by solving the matrix equation 
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1 We again use x to distinguish hypothetical masses (xy, Xx, Xn) from the true masses (My, Mx, Mn) 
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where a = ot\ + «2 and (3 = fii + 02- We substitute p and q into the on-shell constraints 



giving two quadratic equations for p Q and q Q . These give four complex solutions for the pair p a 
and q a . We test each event for compatibility with a hypothetical triplet of masses (xy,Xx, Xn) = 
(xy, Xy — AMy^r + AMxn, Xy — AMyj). If there are any purely real physical solutions where 
(p > and q D > 0), then we consider the mass triplet (xy, Xx, Xn) viable. 

As we scan xy, a solution begins to exist at a value less than or equal to My and then 
sometimes ceases to be a solution above My. Sometimes there are multiple islands of solutions. 
To find the M%c, we can test each bin starting at xy = AMyjy along the path parameterized 
by xy an d the mass differences to find the first bin where at least one physical solution exists. 
This is the lower bound value of M%c for the event. 

Likewise for an upper bound. We begin testing at the largest conceivable mass scale we 
expect for the Y particle state. If a solution exists, we declare this a trivial M%c,UB- If no 
solution exists, then we search downward in mass scale until a solution exists. 

A faster algorithm involves a bisection search for a solution within the window that starts at 
AMyn and ends at our highest conceivable mass. We then use a binary search algorithm to find 
at what xy the solution first appears for M^q,lb or at what xy the solution disappears giving 
M%c,UB- There are rare events where there are multiple islands of solutions. This occurs in 
about 0.01% of the events with UTM and in about 0.1% for hr = 250 GeV. In our algorithm 
we neglect windows of solutions more narrow than 15 GeV. We report the smallest edge of 
the lower-mass island as the lower bound and the upper edge of the larger-mass island as the 
upper bound. Because of the presence of islands, we are not guaranteed that solutions exist 
everywhere between M^c,lb and M^c,ub- With the inclusion of energy resolution errors and 
background events, we also find cases where there are no solutions anywhere along the path 
being parameterized. If there is no solutions anywhere in the domain we make M%c,lb to be 
the largest conceivable mass scale, and we set M^c,ub = AMy/v- 

8.2.1 Comparison to other Mass Shell Techniques 

The variable M%c is a hybrid mass shell technique |168j . In Chapter [4] we reviewed other mass 
shell technique that measure the mass in the case of three new states. Cheng, Gunion, Han 
Marandella, McElrath (CGHMM) |166] describe counting solutions at assumed values for the 



P 2 o ~ (p(Po,qo)) 2 = Xn 
<7o - (q(Po,q )) 2 = Xn 
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mass for Y, X, and N. By incorporating a minimization or maximization, we enhance CGHMM's 
approach because we have a variable whose value changes slightly with slight changes of the 
inputs instead of the binary on-off that CGHMM has with the existence of a solution^] We 
also incorporate knowledge of the added information from other measurements which accurately 
determine the mass differences. Finally, the quantity M$c can form a distribution whose shape 
tells us information about the masses. Because for most events there is only one 'turn-on' 
point below My, the distribution M^c,LB is very similar to the derivative of the figure 8 of 
CGHMM[166j to the left of their peak and M^c,ub is similar to the negative of the derivative 
to the right of their peak. They differ in that there may be multiple windows of solutions; 
also CGHMM's Fig 8 is not exactly along the line of fixed mass differences; and the effect of 
backgrounds and energy resolution are dealt with differently. 

We also hope to show that the use of the distribution's shape enables us to exploit the 
essentially non-existent dependence of the distributions on the unknown collision energies and 
incorporate the dependency on UTM directly. This diminishes the dependence of the measure- 
ment on the unknown model while still allowing us to exploit the majority of the distribution 
shape in the mass determination. 

After studying previous MSTs, we were tempted to use Bayes theorem with a parton distri- 
bution function as a likelihood function as was done in Goldstein and Dalitz j!61] and Kondo, 
Chikamatsu, Kim |162j (GDKCK). They used the parton distribution function to weight the 
different mass estimates of the top-quark mass (My in our topology). We found that such a 
weighting leads to a prediction for My much smaller than the true value. This can be understood 
because the parton distributions make collisions with smaller center-of-mass energies (small x) 
more likely, therefore the posterior will prefer smaller values of My which are only possible for 
smaller values of x. Only if one includes the cross-section for production, i.e. the likelihood of 
the event existing at all, in the Bayes likelihood function will we have the appropriate factor 
that suppresses small values of x and therefore small values of My. This balance therefore leads 
to the maximum likelihood (in the limit of infinite data) occurring at the correct My. Unfor- 
tunately, inclusion of the magnitude of the cross section introduces large model dependence. In 
the case of the top-quarks mass determination, the GDKCK technique gives reasonable results. 
This is because they were not scanning the mass scale, but rather scanning xy (the top quark 
mass) while assuming xn = M^r = and xx = Mx = M\y. The likelihood of solutions as 
one scans xy rapidly goes to zero below the true top-quark mass M top - The parton distribution 
2 I am grateful to Chris Lester for pointing out to me the importance of this feature. 
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suppresses the likelihood above the true Mt op - The net result gives the maximum likelihood near 
the true top-quark mass but suffers from a systematic bias [163 [164 that must be removed by 
modeling [165 . 

8.3 Factors for Successful Shape Fitting 

One major advantage of using the M^c distribution (just as the M2C distribution) is that the 
bulk of the relevant events are used to determine the mass and not just those near the endpoint. 
To make the approach mostly model independent, we study on what factors the distributions 
shape depends. We show that there is a strong dependence on upstream transverse momentum 
(UTM) which can be measured and used as an input and therefore does not increase the model 
dependency. We show there is no numerically significant dependence on the collision energy 
which is distributed according to the parton distribution functions. This makes the distribution 
shape independent of the production cross section and the details of what happens upstream 
from the part of the decay chain that we are studying. We model these effects with a simple 
Mathematica Monte Carlo event generator assuming My = 200 GeV, Mx = 150 GeV, and 
M N = 100 GeV. 

• Effect of Combinatorics Ambiguities 



Just as in the topology in Fig 4.4 studied earlier, where X2 decays via a three body decay, 
the branch assignments can be determined by either distinct OSSF pairs or by studying which 
OSSF pairs have both m\2 < maxm-12. In 90% of the events, there is only one combination 
that satisfies this requirement. This allows us to know the branch assignment of 95% of the four 
lepton events without ambiguity. 

Unlike the three-body decay case, the order of the two leptons on each branch matters. The 
intermediate mass M x = («2 + p) 2 depends on «2 and does not depend on a\. To resolve this 
ambiguity we consider the four combinations that preserve the branch assignment but differ 
in their ordering. The M%c,lb for the event is the minimum of these combinations. Likewise 



the M^c,ub is the maximum of these combinations. As one expects, Fig. 8.2 (Left) shows how 
the combinatorics ambiguity degrades the sharpness of the cut-off at the true mass. Not all 
applications share this ambiguity; for example in top-quark mass determination (pair produced 
with Y = top, X = , N = v) the 6-quark-jet marks a\ and the lepton marks 0.2 ■ 

• Effect Large Upstream Transverse Momentum 
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Ideal Combinatorics Vs RealisticCombinatorics 
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Figure 8.2: (Left) The M^c distributions before (solid) and after (dashed) introducing the 
combinatoric ambiguity. (Right) The M$c distributions with and without UTM. The no UTM 
case {kx = 0) is shown by the solid line; the large UTM case with kx = 250 GeV is shown by 
the dashed line. 

In a similar behavior to M2C, the distributions of the variable M^c show a strong dependence 
on large upstream transverse momentum (UTM). In our case study this is identified as the 



combination of all the hadronic activity. Fig 8.2 (Right) shows the stronger upper-bound cut- 



off in the presence of large UTM. Unlike M2C, in M^c with kx = we still have events with 
non-trivial upper bound values. 



We also tested the distribution for different values of k 2 . In Fig 8.2 (Right) we fixed k 2 = 
(100 GeV) 2 . We also performed simulations with k 2 = (500 GeV) 2 and found the difference 
of the two M3C distributions consistent with zero after 15000 events. In other words, the 
distribution depends mostly on k x and k y and appears independent of k^. 

• The Effects of Detector Energy Resolution 

Compared to M2C1 the information about the extra states gives a stronger set of bounds. 
Unfortunately, the solution is also more sensitive to momentum measurement error. We model 



the finite energy resolution using Eqs(7.13 7.15). The hadronic energy resolution is larger than 



the leptonic energy resolution which will increase the uncertainty in the missing transverse 
momentum. 

Fig. |8.3| shows the effect on the M%c distribution of realistic leptonic energy resolution while 
keeping kx = 0. On the left we show the energy resolution (dashed line) compared to the perfect 
energy resolution (solid line). The energy resolution washes out the sharp cut-off. On the right 
we show M3C.LB with perfect energy resolution plotted against the result with realistic energy 
resolution. This shows that the cut-off is strongly washed out because the events with M%c closer 
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Figure 8.3: The effect of energy resolution on the M%c distribution. (Left) The dotted line 
shows the energy resolution has washed out the sharp cut-off. (Right) M^c,lb with perfect 
energy resolution plotted against the result with realistic energy resolutions. 

to the true value of My (200 GeV in this case) are more sensitive to energy resolution than the 
events with M^c closer to AMyn- The peak in the upper-bound distribution at M^c,ub = 100 
GeV comes from events that no longer have solutions after smearing the four momenta. 

Because the energy resolution affects the distribution shape, its correct modeling is impor- 
tant. In the actual LHC events the energy resolution will depend on the hadronic activity 
in the events being considered. Two events with the same kx = may have drastically different 

resolutions. Modeling the actual detector's energy resolution for the events used is impor- 
tant to predict the set of ideal distribution shapes which are compared against the low-statistics 
observed data. 



• Parton Distribution Function Dependence 

For a mostly model-independent mass determination technique, we would like to have a 
distribution that is independent of the specific production mechanism of the assumed event 
topology. The parton distributions determine the center-of-mass energy ^J~s of the hard collisions; 
but the cross section depends on model-dependent couplings and parameters. The events we 
consider may come from production of different initial states (gluons or squarks) but end in 
the assumed decay topology. The M%c distribution, like the M2C distribution, shows very little 
dependence on the underlying collision parameters or circumstances. 



Fig 8.4 (Left) shows the dependence of the M%c distributions on the parton collision en- 
ergy. The solid line shows the M%c distributions of events with collision energy yfs distributed 



according to Eq(7.10), and the dashed line shows the M^c distributions of events with fixed 
yfs = 600 GeV. Fig 8.4 (Right) shows the difference of these two distributions with 2a error 
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Figure 8.4: The dependence of the M^c distributions on the parton collision energy. The solid 



line shows the collision distributed according to Eq(7.10), and the dashed line shows the collision 
energy fixed at yfs = 600 GeV. 

bars as calculated from 15000 events. The two distributions are equal to within this numerical 
precision. 

• Effects of f T Cuts 



As described in |187| I186| [T5| [16], the Standard Model four leptons with missing transverse 
momentum backgrounds are very strongly suppressed after a missing transverse momentum 
cut. This requires an analysis of what type of an effect a fj* > 20 GeV cut will have on the 



distribution shape. Fig. 8.5 shows that the effect of this cut is dominantly on the smallest M%c 
bins. On the left we see the M%c lb result versus the frj*. Unlike the M2C case in Fig 



7.8 



the M3C solutions in Fig. 8.5 do not correlate with the fj*. The right side of the figure shows 
difference between the M%c,ub and M%c,lb distributions with and without the cut > 20 
GeV. The smallest bins of M^c,lb are the only bins to be statistically significantly affected. 
The left-side suggests this lack of dependence on fij, cuts is somewhat accidental and is due 
to the nearly uniform distribution of M^c solutions being removed by the cut. The stronger 
dependence of the smallest M%c bins on the fx cut means we can either model the effect or 
exclude the first bins (about 10 GeV worth) from the distribution used to predict the mass. We 
will choose the latter because we will find that the background events also congregate in these 
first several bins. 

• Spin correlations 

In our simulation to produce the ideal curves, we assumed each decay was uncorrelated with 
its spin in the rest frame of the decaying particle. Spin correlations at production may affect 
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Figure 8.5: The effect of missing transverse momentum cuts on the M^c distributions. (Left) 
The M^clb result versus the (Right) The difference of the M^c,ub and M^clb distribu- 
tions with and without the cut > 20 GeV. The smallest bins of M^c,lb are the only bins to 
be statistically significantly affected. 

this; however, such spin correlations are washed out when each branch of our assumed topology 
is at the end of longer decay chain. These upstream decays are the source of considerable UTM. 

Some spin correlation information can be easily taken into account. The myi (or 77134) 
distribution's shape is sensitive to the spin correlations along the decay chain |15(J] . The observed 
mi2 (or 77134) distribution can be used as an input to producing the ideal distribution shape. In 
this way spin correlations along the decay chain can be taken into account in the simulations of 
the ideal distributions. 

Spin correlations between the two branches can also affect the distribution shape. To demon- 
strate this, we modeled a strongly spin-correlated direct production process. Fig|8.6|(Left) shows 



the spin-correlated process that we consider. Fig 8.6 (Right) shows the M^c upper and lower 



bound distributions from this process compared to the M%c distribution from the same topol- 
ogy and masses but without spin correlations. We compare distributions with perfect energy 



resolution, m v 



m,- 



GeV, M Y = 200 GeV, M x = 150 GeV, and M N = 100 GeV. Our 



maximally spin-correlated process involves pair production of Y through a pseudoscalar A. The 
fermion Y in both branches decays to a complex scalar X and visible fermion v\ through a 
purely chiral interaction. The scalar X then decays to the dark-matter particle N and another 
visible particle i>2- The production of the pseudoscalar ensures that the Y and Y are in a config- 
uration y/2 1 (\ ||) + I II)). The particle Y then decays with X preferentially aligned with the 
spin. The Y decays with X* preferentially aligned against the spin. Because X is a scalar, the 
particle N decays uniformly in all directions from the rest frame of X. The correlated directions 
of X causes the two sources of missing transverse momentum to be preferentially parallel. The 
resulting greater magnitude of missing transverse momentum increases the cases where M^c has 
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Figure 8.6: Effect of a spin correlated process on the M^c distributions. Modeled masses are 
My = 200 GeV, M x = 150 GeV, and M N = 100 GeV. The thick black and thick blue lines 
show the distributions of the uncorrelated lower bound and upper bound M^c- The dotted red 
lines show the distributions of the spin correlated process. 

a solution closer to the endpoint. For this reason the spin correlated distribution (red dotted 
distribution) is above the uncorrelated distribution (black thick lower bound distribution and 
blue thick upper bound distribution). The upper bound distribution is statistically identical 
after 25000 events. The lower bound distribution clearly has been changed, but not as much as 



the M2C distribution in Fig. 6.3 This is due to the subsequent decay of the X particle which 
lessens the likelihood that the two iVs will be parallel. For the remainder of the chapter we 
assume no such spin correlations are present. 

• Backgrounds 

The Standard Model backgrounds for four-leptons and missing transverse momentum are 
studied in |187] 1186] . In the two previous chapters we summarized the SM backgrounds and the 
dominance of SUSY backgrounds for this channel (see also [I5][l6]). As was mentioned earlier, 
the Standard Model backgrounds for four leptons with missing transverse momentum are very 
strongly suppressed after a missing transverse momentum cut. 

To improve the quality of the fit, a model for backgrounds can be created based on as- 
sumptions about the origin of the events and wedge-box analysis like those described in Bisset 
et.a/[160 and references therein. We preformed such a model in the previous chapter, and found 
the distribution shape isolated the correct mass of x° to within 1 GeV with versus without the 
background model. Although the quality of the fit without modeling the backgrounds is de- 
creased, we find that the mass of the LSP associated with the best-fit is tolerant to unknown 
backgrounds. In the SPS la example studied in the next section, SUSY background events form 
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about 12% of the events. We again see less than a 1 GeV shift in Af^° with versus without the 
background events. As such, we do not try to model the background in this M%c study. 

8.4 Estimated Performance 

With an understanding of the factors affecting the shapes of the M%c distributions, we combine 
all the influences together and consider the mass determination performance. We follow the 



same modeling and simulation procedures used in Sec. 7.2 except now we include an on-shell 
intermediate state and calculate M%c- We use HERWIG [T71 [181 EH] to generate events according 
to the SPS la model [3 . This is an mSUGRA model with m = 100 GeV, m 1/2 = 250 GeV, 
A = —100 GeV, tan/3 = 10, and sign(/i) = +. We initially assume the mass differences 
AMyn = 80.8 GeV and AMxn = 47.0 GeV have been previously measured and take them as 
exact. We later show how the distribution shape with the mu endpoint also solves for the two 
mass differences. 

Like M2C1 t ne distributions is able to be well-predicted from observations. When we 
are determining masses based on distribution shapes, the larger the area difference between two 
distributions representing different masses, the more accurately and precisely we will be able to 
tell the difference. Unfortunately, the M^c distribution is sensitive to the finite momentum res- 
olution errors and combinatoric errors which have the effect to decrease the large area difference 



between the distributions of two different masses shown in Figs. 8.1 

Just as in Chapter [7j we model the distribution shape with a simple Mathematica Monte 
Carlo event generator, and compare the predicted distribution shapes against the HERWIG data 
modeling the benchmark point SPS la. We again use the observed UTM as an input to the 
Mathematica simulated ideal distributions. By modeling with the Mathematica, which does 
not use the SUSY cross sections, and comparing to more realistic HERWIG generated data, we 
hope to test that we understand the major dependencies of the shape of the M$c distributions. 
The Mathematica event generator produces events based on assumptions uniform angular dis- 
tribution of the parents in the COM frame, the parent particles decay with a uniform angular 
distribution in the rest frame of the parent. The particles are all taken to be on shell, > 
is simulated by boosting the event in the transverse plane to compensate a specified kj*. 



Figure 8.7 shows the performance. The left side of Fig. 8.7 shows the M%c lower bound and 
upper bound counts per 5 GeV bin from the HERWIG generated data, and it shows the predicted 
ideal counts calculated with Mathematica using the observed UTM distribution and assuming 
My-o = 95 GeV. The upper bound and lower bound show very close agreement. The background 
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Figure 8.7: Fit of ideal M%c{M^°) distributions to the HERWIG generated M^c distributions. 
Includes combinatoric errors, backgrounds, energy resolution, and fx cuts. (Left) The observed 
HERWIG counts versus the expected counts for ideal = 95 GeV. (Right) The x 2 fit to ideal 
distributions of A% = 80, 85, 90, 95, 100, 105, 110 GeV. The correct mass is = 96.0 GeV. 

events are shown in dotted lines and are seen accumulating in the first few bins. These are the 
same bins dominantly affected by fx cuts. For this reason we excluded these first two bins 
from the distribution fit. The right side of Fig 8.7 shows the x 2 fit of the HERWIG simulated 



data M3C distribution to the ideal M^ciM^) distribution with taken as the independent 
variable. Ideal distribution shapes are calculated at values of = 80, 85, 90, 95, 100, 105, 110 
GeV. The x 2 fitting procedure is described in more detail in Appendix [dJ All effects discussed 
in this chapter are included: combinatoric errors, SUSY backgrounds, energy resolution, and 
fx cuts. Our ideal curves were based on the Mathematica simulations with 25000 events per 
ideal curve. Despite the presence of backgrounds, the x 2 1S n °t much above 1 per bin. 



The particular fit shown in Fig 8.7 gives = 98.6±2.2 GeV where we measure uncertainty 



by using the positions at which x 2 lS increased by one. We repeat the fitting procedure on nine 
sets each with ~ 100 fb" 1 of HERWIG data (~ 1400 events for each set). The mean and standard 
deviation of these nine fits give = 96.8 ± 3.7 GeV. After 300 fb -1 one should expect a \/3 
improvement in the uncertainty giving ±2.2 GeV. The correct mass is = 96.05 GeV. 

Our technique also enables a combined fit to both the mass differences and the mass scale. 



The mu endpoint in Eq(4.3) constrains a relationship between the three masses. Gjelsten, 



Miller, and Osland estimate this edge can be measured to better than 0.08 GeV \157\ 1159] using 
many different channels that lead to the same edge, and after modeling energy resolution and 
background. In the next several paragraphs we show that by combining this edge with the fits to 
the M3C upper bound and lower bound distribution shapes, we can constrain all three masses. 
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Figure 8.8: Combined constraint from fitting both maxm« and M^c with the mass difference 
as free parameters. We parameterized the difference from the true values in the model by 
AMyn = 80.8 GeV + 5AM YN and AMxn = 47.0 GeV + SAM XN - We shown the 1,2,3a 
contours. 

We first numerically calculated the effect of errors in the mass differences. We used 300 fb _1 
of events (about 3600 signal and 450 background) including all the effects discussed. We pa- 
rameterize the error from the correct mass difference in the model by the variables 5 AMyn an d 
SAMxn so that mass differences are given by AMyn = 80.8 GeV + 5AMyn an d AMxn = 
47.0 GeV + 5 AMxn ■ We calculated the Xm a t 8 points surrounding the correct mass difference 
by amounts 5AMyn = if GeV and 5AMxn = if GeV. The minimum Xm 3C a t each of the 
9 points gives the value of for each mass difference assumed. The position of the minima 
can be parameterized by a quadratic near the true mass difference. The resulting fit 

M n = 96.4 + 1.9 (5AM X n) 2 + 2-5 5AM YN SAMxn + 3-2$AMxn -3.8 (5AM YN ) 2 -8.3 5AM YN 

(8.15) 

shows in units of GeV how the mass is affected by small errors in the mass difference. 

The Xm 3C , a t these 9 different values for the mass difference provides another constraint on 
the mass differences. Fitting the Xm to a general quadratic near the true mass difference gives 



x 2 Mic = 162 + 38 {5 AMxn) 2 - 8 5AM Y n5AM X n - 5 6AM X n - 25 5AM YN - 



3.16) 



The Xm described by Eq(8.16) shows a sloping valley. The sides of the valley constrain 



5AMxn &s seen by the large positive coefficient of (6 AMxn) 2 ■ The valley slopes downward 



119 



along 5AMyn as can be seen by the large negative coefficient of 5AMyn which leaves this axis 
unbounded within the region studied. 

The unconstrained direction along AMyn can by constrained by the mass relationships 
given by the endpoint maxmu or by Mt2 as described in Chapter [5] Here we work with 
max.mii to provide this constraint. We calculate the Xmaxmu usm S Sm.ax.mu = 0-08 GeV, and 



Eq(4.3) with My = AM Y n + SAM Y n + A% and M x = AM X n + SAM X n + where we 
use My.° from Eq(8.15). This x maxrnn constrains a diagonal path in (5 AMyn, SAMxn)- The 
value of the Xmaxm;; at the minimum is a constant along this path. The combined constraint 
^M 3C t° Xmaxmu leads to the a minimum at SAMyn = 0.18 GeV and SAMxn = 0.25 GeV 



where = 95.7 GeV as shown in Fig. 8.8 We have shown the contours where x 2 increases 
from its minimum by 1,2 and 3. The uncertainty in the mass differences around this minimum 
is about ±0.2 GeV. The bias from the true mass differences is due to the unconstrained Xm 3C 
along SAMyn- We can use modeling to back out the unbiased mass differences. Propagating the 
effects of uncertainty in the mass differences, we estimate a final performance of = 96.4±2.4 
GeV after 300 fb -1 with about 3600 signal events amid 450 background events. We find the mass 
differences (without bias correction) of M^-M^o = 81.0 ±0.2 GeV and Mf -M^o = 44.3 ±0.2 
GeV. This is to be compared to the HERWIG values of A% = 96.0 GeV, M%> - = 80.8 GeV, 
and Mr - Mc° = 44.3 GeV. 

l R X-l 

How does this performance compare to other techniques? Because SPS la is commonly 
used as a test case, we can approximately compare performance with two different groups. The 
technique of [156, 1157] 1158} H"59] which uses edges from cascade decays determines the LSP mass 
to ±3.4 GeV with about 500 thousand events from 300 fb -1 . The approach of CEGHM [W] 
assumes a pair of symmetric decay chains and assumes two events have the same structure. 
They reach ±2.8 GeV using 700 signal events after 300 fb™ 1 , but have a 2.5 GeV systematic 
bias that needs modeling to remove. Both techniques also constrain the mass differences. By 
comparison we find ±3.7 GeV after 100 fb -1 (1200 signal, 150 background) and estimate ±2.4 
GeV after 300 fb -1 (3600 signal, 450 background) and propagating reasonable uncertainties in 
the mass differences. The uncertainty calculations differ amongst groups. Some groups estimate 
uncertainty from repeated trials, and others use the amount one can change the mass before 
X 2 increases by one. Without careful comparison under like circumstance by the same research 
group, the optimal method is not clear. What is clear is that fitting the M^c and M^c,ub 
distributions determines the mass of invisible particles as well if not better than the other 
known methods in both accuracy and precision. 



120 



Chapter Summary 



In this chapter, we have extended the constrained mass variable to the case with three new 
particle states. We assume events with a symmetric, on-shell intermediate-state topology shown 



in Fig. 4.3 We can either assume that we have measured the mass difference between these 
new states through other techniques, or combine our technique with the maxm^ edge to find 
both mass differences and the mass scale. The new constrained mass variables associated with 
events with these three new particle states are called M^c,lb and M^cub-, and they represent 
an event-by-event lower bound and upper bound (respectively) on the mass of the third lightest 



state possible while maintaining the constraints described in Eqs(8.2 8.7). We have shown that 
most of the Mic distribution properties described in the previous chapter carry through to M^c- 
The additional particle state and mass difference enable a tighter event-by-event bound on the 
true mass. The M%c distribution is more sensitive than the Mic distribution to the momentum 
and energy resolution errors. Studying the performance on the SPS la benchmark point, we 
find that despite the energy resolution degradation, we are able to determine to at least the 
same level of precision and accuracy or possibly even better precision and accuracy than that 
found by using cascade decays or by using other MSTs. 
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Chapter 9 

Discussions and Conclusions 



In this thesis, we have started in Chapter [2] with a study of principles that in the past have 
successfully predicted the existence and the mass of new particle states. We showed astrophysical 
evidence for dark matter and discussed properties that make supersymmetry an attractive theory 
for the dark-matter candidate. We discussed the grand unification of the three gauge couplings, 
and we introduced the mass-matrix unification suggested by Georgi and Jarlskog. In Chapter 
[3] we showed how the mass unification is being quantitatively challenged by tighter constraints 
on the parameters, but that the mass unification hypothesis leads to a favored class of tan/? 
enhanced threshold corrections where tan (3 is large and the sign of the gluino mass is opposite 
that of the wino mass. Chapter [4] turns to the task of measuring the masses of new particle states 
if they are produced with dark matter at hadron colliders. We observed that determining the 
sign of M3, determining spin, and determining the entire sparticle mass spectra are all facilitated 
by a model-independent measure the mass of the dark-matter particle. We discussed the current 
techniques to constrain the masses of the new particle states. We developed ideas of how to 
use Mt2 in Chapter [5] which enables us to find new constraints between particles in symmetric 
cascade decays, and we argued that Mt2 is a better variable to extract this information than 
Mct- Chapters [6] and [7] discussed the constrained mass variables M2C and M2c,ub which uses 
a previously measured mass difference to find the maximum and minimum mass of the dark- 
matter particle on an event-by-even basis. This technique benefits from an M2C and M2CUB 
distribution that can be determined from experimental observables and the unknown dark- 
matter mass. The distribution shape can be fit to the bulk of the data to estimate the mass of 
the dark-matter particle more accurately than end-point techniques. Chapter [8] discussed the 
benefits and down- falls of the constraint mass variable M^c- Unlike the edges- from-cascade- 
decays method |156| 1157, 158, 1159] and the mass-shell technique of CEGHM |167] each of which 
require the existence of four new particle states, the constrained mass variables only require two 
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new particle states for M2C and three new particle states for M$c- Although proper comparisons 
to other techniques are difficult, it is clear that the constrained mass variables can determine the 
mass of invisible particles at least as well if not better than the other known model-independent 
methods in both accuracy and precision. 

From here what are the directions for continued research? The origin of the Yukawa mass 
matrices remains unknown. The Georgi-Jarlskog relations appear to hint at a solution, but only 
decrease the number of unknowns slightly. The fits to the GUT-scale Yukawa parameters is a 
starting point for more research on explaining this structure. There are relatively few studies 
on how to measure the sign of the gluino mass parameter relative to the winos mass parameter; 
there are even fewer that do not make use of simplifying assumptions about the underlying 
supersymmetry model. The results of Chapter [5] need to be extended to include the effects of 
combinatoric ambiguities. In Chapters [7]{8] there are several distribution properties which should 
be better studied. We would like to understand the physics of the invariance of the constrained 
mass variable distributions to changes in y 7 ^ and choice of k 2 when ^ / 0. In Chapter [8] it 
would also be nice to classify the conditions which determine how many islands of solutions exist 
for M3C. The M%c method could be tested on top-quark data to remeasure the masses of the 
t, W and v. In Chapters [6]{8] we have considered two constrained mass variables M^c and M^c 
where the decay chains are symmetric, but we may find greater statistics for asymmetric events. 
For example in many models one expects more events where X2 ^ s produced with xi 01 X ■ 
The concepts of the constrained mass variable could be extended to deal with these asymmetric 
event topologies or to cases with both LSPs and neutrinos. 

In conclusion, the thesis contributed to predictions of mass properties of new particle states 
that may be discovered at the LHC and developed the tools to make precision mass measurements 
that may help falsify or validate this prediction and many others. Modern physics is rooted in 
an interplay between creating theoretical models, developing experimental techniques, making 
observations, and falsifying theories. The LHC's results will complete this long-awaited cycle 
making experimental observations to constrain the growing balloon of theories of which this 
thesis is a small part. 
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Appendix A 

Renormalization Group Running 



This appendix details our use of the Renormalization Group Equations (RGEs) to relate the 
reported values for observables at low-energy scales to their values the GUT scale. The one-loop 
and two-loop RG equations for the gauge coupling and the Yukawa couplings in the Stan- 
dard Model and in the MSSM that we use in this thesis come from a number of sources 
|113] [40 [114 [110 . We solve the system using the internal differential equation solver in Mathe- 
matica 6.0 on a basic Intel laptop. We considered using SoftSUSY or SuSpect to decrease the like- 
lihood of programming errors on our part, but neither met our needs sufficiently SoftSUSY[191 
has only real Yukawa couplings, and SuSpect |183j focuses primarily on the third generation RG 
evolution. 

We assume a model for the neutrino spectrum that does not impact the running of the quark 
or charged lepton parameters below the GUT scale. This assumption is not trivial; there are 
neutrino models where Y u will affect the ratio yb/Ur by as much as 15% |192j . Furthermore 
due to unitarity constraints, the UV completion of the effective dimension-five neutrino mass 
operators must take place before the GUT scale |193j . One way to have a UV completion below 
the GUT scale while making the RGE running independent of the specific model is to set-up a 
non-degenerate right-handed neutrino spectrum. One then constructs the couplings such that 
large Yukawa couplings in Y v are integrated out with the heaviest right-handed neutrino mass 
which is at or above the GUT scale. One can then have the remaining right-handed neutrinos 
considerably below the GUT scale which couple to much smaller values Y u that do not affect 
the running of Y e or Y u . The rules for running the neutrino parameters and decoupling the 
right-handed neutrinos for the Umns are found in Refs. |113||194] |195] |196] |197| . For our 
work finding the quark and lepton masses and mixings at the GUT scale, we assume the effects 
of the neutrinos on the RGE equations effectively decouple at or above the GUT scale. 

Our running involves three models at different energy scales. Between Mx to an effective 
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SUSY scale Ms, we assume the MSSM. Between Ms to Mz, we assume the Standard Model. 
Below Mz, we assume QCD and QED. We do not assume unification, rather we use the low- 



energy values in Table 3.1 and Eq(3.1) provide our boundary condition. The uncertainty in the 
final digit(s) of each parameter is listed in parenthesis. For example 123.4 ± 0.2 is written as 
123.4(2). 

First, we use QCD and QED low-energy effective theories to run the quark and lepton masses 
from the scales at which the particle data book reports their values to Mz- The Vqkm mixing 
angles are set at this scale. Next, we run using the Standard-Model parameters, including the 
Higgs VEV and Higgs self coupling from Mz to an effective SUSY scale Ms- At an effective 
SUSY scale Ms, we match the Standard Model Yukawa couplings using the running Higgs VEV 
v{Ms) onto MSSM Yukawa couplings. In the matching procedure, we convert the gauge coupling 



constants from MS to DR as described in [191J, and we apply approximate supersymmetric 
threshold corrections as parameterized in Chapter [3] Details of the RG equations and the 
matchings are described in the following sections. 

A.l RGE Low-Energy SU(3) c xU(1) E m up to the Standard Model 

Below Mz we use QCD and QED RG Equations to move the parameters to the scale reported 
in the Particle Data Book. We run the light-quark masses between = Mz to fii = 2 GeV 
with the following factors: 

Vu,d,s = (1-75 ± 0.03) + 19.67 («i 5) (M z ) - 0.118) (A.l) 
r/ c = (2.11 ± 0.03) + 41.13 (af\M z ) - 0.118) (A.2) 
i] b = (1.46 ± 0.02) + 8.14 (af\M z ) - 0.118) , (A.3) 



that are calculated with the 4-loop RGE equations using Chetyrkin's RUNDEC software for 
Mathematica |198j . These factors are used as m^rrib) /rjb = m&(Mz). The uncertainty in the 
first term is an estimate of the theoretical error from neglecting the five-loop running give by 



5r]/rj ~ exp(C(l)(a s ) 5 log —^-) — 1. RUNDEC also converts the top quark's pole mass to an MS 
running mass, and applies the small threshold corrections from decoupling each of the quarks. 
The different rj factors take the parameters from a six-quark effective model at [i = Mz to their 
respective scales. 

The leptons can be converted from their Pole mass into running masses at Mz by a similar 
running factor: 

rj e = 1.046 rjf, = 1.027 rj T = 1.0167. (A.4) 
125 



Parameter 


At Scale fi = M z 


m t 


169.6(2.3) GeV 


m c 


599(50) MeV 


m u 


1.54(50) MeV 


m b 


2.87(6) GeV 


m s 


57(12) MeV 


m d 


3.1(5) MeV 



Table A.l: The MS values for the running quark masses at Mz- 
We did not incorporate the lepton masses as error sources. 



The input values in Table 3.1 lead to the values at Mz in Table A.l where we have propagated 
uncertainty from both the strong coupling constant uncertainty and the uncertainty in the mass 
at the starting scale. These values for the parameters at Mz agree with those calculated recently 
elsewhere in the literature |199| . 

The gauge coupling constants g\ and are determined from egM, sin 2 9w at the scale Mz 

by 



gt(M 2 



5 e 2 EM (M 2 



gl{M z ) 



-EM 



{M z ) 



(A.5) 



3 cos 2 6w ^' sin 2 6w 

A.2 RGE for Standard Model with Neutrinos up to MSSM 

We found using the two-loop vs the one-loop RG equations for the gauge-coupling produced a 
shift in the low energy parameters of more than la. For this reason, we always use two-loop RG 
equations for the gauge couplings. For the Yukawa couplings we find the two-loop vs the one-loop 
RGE equations shift the low-energy parameters much less than the experimental uncertainty. 
Therefore, we perform our minimizations using the one-loop RG equations for the Yukawa and 
checked the final fits against the two-loop RG equations. 

We reproduce here for reference the two-loop gauge coupling RG equations and the one-loop 
Yukawa coupling RG equations as this level of detail is sufficient to reproduce our results. We 
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define t = logM. The equations are compiled from [196 , (1 13j . |110| . 



d _ 1 41 3 
It 91 ~ 16V 2 Yo 91 



9l 



(16vr 2 )< 



d _ 1 -19 3 

Jt 92 ~ m^~^ 92 



d -7 3 

dt 93 ~ 16^ 93 



(16vr 2 ) 2 VlO 

i 

9% ( I I 



g\ + | g\ + 12 g\ - ^Tr(Y u Y^) - ^Tr(Y d Y^) - l Tr(Y £ Y^.7) 



(16vr 2 ) 2 VlO 



d y V 1 / — ,„, , „ 3 



dt 16vr 2 V 2 
d yv 1 ( _„ 3. 



dt 16vr 2 V 2 



where 



16vr 2 V 2 

-F e = (-IG e + IT+ -Y e Y e ^ 

dt 16vr 2 V 2 



r 1 " — 17 n 2 , 9 „2 , q 2 /^d _ 1 „2 , 9 „2 , o n 2 

G — 20 9l + 4 52 + ° 93 G — 4 5l + 4 52 + 8 53 
rje _ 9 n 2 , 9 „2 ™ _ 9 _2 , 9 J2 

4 5l ' 4 ^2 u — 20 5l "T 4 52 



2Tr(y d y dt ) 


) 


(A.8) 


+ yyt) 




(A.9) 


3 yd Y d ^ 

2 y 


| .y u 


(A.10) 


_ 3 ye ye A 

2 y 




(A.11) 


y M y"^ 

2 J 


| .Y d 


(A.12) 


3 y^ y^A 
2 ) 


i .y e 


(A.13) 



9 „2 i 9 „2 nv — 9 „2 , 9 „2 

(A.14) 



The Higgs self-interaction A (following the convention —X(H^H) 2 C C) obeys the RG equa- 
tion: 

16tt 2 |a = i2A 2 -(^ 2 + 9 9 i)A + ^(A 5 4 + 2 5 2 5 2 +5 4 )+4TA (Al5) 
— 4Tr(3 y u y"ty«y«t _|_ 3 ydydfydydf _|_ yeyetyeyef\ (A. 16) 

Between Mz and M5, we run the Standard-Model VEV using |40] : 

.2 '/ .. .. f-K 1 .;. , ..2, 



16 ^ V = " U ( 5 51 + ^ " T ' ' (A ' 17) 

167r 2 4^ = (--Y e Y^K--KY e Y e ^ + -Y u Y^K + -KY v Y v ^ 
dt \ 2 2 2 2 

+1 (tt(2 y v y^ + 2 y e y et + 6rrt + 6y d y dt ) - 3 5 | + 4a) « 
1 M ^= 1^2 ((^n'^ + Miiiirtr) . (A.18) 

We normalize g\ to the SU(5) convention: g\ = |<7y. 
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A.3 RGE for the MSSM with Neutrinos up to GUT Scale 



The MSSM RG equations are in the DR scheme. To convert MS running quark masses to DR 
masses at the same scale, we use 

*- = = Ba =fi-i4-S(A) , V (A.w) 



mwB y 3 4vr 2 72 v 4vr 2 ' 

No corrections were used for switching lepton running masses from one scheme to the other. 
The gauge coupling constants are related via [191J 

«3E5 = «sk + (A - 20) 

where 03 = g\j^ and 02 = g\j^. We define t = logM. The RG equations are: 
d 1 33 3 

+ + T 9 * + f 9 * ~ t^XTy*') ~ l 4^ YdYd] ) ~ yTr(y T A^) 

d 1 3 

dt 92 - 16^ ff2 



+ , 9 \, (-gl + 2bgl + 2Ag 2 - 6Tr(y"y ut ) - 6Tr(Y d Y d ^) - 2Tr(Y e Y et ) ) (A.23) 
(16tt z ) z \5 ' 



d — 3 3 

dt 93 = 16^ 93 



+ 9 \ ( — g\ + %g\ + 14 ff | - 4Tr(y"y ut ) - 4Tv(Y d Y d ^) ) (A.24) 



(16^ 2 ) 2 V 5 

jY u = (^-IG U + I3Tr(Y u Y u ^)+ITv(Y u Y u ^)+3Y u Y^ + Y d Y d ^.Y\A.25) 

jY u = ^— I G u + 13 Tr(y M y"^) + I Tr(Y u Y v ^) + 3Y U Y U ^ + Y e Y e ^ .Y v (A. 26) 

^-Y d = ^-^(^-IG d + I3Tr(Y d Y d ^)+ITr(Y e Y^)+3Y d Y d ^ + Y u Y u ^.Y d (A.27) 
jY e = (-IG e + I3Tr(Y d Y dt ) + ITr(Y e Y et ) + 3Y e Y et + Y U Y V ^ .Y e (A.28) 



where 



_ 13 „2 , o n 2 , 16 „2 _ 17.2 ,o-2 1 16 „2 

G e = | 5 f + 3 3 1 G^ = | 5 2 + 3 5 2 . 



(A.29) 



16tt 2 ^k = ( Y e Y e ^K + K Y e Y et + YT^k + ^"7"* 
dt \ 

+1^2 IrfFTt + 6Y u Y ut ) - - 65^ k (A.30) 
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Figure A.l: The impact of RG running parameter ratios with Ms = 500 GeV. These ratios 



determine x defined in Eq. A. 34 If Ms = Mz, all three are degenerate at small tan/?. 



j t M RR = ^2 (2 (r 4 y v ) t M RR + 2 M RR y v 



(A.31) 



We normalize g\ to the SU(5) convention: g\ = |<7y. 

A. 4 Approximate Running Rules of Thumb 

At tan/3 < 10 and when Ms ~ Mz, the running of parameters from Mz to Mx obey the 
following relationships exploited in the RRRV |L study: 



flo(M x ) 
f)o{M z ) 



po(M X ) 
po{M z ) 



Xo(M x ) 
Xo(M z ) 



1 



Aq(Mx) 
Ao{M z ) 



X 



(m u /m c ) (M x ) _ (m d /m s ) (M x ) 
(m u /m c ) (M z ) 

where 



(m d /m s ) (Mz) 



1 



(m s /m b ) (M x ) 
(m s /m b ) (M z ) 



(m c /m t )o{M x ) 
(m c /m t )o(Mz) 



X ~ exp 



to 



-Vt 
16vr 2 



dt 



(A.32) 



X 3 (A.33) 



(A.34) 



and t a = log Mz and tj = log Mx- Fig. A.l shows the two-loop results reflecting how the running 
of these ratios change when we consider Ms = 500 GeV and larger values of tan j3. If we had 



plotted the Ms = Mz case, the curves in Fig. A.l would have been completely degenerate at 



small tan (3. As a consistency check, we compare our results to Fusaoka and Koide jl 10] . They 
find x ~ 0.851 at tan/3 = 10 with their choice of mt{Mz) = 180. Our code also gives x = 0.851 
at tan (5 = 10 if we select Ms = Mz and omit the conversion to DR scheme. 
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Appendix B 

Hierarchical Yukawa couplings and 
Observables 



In this appendix, we find the masses and unitary matrices associated with a general hierarchical 
Yukawa coupling matrix. We use this expansion to find general expressions for Vckm to the 
order needed for comparison to current experiments. Consider a general hierarchical Yukawa 
matrix 

AT! *?s\ fO(e A ) 0(e 3 ) 0(e 3 )\ 

Y u = y t \Y£ Y 2 \ Y 2 l] « Y u = y t 0(e 3 ) 0(e 2 ) 0(e 2 ) (B.l) 
Wi 1 / \0(e 3 ) 0(e 2 ) 1 ) 

where e is a small parameter. Diagonalization of Y u Y u ^ leads to the diagonal matrix from 

jjuyuyu] = \D U \ 2 where Y u = D U U^. The matrix \D U \ 2 gives the square of the mass 

eigenstates. 

\D%/y t \ 2 = l + Y 2 u 3 Y 2 u * + Y-Y-* + 0(e 5 ) (B.2) 
\D2/Vt\ 2 = Y&Yg + Y&Yg + Y? 1 Y£-YgY%Y&-Y&Y%'Y£' + 0(e 7 ) (B.3) 

\ru* \ru yu yti yw \ru* \ru \ru* yu yw 

I n« /„, |2 _ y« y«* Ml 1 12 *21 1 11 J 12 J 21 , -* 12*12 J 21 *21 , AyaO\ fu A\ 

\ U ZlVt\ - ^11 ^11 77^ y^ I y« y U * <~ U l £ J I 13 - 4 ) 

2 22 2 22 2 22 2 22 

The resulting approximate expressions for are 

(t/£t)i3 = Y« + y« Y 3 - + 0( e 6 ) (B.5) 

([/«t )2 3 = yu + yu yu* + 0(£6) (B 6) 

([/ £ t) 33 = ! _ 1 yuyu* + ^6) (R?) 
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These results have been checked numerically including the phase. 



If we assume the Y d Yukawa matrices follow the same form as Eq. B.l we can obtain parallel 
expressions for Iff. The CKM matrix is then given by Vckm 
expressions for the individual components: 
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In the above equations, we have included enough terms to achieve experimental accuracy for 
the Yukawa couplings. 

For completeness, we include the diagonal entries 

\rd \rd* xrd* v u v u V u * 

Y 1 12 12 i 12 1 12 2 12 2 12 fB 20"l 

ud o vd v d * v d * vu 9 V u * ^ ' ' 

z 2 22 2 22 2 22 2 22 22 2 22 

v^d "v/"d* x^rf \ru* v u \ru* 

y 1 _ 12 12 _|_ l 2 12 _ 2 12 I 12 /g 21) 

cs o vd* vd vu* o V"" V u * ^ ' ' 

^ 2 22 2 22 2 22 2 22 ^ 2 22 2 22 

% = i + y 2 <|y--^3_-^3„ (B.22) 

To express the Wolfenstein ^4 and A terms to experimental accuracy, one needs the lengthy 
expression for V us . Because the Wolfenstein parameters are redundant with the Vckm elements 
from which they are based, we will omit explicit expressions for A and A. The Wolfenstein 
parameters p + ifj represent a phase convention independent expression for the CP violation and 
can be expressed compactly: 

. \rd* \ru* \*u* \su* 

\^d* v u * J 23 12 _i_ 1 12 J 23 

1 1 1 1 1 ^ \^U* ' VU* 

P + iV = - J y d yu*^ ~ 2 \ ■ (B.23) 
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The Jarkslog CP Invariant is given by 
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Appendix C 

Verifying M^ 2 in Eq(5.4) 



We derived the Mt2 side of Eq(5.4) by following the analytic solution given by Barr and Lester 
in |170j . In this appendix, we outline how to verify that Mt2 is is indeed given by 

M T2 (x = 0, a, (3, f T = -or - Pt) = 2(«t • Pt + \a T \ \Pt\) (C.l) 

when a 2 = P 2 = and p 2 = q 2 = x 2 = and qt = 0. To do this we note that Mt2 can also 
be defined as the minimum value of (a + p) 2 minimized over p and q subject to the conditions 
p 2 = q 2 = x 2 (on-shell dark-matter particle state), and (a + p) 2 = (/3 + q) 2 (equal on-shell 
parent-particle state), and (a + (3-\-p + q + g)x = (conservation of transverse momentum) [T5] . 



The solution which gives Eq(C.l) has pr = —Pt and qT = —olt with the rapidity of p{q) 
equal to the rapidity of a (/?). We now verify that this solution satisfies all these constraints. 
Transverse momentum conservation is satisfied trivially: {a + P+p + q)x = {a + P — a — P)t = 0. 
The constraint to have the parent particles on-shell can be verified with 2|cty | Iprl — %Pt • O-T = 
2\p T \\q T \ - 2q T ■ (3 T = 2\Pt\\Pt\ + 2a T • Pt- 

Now all that remains is to show that the parent particle's mass is at a minimum with respect 
to ways in which one splits up p and q to satisfy pt + Qt = Pt while satisfying the above 
constraints. We take p and q to be a small deviation from the stated solution pt = ~Pt + $T 
and qT = —olt — &T where 5t is the small deviation in the transverse plane. We keep p and q on 
shell at x = 0. The terms p a , p z , q a , q z are maintained at their minimum by keeping the rapidity 
of p and q equal to a and p. The condition (jp+a) 2 = (q + P) 2 is satisfied for a curve of values for 
5t- The deviation tangent to this curve near | 5t \ = is given by <5r(A) = A z x (otI/^tI + Pt\&t\) 
where x is a cross product, z denotes the beam direction, and we parameterized the magnitude 
by the scalar A. Finally, we can substitute p and q with the deviation <5y(A) back into the 
expression (a + p) 2 and verify that Holt ■ Pt + \®-t\ \Pt\) at A = is indeed the minimum. 
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Appendix D 

Fitting Distributions to Data 



In order to determine Mjy (which is M^° in our case studies) we perform a \ 2 fit between 
ideal constrained mass variable distributions and the HERWIG data or other low-statistics data. 
Because the mass difference is a given, it does not matter whether we work with Mn or My as 
the independent variable. 

First we illustrate the procedure if there is only one distribution being fit as in Chapter |6| 
We work with the constrained mass variable M2C using My as the independent variable. To 
do this we define a x 2 distribution by computing the number of events, Cj, in a given range, j, 
(bin j) of M%c- The variable N is the total number of events in bins that will be fit. Assuming 
a Poisson distribution, we assign an uncertainty, a~j, to each bin j given by 



a 



1 



[N f(M XCj ,My) + C j ). (D.l) 



Here the normalized distribution of ideal events is f{MiCi My), and the second term has been 
added to ensure that the contribution of bins with very few events, where Poisson statistics does 
not appljj^J have a reasonable weighting. Then x 2 1S given by 



The minimum x 2 (My) is our estimate of My. The amount My changes for an increase of x 2 
by one gives our 1 a uncertainty, 5 My, for My (200 . As justification for this we calculate ten 
different seed random numbers to generate ten distinct groups of 250 events. We check that the 
My estimates for the ten sets are distributed with about 2/3 within 5 My of the true My as 
one would expect for 1 a error bars. One might worry that with our definition of x 2 ; the value 
of x 2 P er degree of freedom is less than one. However this is an artifact of the fact that the 
bins with very few or zero events are not adequately described by Poisson statistics and if we 



1 By this we mean that N f(M2Cj,My) has a large percent error when used as a predictor of the number of 
counts Cj when iV jiM^Cj, My) is less than about 5. 
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remove them we do get a reasonable \ 2 P er degree of freedom. The determination of My using 
this reduced set gives similar results. 

In Chapters [7] and [8] we determine the mass with two distributions (the upper and the lower 
bound distributions). We treat these two distributions separately and add the resulting x 2 to 
form a final x 2 ■ In these chapters we choose M-^° as the independent variable. We illustrate 
with Mic,lb and M2c,ub, but the procedure is the same for M%c- 

First we update the definitions for this case. We define Nlb as the number of M 2 c events in 
the region to be fit, and likewise Nub is the number of M^c.xjb events in the region to be fit. The 
M2C of the events are grouped into bins, Cj, in a given range, j. The variable /lb(M2c, M^o) 
is the normalized M2C distribution of ideal events expected in bin j as calculated with an 
assumed M^°, the measured M_, the observed mu distribution, the observed UTM distribution, 
and the appropriate detector simulator. We likewise define the upper bound distribution to 
be fuB(M2c,UB, We also define the background distribution for lower-bound and upper- 

bound distributions to be $b,lb(M2c) and /b,c/b(A^2C,[/b) and the fraction of the total events 
we estimate are from background A. 

Again, we assign an uncertainty, Cj, to each bin j given by 

°LBj( M %) = \ {Nlb ((1 - X)fL B (M 2Cj , M%>) + A/ B)LB (M 2Ci )) + Cj) , (D.3) 

and likewise for the upper-bound distribution. The second term has been added to ensure an 
appropriate weighting of bins with very few events that does not bias the fit towards or away 
from this end-point. In bins with few counts, normal Poisson statistics does not appljj^j 
The x 2 is given by 

Cubj - N UB (1 - A) f UB (M 2 c,u Bj ,M ri ) - N UB A fB,UB{M 2C ,UB j , M^)\ 

a UB,j J 

We calculate ideal distributions for = 50, 60, 70, 80, 90 GeV. We fit quadratic interplant 
through the points. The minimum x 2 (-^x?) °^ the interplant is our estimate of The 
amount changes for an increase in \ 2 by one gives our 1 a uncertainty, <5M^°, for My.o [200 . 

The M3C fits in Chapter [8] were performed following this same procedure. In this case 
we choose = 80,85,90,95,100,105,110 GeV. For the M%c studies we did not create a 




2 By this we mean that N f(M2Cj,M-£°) has a large percent error when used as a predictor of the number of 
counts Cj when iV f(Mzcj, M-^°) is less than about 5. 
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background model and therefore fixed the background parameter to A = 0. The studies involving 
nine different choices of (AMyn, AM^jv) were done using = 90, 100, 110 GeV. 

Some subtleties for which we should check. If the bin size is not large enough, the artifi- 
cially large variations in the distribution sometimes bias the fits to place the endpoint near a 
fluctuation. The bin size should be made large enough so that this does not happen. This can 
be checked for by testing if the results are invariant with respect to changing the bin size by a 
small amount. Because of the larger number of lower bound events, the optimal bin size may 
be different for the upper-bound and the lower-bound distributions. 
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Appendix E 

Uniqueness of Event Reconstruction 



In the Chapters [6] and [8] we claim that the events near an endpoint of M2C and M$c distributions 
(events that nearly saturate the bound) are nearly reconstructed. This appendix offers a proof 
of the claim. To prove uniqueness, we need to establish that as M%c or Mic of an event (lower 
bound or upper bound) approach the endpoint of the distributions, the solutions with different 
values of q and p approach a common solution. 



We begin with M%c- Sec. 8.2 shows that there are at most four solutions given Mat, Mx 



and My formed by the intersection of two ellipses in (p ,Qo) defined by Eqs(8.13 8.14) as shown 



in Fig E.l. Consider the case that an event has a lower bound M^c near My. We are guar- 
anteed that a solution occurs at the true mass scale when we choose the correct combinatoric 
assignments. As one varies the mass scale downward, the two ellipses drift and change shape 
and size so that four solutions become two solutions and eventually, at the value of M%c for the 
event, become one single solution. When the disconnection of the two ellipses occurs near the 
true mass scale, the value of M$c will be near the endpoint. The unique solutions for p and q 
given at M^c are nearly degenerate with the true values of p and q found when one uses the true 
masses to solve for p and q. The closer M$c is to the endpoint the closer the two ellipses are to 
intersecting at a single point when the true masses are used and to giving a unique reconstruc- 



tion. The example pictured in Fig. E.l show an event with M%c within 1% of the endpoint and 
where the p and q are reconstructed to within 4%. This shows that for M%c events that are near 
the endpoint allowing for any choice of combinatorics then nearly reconstruct the true values for 
p and q. If there are combinatoric ambiguities, one need to test all combinatoric possibilities. 
If the minimum combinatoric option has a lower bound at the end-point, the above arguments 
follow unchanged. The above arguments can be repeated to show M 3 c,ub near the endpoint 
also reconstructs the correct p and q. 

Next we turn to Mic- For every event the lower-bounds satisfy M2c(My — M^r) < M^ciM-Y — 
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Figure E.l: Shows the ellipses defined for p Q and q Q in Eqs(8.13 8.14 ) using the correct mass scale 
for an event that nearly saturates the M^c endpoint. For this event, the M%c hes within 1% of 
the endpoint and reconstructs p and q to within 4%. Perfect error resolution and combinatorics 
are assumed. 

Mn,Mx — Mn). With M2C the propagator (p + c^) 2 , which we can equate with x\i 1S no * 
fixed. The kinematically allowed values for xx are M'L < x\ ^ assuming the visible 



states «i and 0,2 are massless. Eq(8.12) shows that p and q solutions are linear in x\ with 
no terms dependent on xx alone or other powers of xx- Including x\ as a ^ ree parameter 
in Eqs( 8.13|8.14 ) leads to two ellipsoids (or hyperboloids) in the space {p ,Qo,Xx)- We will 
assume without loss of generality that these are ellipsoids. Again, at the true mass scale we 
are guaranteed the two ellipsoids intersect at an ellipse. Now as one varies the mass scale the 
two ellipsoids drift and change shape and size. The M2C value then corresponds to the mass 
scale where the two ellipsoids are in contact at one point. As we select events with a value of 
M2C that approaches the true mass scale the intersection of the two ellipsoids shrink to a point 
giving a unique reconstruction of p and q. The combinatoric ambiguities for M2C are avoided 
by selecting events with two distinct OSSF pairs. Events that saturate the upper bound of M2C 
also reconstruct p and q by the same logic as above. 
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Appendix F 

Acronyms List 



• ACDM A(Cosmological Constant) Cold Dark Matter 

• CKM Cabibbo-Kobayashi-Maskawa 

• CMB Cosmic Microwave Background 

• DR Dimensional Reduction (Renormalization Scheme) 

• GIM Glashow-Iliopoulos-Maiani 

• GJ Georgi-Jarlskog 

• GUT Grand Unified Theory 

• ISR Initial State Radiation 

• KK Kaluza Klein 

• LHC Large Hadron Collider 

• LKP Lightest Kaluza Klein particle 

• LSP Lightest supersymmetric particle 

• MS Minimal Subtraction (Renormalization Scheme) 

• MACHO MAssive Compact Halo Objects 

• OSSF Opposite-Sign Same Flavor 

• QFT Quantum Field Theory 

• REWSB Radiative Electroweak Symmetry Breaking 

• RG Renormalization Group 

• RGE Renormalization Group Equations 

• SM Standard Model 

• SSB Spontaneous Symmetry Breaking 

• SUSY Supersymmetry 

• UED Universal Extra Dimensions 

• UTM Upstream Transverse Momentum 

• VEV Vacuum Expectation Value 
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